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During the past decade the interaction of light with multi-atom ensembles has attracted a lot 
of attention as a basic building block for quantum information processing and quantum state 
engineering. The field started with the realization that optically thick free space ensembles can be 
efficiently interfaced with quantum optical fields. By now the atomic ensemble - light interfaces 
have become a powerful alternative to the cavity-enhanced interaction of light with single atoms. 
We discuss various mechanisms used for the quantum interface, including quantum nondemolition 
or Faraday interaction, quantum measurement and feedback, Raman interaction, photon echoe 
and electromagnetically induced transparency. The paper provides a common theoretical frame for 
these processes, describes basic experimental techniques and media used for quantum interfaces, 
and reviews several key experiments on quantum memory for light, quantum entanglement between 
atomic ensembles and light, and quantum teleportation with atomic ensembles. We discuss the two 
types of quantum measurements which are most important for the interface: homodyne detection 
and photon counting. The paper concludes with an outlook on the future of atomic ensembles 
as an enabling technology in quantum information processing. Published in Reviews of Modern 

Physics 82, 1041 (2010) 

|http : //link . aps . org/doi/10 . 1103 /RevModPh ys . 82 . 1041| 
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I. INTRODUCTION 

A. History and motivation 

Quantum features of atom-light interaction have been 
among the central issues in physics since the early days 
of quantum mechanics. Starting in the 1960s, with the 
development of quantum optics - the field where sec- 
ond quantization of light is central - quantum electro- 
dynamics became part of optical and atomic physics. 
For decades after that the inherently quantum features 
of atom-light interaction have been studied primarily 
within the framework of cavity Quantum Electrodynam- 
ics (QED) where light can be efficiently coupled to a few 
atoms or even to a single atom. Despite the spectac- 
ular progress achieved in this direction, the complexity 
and technical challenges associated with an atom strongly 
coupled to a high-finesse cavity were calling for alterna- 
tive approaches. 

A new approach to the matter-light quantum interface 
came with the realization of the fact that a large collec- 



tion of atoms - an atomic ensemble - can be efficiently 
coupled to quantum light if a collective superposition 
state of many atoms can be utilized for the coupling. 
The simplest example of such coupling is presented in 
Fig. [TJi. A collection of atoms in the ground state is illu- 
mi nated with two modes of light a + and a_ . As shown 
by |Kuzmich et al.\ (Tl997), if the modes possess quan- 
tum correlations (entanglement) and light is absorbed 
by the atomic ensemble, the quantum correlations can 
be mapped on the collective superposition of the two fi- 
nal states of the atoms. The strong coupling condition 
in this case amounts to the requirement of large resonant 
optical depth d of the atomic ensemble. It later turned 
out, that the requirement d ^> 1 is the most significant 
requirement for all types of the quantum interface be- 
tween atomic ensembles and light known up to now. The 
experiment demonstrating that a quantum feature of ra- 
diation (squeezing) can be transferred onto atoms via the 
process shown in Fig. |TJl was performed by Hald et ( 
(1999). This appro ach has been fu rther developed using 
photon echo ideas ( Moiseev[ 2QQ3| ). 

A natural next step was to utilize long lived atomic 
ground states for the interfac e via a Raman interaction 
in a A-scheme, as proposed by Kozhekin et al.\ ( |2QQQ[ ) for 
storage of squeezed states. The Raman process together 
with Electromagnetically Induced Transparency (EIT) 



( [Boiler et a/.[|1991[|Fleischhauer et al.||2005||Lukin[|2003 ) 
have soon become important routines for quantum inter- 



faces. After Hau et al. (1999) demonstrated that EIT al- 



lows for very slow propagation of light through an atomic 
ensemble it was quickly realized that reducing the group 
velocity to zero would enable an atomic memory for light 
( Fleischhauer and Lukin 2000 Lukin et a/. , [2000) and 
the first experimental demonst rations of this for classical 
pulses have been presented by Liu et al. (2001); Phillips 
et al.\ (12001). 



Quantum nondemolition (QND) measurement (Bra- 
ginsky and Khalili, 1996) based on light-matter interac- 



tion has emerged as a powerful tool for quantum state en- 
gineering, first in the cavity QED setting and then in the 
atomic ensemble context as an efficient method for gen- 



eration of spin squeezing (Kuzmich et al., 1998). Shortly 



thereafter QND interaction with atomic ensembles has 
become one of the main instruments for the quantum 
interface. 

The process depicted in Fig. [T^i is a rudimentary ex- 
ample of one of the main routines for atoms-light quan- 
tum interface: the quantum state transfer from light to 
atoms, or the quantum memory for light. The ability 
for mapping, storing, and retrieving quantum states of 
light - the natural long-distance carrier of information 
- onto the material storage medium is one of the major 
enabling procedures in quantum information processing. 
In this review we cover various approaches to the quan- 
tum memory, including the Raman process, EIT, photon 
echo, and the QND measurement and feedback. We re- 
view methods which provide a long term quantum mem- 
ory with the fidelity better than any classical procedure 



can achieve, as well as approaches which allow to preserve 
entanglement in the process of storage and retrieval. The 
interface can be implemented either via interaction only 
or by the teleportation-like procedure involving gener- 
ation of entanglement, Bell measurement on light and 
quantum feedback onto atoms. In this article we will 
discuss both approaches. 

The second most important routine for the quantum 
interface is generation of entanglement between light and 
atoms. The light /atoms entanglement in turn enables 
generation of entanglement between remote atomic en- 
sembles, as well as atomic teleportation and entangle- 
ment swapping protocols. Furthermore the light-atom 
entanglement also allows for quantum memory through 
light-atom teleportation. 

The quantum interface can be formulated either in 
the Schrodinger or in the Heisenberg picture. For ex- 
ample, the transformation of a quantum state of light 
into a quantum state of atoms in the Schrodinger picture 
U\$l l^a) —> |0l^a) corresponds to the operator trans- 
formation JJ^ajJJ = aL in the Heisenberg picture. The 
two pictures are equivalent, and we will mostly use the 
Heisenberg picture throughout this review. 

In this review we discuss protocols based on both ho- 
modyning of light and photon counting. The most dra- 
matic difference between the two approaches is that a 
single homodyne measurement does not necessarily dis- 
tinguish between a vacuum and a non-vacuum state, 
whereas an ideal photon counter does. This makes ideal 
photon counting insensitive to losses if the protocol is 
conditioned on a click of the detector. This feature is 
important as an elementary purification mechanism, but 
it also makes protocols which use it probabilistic. On 
the other hand, homodyning always yields a measure- 
ment result and is thus deterministic. Another difference 
between the two measurements is that photon count- 
ing yields a discrete variable result, whereas homodyning 
yields a continuous variable outcome. From a practical 
perspective detectors used for homodyning are almost 
perfect in their quantum efficiency and dark current, 
whereas photon counters usually are less than perfect 
(although the progress in their development driven by 
quantum information applications is remarkable). The 
distinction between the two approaches is, however, not 
strict, as discussed later in this review. For example, 
photon counting can be used as a deterministic charac- 
terization of a protocol if the absence of a photon count 
(the vacuum contribution) is included in the analysis. 
Continuous variable outcome of a homodyne measure- 
ment can be "digitized" if suitable superposition states 
are employed. Various figures of merit are used for char- 
acterization of quantum interfaces, as discussed in Sec. 

EAl 

Probabilistic protocols based on generation of a single 
collective atomic excitation of an atomic ensemble follow- 
ing detection of a photon emitted by the ens emble have 
been actively developed in the recent years (|Chaneliere 
TFal] [20071 IChen et al.\ [20061 |Chou et al.\ |200"4| |2007| 



Kuzmich et al. , 2003 ; Matsukevich et al. , 2006 ; Matsuke- interaction for reasons discussed below. 



vich and Kuzmich[ |2004[ |van der Wal et a/J[2003] ). This 



approach has been motivated by a proposal for a quan- 
tum repeater with atomic ensembles ([Duan et al. 2001). 



The research on quantum repeaters deserves a separate 
review paper and will only be rather briefly discussed in 
this article. 

The requirements for the atomic memory may differ 
depending on the particular application. A distributed 
quantum computer network requires the complete set of 
memory capabilities: mapping of the light state onto 
memory, storage and operations on the memory state, 
and retrieval of the memory state back onto light for fur- 
ther processing. Applications in quantum communica- 
tions, which involve local operations on stored states and 
classical communications between partners, often only re- 
quire a measurement of the memory state in a specific 
basis, i.e., no full retrieval of the quantum state of the 
memory back onto light is necessary. Yet other propos- 
als, such as linear optics quantum computing which uses 
off-line entanglement resources require only the retrieval 
of the atomic state onto light but it has to be rather 
efficient and with high fidelity ( Menicucci et al. t 2006 ) . 




FIG. 1 Elementary level schemes: (a) A simple absorption 
scheme with the quantum fields a± mapped onto the atomic 
states a a, (b) Beam splitter type interaction - basis for Ra- 
man and EIT memory schemes, (c) Parametric gain type in- 
teraction - basis for entanglement schemes, (d) Quantum non- 
demolition (or Faraday) interaction - basis for entanglement, 
memory and teleportation schemes 



B. Elementary level schemes 

Naturally, the atomic levels used for storage of a quan- 
tum state should be long lived, with particular require- 
ments for the lifetime depending on applications. For 
example, for a memory used in a long distance commu- 
nication protocol, the memory lifetime usually should be 
longer than the time required for classical communica- 
tion over this distance. For a few hundred kilometers this 
time is of the order of 10 -3 sec. Short distance applica- 
tions may require shorter memory time but one should 
keep in mind that a low loss fiber loop can be a strong 
competitor for a short term atomic quantum memory for 
light (Pi ttman and Franson| [2002). A 5 km fiber loop 
can, in principle, store a photon for 25 /isec with only 
20% losses at the telecom wavelength. However even for 
short storage times atoms have the important advantage 
of being able to provide on-demand retrieval and may in 
addition be advantageous if a nontrivial operation has to 
be performed on the stored quantum states. 

Preferably the optical atomic transitions used for cou- 
pling light to the storage ground states of the atoms 
should be strong in order to have a large bandwidth of 
the memory. Therefore strong dipole allowed transitions 
are typically used for the interaction, but some experi- 
ments, in particular in solid state systems, compensate 
for weak optical transition by having a large number of 
atoms. Figs, [lib), (c), and (d) present the atomic level 
schemes typically used for the interface. The (b) and (c) 
parts of the figure present the A - scheme used in the Ra- 
man and EIT memory schemes as well for entanglement 
generation. The (d) part shows the so-called Faraday 
interaction which is sometimes also referred to as QND 



In (b) the atoms are prepared in the ground state 
coupled to the quantum field cll (thin line), with the 
bold line showing a strong coupling field. With a large 
detuning from the optically excited states, the interac- 
tion Hamiltonian for such a system, after adiabatic elim- 
ination of the excited state, can be cast in the form 
H = Xbs^l^a + H.C. In quantum optics this Hamilto- 
nian is often referred to as the beam-splitter Hamiltonian 

xmhardt (2003). The interface mixes the input atom 
and light states as a "beamsplitter" and the "reflection 
coefficient" of unity corresponds to a perfect state swap- 
ping between light and atoms. The detailed derivation 
of this Hamiltonian for the light-atoms interaction will 
be given later in the article, but the intuitive picture is 
obvious - if a single photonic excitation is removed (an- 
nihilated) from the field cll, a single collective atomic 
excitation a\ is created. If this process is efficient, it 
works as the Raman-type quantum memory for light in- 
troduced for atomic ensembles by Kozhekin et al. (2000) 
and described in Sec. IV. Dl The same level scheme can 
be used for the EIT-based memory w here the fields are 
resonant with the optical transition (Fleischhauer and 



Lukin 2000 



Lukin et al., 2000), although in this limit 



it is essential to account for spontaneous emission and 
the effective beam splitter Hamiltonian is less applicable 
(see the detailed theory in Sees. II). EIT based quantum 



memory experiments are also described in Sec. V.D 



Part (c) of the figure shows the same atomic structure 
but now the fields are arranged in a way which can be 
used for the atoms-light entangling interaction with the 
Hamiltonian H = xp^l&a +H.C. The Hamiltonian is for- 
mally identical to the parametric gain interaction Hamil- 
tonian |Leonhardt] (|2003[) which has been a workhorse for 



studies of entangled and non-classical states of light since 
the 1960s. The important new feature in the present case 
is that two entangled operators belong to a light mode 
and an atomic mode respectively. This kind of entan- 
glement has been used for unconditional light-to-atoms 



teleportation experiment described in Sec. [VI] and its 
probabilistic version discussed in Sec. |IV.C| is the basis 
for the repeater protocol of Duan et al. (2001). 



To complete the discussion of elementary level schemes 
used for basic interface routines we consider the four-level 
scheme shown in figure Fig. Illd). The Hamiltonian for 
this interaction can be obtained by combining the beam- 
splitter Hamiltonian (Fig. [I] (b)) and the parametric 
entangling Hamiltonian with equal coupling constants, 
Xbs = Xp (Fig. [l] (c)) provided that the two quantum 
fields (thin lines) belong to the same mode: H ~ \PlPa 
where the canonical operators for light and atoms which 
obey the canonical commutation relation [X, P] = i have 
been introduced. This interaction allows for a QND mea- 
surement of the atomic operator Pa by means of de- 
tection of the light operator Xl,. As discussed in de- 



tail in Sections IV A,B the QND measurement projects 
atoms into an entangled state. The same interaction is 
often called the quantum Faraday interaction because in 
case of magnetic levels it leads to polarization rotation of 
light. The QND-Faraday interaction of atoms and light 
followed by the measurement on the light and the feed- 
back conditioned on the measurement applied onto atoms 
was used to demonstrate quantum memory for light as 
described in Section [Vj 

As shown theoretically and experimentally in 



(Wasilewski et al., 2009), a more general combination 



of the beam-splitter Hamiltonian and the parametric 
entangling Hamiltonian with unequal weights performs 
both those operations simultaneously. 

To summarize, the basic features of most quantum in- 
terface protocols to date can be understood by analyz- 
ing simple 3- or 4-level atoms. Besides the condition 
d ^> 1 mentioned above, another unifying feature for all 
approaches which use multi-atom ensembles is the pos- 
sibility to initialize the ensemble, e.g., by optical pump- 
ing, in one of the ground substates. Choosing suitable 
atomic transitions, and polarizations and frequencies of 
the quantum and classical coupling fields one can choose 
between various routines, such as memory, entanglement 
and Faraday interaction, as shown in Fig. [I] In the fol- 
lowing theoretical description we will provide a unified 
approach to all these types of interfaces. 

In the literature the interface protocols are often di- 
vided into those for states of continuous variables and 
those for discrete variables, or qubits. The former are 
usually based on the Faraday interaction and are de- 
scribed most conveniently in terms of X, P operators 
measured by homodyne detection, while the latter are 
commonly based on the A-type interactions in combina- 
tion with counting of single photons and are most easily 
described in terms of a, a) representations. One of the 
goals of this review is to show that continuous and dis- 



crete variable protocols can in many aspects be treated 
on equal footing, and that the choice of variables is de- 
fined by the convenience of description and the type of 
measurements involved. E.g., in the ideal limit a memory 
protocol which is most conveniently described by X, P 
operators could be used to store a single photon with 
perfect fidelity. On the other hand, the state of the 
atomic memory in protocols which use a, a) represen- 
tation can be conveniently analyzed by atomic tomog- 



raphy in the X, P basis ( Fernholz et al. 



2008 Sherson 



et al. 2006b) (for a recent review on quantum tomogra- 
phy see ( |Lvovsky and Raymer 2009). The type of er- 
rors which appear under non-ideal conditions is of course 
different for different protocols, but irrespectively of the 
specific protocol, the condition d ^> 1, leads to fewer er- 
rors. Which protocol to use, for a given optical density, 
therefore depends on the specific application and a de- 
tailed analysis of the imperfections should be made in 
each situation. 

This review strives to provide a coherent picture of the 
work on the quantum interface between light and atomic 
ensembles using various approaches, atomic media, and 
protocols. It includes the discussion of major experimen- 
tal achievements to date and concludes with the analysis 
of the current limitations and future goals. 



II. THEORETICAL BACKGROUND 

A. Description of light and atoms 

Harmonic oscillators Throughout this review we shall be 
dealing with single modes of the electromagnetic field and 
collective spin excitations of atomic ensembles which can 
be well approximated by harmonic oscillators with canon- 
ical position and momentum operators X n and P n , where 
n refers to the mode number. In most cases we shall 
omit the hats on the operators in what follows. These 
canonical operators are dimensionless with the standard 
commutator 



[^■n^MnJ — IVmn- 



(i) 



The harmonic oscillators can also be described in terms 
of the annihilation operators 



dn = -i=(X n +iP n ), 



(2) 



which have commutation relation \a ni a, 



ml ~ u nm- 



Instead of labeling by a discrete number n the modes 
can be denoted by a continuous parameter, e.g., the po- 
sition vector r with the commutation relations 



[X(r),P(f')} = i6(r-r'), 
[a(f),a\f')] = 6(f-r'). 



(3) 
(4) 



In some cases we deal with the storage or transfer of a 
set of n modes. Such discrete modes can be constructed 



from the continuous modes by introducing a complete 
orthogonal set of mode functions {u m (f)} satisfying 



/ dr u* m {r)u n {f) = 6 mn 



(5) 
(6) 



If we now define the discrete annihilation operators 

a<m= dr u* m {r)a{r) (7) 

they have the appropriate commutator [an,^] = 5 mn . 



Light Light beams travelling in the ^-direction can be 
described (in c.g.s. units) in the paraxial approximation 
by a quantized electric field 



E(f) 



m,a,k 



^^e a u m (f ± ]z)e ikz a Ljmak ^R.C. (8) 



Throughout this review we set ft — 1. /is the length 
of the quantization volume, and the sum is over the 
polarization a, the transverse mode number ra, as well 
as the longitudinal wave vector k. The mode functions 
Um(r_\_;z)i where f± = (x,y) describes the transverse 
profile of the beam, form a complete orthogonal set in 
the plane transverse to the propagation direction 



/■ 



d fj_u^(f ± ;z)w m /(f ± ;z) = 8 n 



(9) 



In the above we have assumed that the fields belong to 
a narrow frequency band such that for all modes the fre- 
quency under the square root is ujq. Secondly, the field 
should in general be expanded into a complete set of 
modes Uk(r), but in the paraxial approximation, we have 
assumed that we can factor out a polarization vector e a 
as well as ignore the k dependence of the transverse mode 
function u rn (f±; z). 

Instead of using longitudinal wave vectors we use a 
slowly varying position space annihilation operator de- 
fined by 



a L 



>™W = \ 7 J2 eiik ~ k0)Z+iUJOta ^™*k, (10) 



where c is the speed of light and fco = ooo/c. In the con- 
tinuum limit / — > oo this operator has the commutation 
relation 

[«L,mcr(^),ttl / , m / a /(^ / )] = Cfim,m>5*,*'5(z- z). (11) 

Note that we have chosen here a normalization with 
c appearing in the commutator. With this normaliza- 
tion (i) the travelling fields aL jT na( z ^t) — aL,mcr(z — 
ct) considered below have the commutation relation 



appropriate for operators which are a function of t 

(kmafe^at.m'a'fe*')] = ^m,m^a,a^(t- 1')), and (h) 

with this normalization aL, m a(^)^L,ma(^) describes the 
flux of photons in mode ra with polarization a at position 
z. 

In terms of this operator the electric field is given by 



27TCo>o 



E{f) = y^^e <TUm (f ± ;z)e i ( fc ^-^ t )a L , mff (z)+H.C. 

m,(j 

(12) 
Below we mainly deal with a single transverse field mode 
and a single polarization, so that the sum in the expres- 
sion above can be omitted. For an introductory textbook 



to continuous mode quantum optics we refer to ( [Loudon 
2004b. 



We first discuss the theory for atoms with t wo 

10} 



Atoms 

stable ground states 



II.D 



we 



and |1). In subsec. 
show how one can in many cases reduce the description 
for multilevel atoms to two state atoms. The two ground 
states are conveniently described in terms of angular mo- 
mentum operators. We shall here use the x as the quan 



tization axis for consistency with I Julsgaard et al. (2001 



|2004a[ ) and |Sherson et al] ( |2006c| ) . The angular momen- 
tum operators describing the rath atom are 



3x,m = 2(|0U0|-|1U1|), (13) 

j+,m = jy,m + ijz 1 m = \0)m(Mi ( i4 ) 



where j+, m is the operator which raises j Xjin by unity. 

We shall be interested here in collective variables for 
an ensemble containing many atoms. In the simplest case 
such collective operators are given by the total angular 
momentum operators J\ = ^ m ji,m (with / = x,y,z), 
which fulfil the standard angular momentum commuta- 
tion relation 



[Jy, J z\ — iJx- 



(15) 



The collective state with all atoms in state |0), then cor- 
responds to the state \J = A^^/2, M x = Na/2) with total 
angular momentum quantum number J = Na/2 and an 
eigenvalue of J x equal to iV^/2, where Na is the number 
of atoms. If we consider a large number of atoms and 
only weakly perturb the system (only change the state 
of a few atoms), we can approximate the J x operator by 
its expectation value J x ~ (J x )- For readers who feel 
uneasy about replacing an operator by its mean value, a 
more rigorous formulation can be made b y using the so 
called Holstein-Primakoff transformation (Holstein and 



Primakoff 



1940 



Kittel 



1987) or in terms of a Wigner 



group contraction (|Arecchi et al> 1972). Without loss 



of generality we can assume the expectation value (J x ) 
to be positive and we can then introduce new canonical 
position and momentum operators by 



X A 



Jy 



VW 



Pa 



J z 



vW 



(16) 



From (|15| we immediately see that these operators sat- 
isfy the standard commutation relation for position and 
momentum (HI). The collective annihilation operator is 
then 



a A 



^ + tfA.E ro J+,m Em|0)m(l| 



V2 



V^W) v^W) 



(17) 



To get a feeling for this operator, let us consider the 
action of the creation operator a/ A . If we apply this oper- 
ator to the initial state, where all atoms are in state |0), 
we create a symmetric superposition of one atom being 
flipped 



:^|0,0,....,0,l m ,0,...,0). 



^|0,0,0,...,0) 



(18) 
Here |0, 0, ...., 0, l m , 0, ..., 0) is the state, where all atoms 
except the rath atom are in state |0). 

The collective operators introduced above are conve- 
nient for describing the entire ensemble. We shall, how- 
ever, also be dealing with situations where we need to 
consider collective operators which do not involve all 



atoms with an equal weight (Kuzmich and Kennedy 
2004). In the literature such situations are often de- 
scribed by dividing the ensemble into small box es and 
constructing collective operators for each box (Fleis- 
chha uer and Richter[ |1995| |Raymer and Mosto wski , 



1981) 

Here we shall use a slightly different formalism 



(S0rensen and S0rensen l 2008). For a collection of atoms 



at positions r\ , ?2 , 
tion function 



? tn a we define the density distribu- 



»(*=) = E*( f 



»)• 



(19) 



In condensed matter physics this density distribution 
function may be used to describe s cattering from struc- 
tures ( |Chaikin and Lubensky , 1995): averaging this den- 
sity distribution over the random positions of the atoms 
gives the average number density of the atoms n(f) = 
(n(r)), whereas higher order correlations like (n(f)n(r f )) 
describe the correlation responsible for Bragg scattering 
(in this context the classical functi on ([l9|) is sometimes 
refereed to as the density operator ( |Chaikin and Luben- 
|sky[|!995| ), but to avoid confusion with quantum mechan- 
ical operators we shall avoid this terminology). Similarly 
we may introduce continuous atomic spin operators by 



fk(r) = Yl S ^~ ^m)jk,m, 



3k 



(20) 



where k = x,y,z,+,—. A position dependent atomic 
annihilation operator can then be introduced by 



a>A(r) 



1 



VW*W) 



3+(f), 



(21) 



where the average spin density (j x (r)), which we assume 
to be positive, is the quantum mechanical expectation 



value with respect to the internal state averaged over the 
random (classical) position of the atoms. The annihila- 
tion operator has the commutation relation 



[a A (r),a\(r , )}=5(r-r , )i L 



3x(r) 



0*(*)> 



5(r-r'), (22) 



where we in the last step have assumed that the fluctua- 
tions in the mean spin are much smaller than its average. 
We shall use this approximation throughout this article. 



The operators (21 ) will be convenient for describing the 



spatial dependence of various operators. To relate them 
to single mode operators such as Eq. (17), we can intro- 



duce a normalized set of mode functions u n (f) fulfilling 
the orthogonality and completeness relations in Eq. (|6|. 
We can then construct single mode operators as in Eq. 
^\. In particular let us consider the normalized mode 
w sym (f) = \f{jx(f)) / {J x ). If we use this mode to con- 
struct a collective operator, we see from Eq. (|2lt that 



this produces the symmetric operator defined in Eq. ( 17 ) . 
The j x (r) operator will appear in the Hamiltonians 
below, although sometimes it is more convenient to have 
expressions which only involve the annihilation operator 
dA(f)' To find the equations of motion we need to take 
the commutator of a^(r) with the Hamiltonian, but from 
the commutation relation 



[a>A(r)J x (r')] 



S(f— f') 
S(f— r')aA{f) 



^2S(r-r m )j +Ji 



(23) 



we see that we get the same result if we make the re- 
placement 



3x(r) -+ --a\{r)a A {r) 



(24) 



and use the commutation relation in Eq. ( |22| ) (the first 
term is included to ensure that j x has the right value 
in the vacuum state of a a, where all atoms are in state 
|0)). This replacement holds even as an exact opera- 
tor identity in the framework of the Holstein-Primakoff 
transformation (Holstein and Primakoff, 1940). 



B. Interaction of light with model atoms 

We consider atoms with stable ground states denoted 
by \g m ) and excited states denoted by |e m ). The Hami- 
tonian H = H^ + Ha + ^int describing this system is 
the sum of the field energy H^, the atomic energy Ha, 
and the interaction Hamiltonian Hi nt . Below we first con- 
sider the atomic and interaction parts of the Hamiltonian 
and derive an effective interaction Hamiltonian involving 
only the ground states of the atoms. We then include 
the Hamiltonian H^ responsible for the propagation of 
the light field, and derive coupled equations of motion 
for the light and atomic operators. 



Interaction with a single atom To describe the atomic part 
and the interaction let us first consider only a single atom 
at location r. In a rotating frame with respect to the 
laser frequency the atomic Hamiltonian is given by Ha = 
X) m A m |e m )(e m |, where A m is the detuning of the mth 
excited state with respect to the laser frequency. 

In the dipole approximation the interaction between 
light and atoms is described by the Hamiltonian iJ int = 
—E-D, where D is the electric dipole operator for the 
atom. In Appendix |IX. A we describe the adiabatic elimi- 
nation of the excited state and derive an effective ground 
state Hamiltonian 



H Lt 



^2 V m',m{r)\9m'){9r) 



(25) 



where the coupling matrix V m i m is given by 



Vrr. 



.(f) 



E 



(E(-\r) 



D {+ ) ,MD { -), 

m,m" I \ m' ,\ 



(26) 
Here superscripts (+) and (— ) refer to the positive and 
negative frequency components of the electric field and 
dipole operators (the positive frequency part, is the 
part of the operators which removes an excitation, e.g. 
D^] m = {9m>\D\e m )) JLoudon||2004[. (Note that H' h 
is not the same as H nt 



-E • A since H[ nt 



int 

contains a 



contribution from Ha as discussed in Appendix IX. A ) 



Interaction with many atoms Let us now consider the sit- 
uation, where we have many atoms. We are mainly in- 
terested in describing the interaction of a weak quantum 
field with an ensemble driven by a strong classical field. 
In this subsection we only describe the interaction be- 
tween the atoms and the forward propagating quantum 
fields and ignore the spontaneous emission due the cou- 
pling to all other e.-m. modes. We shall include the 
decoherence caused by spontaneous emission in subsec. 
II. C| To simplify the theory we shall only derive equa- 
tions of motion for the initially fully polarized ensem- 
ble (J x (r,t)} ~ n(r)/2 (we ignore the difference between 
the density distribution n(f) and its average value n{r)). 
The atomic operators defined in Eq. ( [21] ) are, however, 
well behaved annihilation operators even for an ensemble 
which is not fully polarized provided that the fluctuations 
of the mean spin are small. This will for instance be the 
case if an ensemble containing many atoms is prepared 
with imperfect optical pumping. 

The interaction Hamiltonian can be obtained by sum- 
ming the single atom Hamiltonian (25) over all atoms. 



This sum may be replaced by an integral by introducing 
the continuous atomic annihilation operator defined in 



Eq. (21) 



H = I d 3 - 



Jd 3 r (v^)a A (f0^oi(rO+H.C.) 

+ n(r)V o + a\(r)a A (r)(V u - V 00 ) 



(27) 



where we have used the replacement (24) for the spin 
operator j x (r)- In Appendix |IX.B we use this general 
Hamiltonian to derive Hamiltonians for the three differ- 
ent model systems in Fig. [l](b)-(d). 

For some situations the 3D Hamiltonians derived in 
Appendix |IX.B can be reduced to one dimension. Let us 
perform this reduction for the beam splitter interaction 
Hamiltonian in Eq. ( [95] ) ? corresponding to the level con- 
figuration as shown in Fig.[2Va), where both classical and 
quantum field is traveling in the z direction. Lets us now 
assume that the density is independent of the transverse 
coordinate n(f) = n(z), and that the classical laser field 
is also constant transverse to the propagation direction. 
Since the mode functions u m (f±; z) form a complete set 
in the plane, we can expand the atomic operator a A (f) 
in the same set 

a ^) = Yl u ™{r±\ z)a Ayrn (z), (28) 



o>A,m(z) = / d 2 ?± u* m (r ± ;z)a A (r). (29) 



where 



These new operators will then have the appropriate com- 
mutation relation [a A ^ rn (z),a ( Arn/ (z f )] = S(z - z')5 m ^. 



We insert this expansion into the Hamiltonian (95), and 
then use the orthogonality relation ^ integrate over 
the transverse coordinate to obtain the one dimensional 
Hamiltonian 



Hbs 



= I dz 



■\V(z,t)\ 2 
4A 



J2 a Am( Z ) a A,m(z) 



1*7(2)1 



A J2 a L,m( z ) a L,m(z) 

m 



2A 



(30) 



where the the coupling constant g(z) and slowly varying 
resonant Rabi frequency Q(z,t) are given by 



9{z) 



27rujn(z) 



Do, 



(31) 



Q(z,i) = 2D { el } • (^ (+) ) exp(-i(k z - u Q t)) 



and the dipole element for the |0)-|e) transition with 
the polarization of the quantum field e q is given by 



Do — ^e~o ' ^q- The Hamiltonian above consists of three 
terms: the first line is the AC- Stark shift of the atomic 
ground state, the second is the index of refraction of the 
gas, and the last line, which is the most important for our 
discussion here, describes the exchange of excitations be- 
tween atoms and light. 



\e) 
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FIG. 2 A configuration with dominant population of level |0): 
(a) Beam splitter interaction with the quantum field of the 
single photon Rabi frequency g on the |0) — ► |e) and the strong 
light with the Rabi frequency Q on the |e) — > |1) transition. 
The two fields are in two-photon resonance with a detuning A 
from the excited state |e). Decay due to spontaneous emission 
goes back to one of the ground states at the rates 70(1) or to 
other levels (not shown) summing up to a total rate 7. (b) 
Parametric gain interaction, (c) Faraday interaction: Atoms 
are polarized to |0), the strong field is linearly polarized along 
x and drives the up-transitions with the Rabi frequency Q 
and the quantum field in y polarization couples to the cross- 
transitions with the single photon Rabi frequency g. The 
fields are in two photon resonance with a detuning A from 
the excited states. Decay due to spontaneous emission goes 
back to one of the ground states at the rates J x ( y ) or to other 
levels (not shown) at the rate 7. 

Note that each of the transverse modes of the light 
field in (30) talks to a single transverse mode of the 



atoms, which then couples back to the same transverse 
light mode. Since the dynamics is actually the same for 
all the involved transverse modes, the atomic ensembles 
may in fact be used as a memory for for multiple trans- 
verse modes (Camacho et al\ |2007| |Shuker et al\ [2008 



Vasilyev et al.[|2008| |Vudyasetu et al.||2008| ). A similar 



description of the reduction from three to one dimension 
is also presented in Andre ( 2005 ) . 



We have used here the same central frequency ujq for 
both the quantum and classical fields, which is true for 
degenerate |0) and |1) states. If they are not degenerate 
the energy difference between the two states can be ac- 
counted for by changing the frequency uj' of the classical 
field. In this case, however, an additional phase factor 
exp(z(/c — ko)z) associated with the difference of the k 
vectors appears. We will discuss the effect of this phase 
factor in Sec. III. Gl and Sec. IVIII 



For the parametric gain Hamiltonian (96), cf. 
Fig. |2|b) , the reduction from three to one dimension can 
be achieved using the same procedure as above. The 



only difference is that instead of Eq. ( 29 ) we should now 
define the discrete atomic operator by 



/ 



aA,m(z) = / d 2 f ± u m (f ± ;z)a A (f). 



(32) 



Although the omission of the complex conjugate in this 
expression compared to Eq. (29) may seem of minor 



importance, it actually does play a role for several exper- 
iments as we will discuss in Sec. IV. C Omitting the sum 
over multiple modes, we arrive at the Hamiltonian 



H G = Jdz B?b$ta\(z)a A {z) 



9*(zMz,t) 
2A 



a[(z)a A (z)+li.C. 



(33) 



Finally we derive the one dimensional Hamiltonian for 
the QND (Faraday) interaction as shown in Fig. [2|c), 
where the x-polarized classical field couples the vertical 
transitions and a quantum field in y polarization couples 
diagonal transitions. From the figure we see that the 
Faraday interaction is essentially a combination of the 
beam splitter and the parametric gain interaction taken 
with the same coupling strength. In fact, as discussed 
in Appendix |IX.B| the three dimensional Hamiltonian for 
the Faraday interaction is simply Hp = (Hbs — Hg)/V2- 
To reduce the problem to one dimension we would like 
to introduce atomic operators similar to Eqs. ( p9[ ) and 
( 32 ) , but now there is an ambiguity as to which of the two 



forms one should use. To avoid this ambiguity we assume 
that the mode fu nctions u m (f± ; z) are real [ for Hermite- 
Gaussian modes (Milonni and Eberly 1988), this condi- 



tion can only be satisfied if the Fresnel number is much 
greater than unity F = Wq/XL ^> 1, where wq is the 
beam waist, L is the length of the medium, and A is the 
wavelength of light (Miil ler et al\ 2005| |S0rensen and 
|S0rensen] |2008| )], and that g*(r)Q(r,t) has a constant 
phase, which we take to be zero. With these assump- 
tions the Hamiltonian can be reduced to the simple form 



H F = - fdz 



9*(z)n(z,t) 



Pl(z)pa(z), 



(34) 



We have omitted here the index of refraction of the gas 
for the reasons discussed in Appendix |IX.B| 

We emphasize that the above 3D to ID reduction pro- 
vides a highly simplified treatment of the propagation of 
light through an atomic gas. In particular a more gen- 
eral treatment should include the spontaneous emission, 
the density-density correlation of the atoms and the opti- 
cally induced dipole-dipole interaction of the atoms. For 
the Faraday interaction a detailed study of the reduction 



from three to one d imension is presented by S0rensen 
and S0rensen (2008). In essence this study confirms the 



treatment presented here provided that the gas is ideal 
and that the mo defunct ions u m (f_\_;z) are the solutions 
to the propagation equation including the index of refrac- 
tion. Little work has been done so far on the optically 
induced dipole interactions. 



Equations of motion The addition of the Hamiltonian for 
light Hl to the atomic part of the Hamiltonian H A and 
the interaction H[ nt discussed above allows to describe 
the propagation of light through the ensemble. As shown 
in Appendix |IX.C| the equation of motion is derived by 
introducing a rescaled time r = t — z/c and becomes a 
differential equation in space (z) instead of time t. 

For the beam splitter interaction, following this recipe, 
we find the equations of motion by calculating the com- 
mutator with Hbs 



d , . \g{z)\ 2 , x 9*(z)Q(t) , x 
— a L (z,t) = i^^a L (z,t) +i y y A V J a A (z,t), 



d_ 
dt 



: \m\ 2 



Ma*® 



a A (z, t) = z '"7a' a A(z, i) + i —^ a L (z, t). 



4A 



2A 



(35) 



The above equations can be solved analytically (Gor- 
shko v et al. 2QQ7c| Kozhekin et al.\ |2QQQ| [Ku priyanov 



et all 2005 



Mishina et al, 2007 Nunn et al. 2007). We 



will defer a discussion of the solutions till Sec. Ill.Cl where 
the spontaneous emission is included. 

Similarly we can find the equations of motion for the 
parametric gain: 



, x .g*(z)Q(z,t) t , x 
9 -a L {z,t)=i y l ^ U \(z,t\ 



d_ 

~dz 
d 



^ Mz,t)\ 2 , N 
dt a A {z,t) = -i — — — a A (z,t) 



(36) 



Again these equations have an analytical solution (Car- 
man et al., 1970 Raymer and Mostowski, 1981). 



For the Faraday interaction the equations of motion 
are much simpler when expressed in terms of x and p 
and read 

d , N g*(z)Q(t) , N 

— X L {z,t) = VA\Z,t\ 



d t N 

—p L (z,t) =0, 

d , N g*(z)tt(t) , N 

FtXA (z,t) = - g -^ PL (z^ 

PA\Z,t) =0. 



(37) 



dt 

Because the two momentum operators are conserved 
quantities these equations describe a QND interaction, 
where one can, e.g., make a measurement of the position 
operator xl after the interaction and thereby obtain a 



QND measurement of the atomic momentum operator 
p Al as will be further explained in Sec. |IV.A" 



The presence of the conserved quantities p A and pl 
makes it straightforward to solve the equations of motion. 
The only z dependence in the above expressions comes 
from the z dependence of the density n(z). We can then 
define the symmetric operators by 



X A 



Pa 



f dz ^n(z)x A (z) 
\Jjdzn(z) 

J dz y / n(z)p A (z) 
\Jjdzn(z) 



(38) 



If the transverse mode function includes all the atoms in 
the ensemble, these operators are equivalent to the sym- 
metric operators defined in Eq. (16). Similar integrated 
operators for the light field are 



X L 



Pl = 



f dt n(t)x L (t) 



'fdtn 2 (t) 
fd tn(t) PL (t) 
fdtfi 2 (t) 



(39) 



where we assume that Q is real. Expressed in terms of 
these variable the equations of motions have the simple 
solutions 



-* L,out — -* L,in5 
^A,out — X Ai [ n + K,Pl, in. 



(40) 



^A,out 



= p. 



A,in- 



The subscript "in" and "out" means input and output 
variables, e.g., X Ljin = X L (z = 0) and X L)OUt = X L (z = 
L) for light and X^in = X A (t = 0) and X A ^ out = X A (t = 
T) for atoms. The coupling constant k is given by 



dt 



2A 2 



dz \g(z)f 



ttujD 2 _ 



j it ^j izniz) . 



(41) 



In the last line we have inserted the expression for the 
coupling constant g(z) (31) and have denoted the dipole 
matrix element by D_ to indicate that it is the coupling 
constant for a— polarized light. 

As we shall see the coupling constant k (and the gener- 
alization of it) will play an important role for characteriz- 
ing the strength of the interaction regardless of the level 
scheme being used. Note, that k only depends on the 
total integrated density. This property can be shown to 
apply also to the A schemes by a s imple rescaling of the 
z coordinate (Gorshkov et al\ 2007c). We shall therefore 



for simplicity only consider a constant density below. 



10 



C. Theory Including Spontaneous emission 



In subsec. |II.B| the spontaneous emission was omitted. 
In this subsection we will give the theory including the 
spontaneous emission and discuss the solutions to the 
equations. 

Instead of repeating the calculations in subsec. 



II.B 



now with a non-zero decay, we note that if the sponta- 
neous emission from each atom is independent, it can be 
accounted for by making the substitution 



A m -► A m ± i 



In 



(42) 



where we have used the expression for the coupling con- 
stant g (31) and have introduced the spontaneous decay 



rate 70 = 4:UJ 3 \Dq\ 2 /3c 3 from the excited state |e) into 



the ground state |0) flMilonni and Eberly 1988) and the 
absorption cros s section for an atom ao = 3A 2 7o/27T7 
(Jackson, 1975). 



It is sufficient in our case to solve the operator equa- 
tions of motion (43) as "classical equations" with the op- 



erators replaced by complex functions (Gorshkov et al. 



2007b Raymer and Mostowski, 1981 ). The reason is that 



the equations are linear in the operators. The solutions 
to the operator equations will therefore be of the form 



in all the calculations performed so far. The choice of 
the sign is discussed in the Appendix |IX.D| and the re- 
sults are stated below. It is also important to note that 
the quantity (j x (f)) used to define the operator a^(r, t) 
in Eq. ( |2T| ), is not necessarily a constant. To ensure 
that the operator a^(r, t) has the right normalization we 
should always normalize by the time dependent expecta- 
tion value (j x (r,t)). The time derivative of (j x (r,t)} will 
introduce extra terms as discussed in Appendix |IX.D| 



Beam splitter interaction As discussed in Appendix |IX.D 
for the beam splitter interaction we should use the minus 
sign in the substitution (42), and the time derivative of 
{j x (r)) can be neglected because the strong classical field 
talks to an almost empty level. For a constant atomic 
density the equations of motion then become 



a L (z,t) 



m 



2 , N ig*tl(t) , N 
Tja L {z,t) + _ , a A (z,t), 



d_ 
dz 

d , , iMt)\ 2 , , igSl*(t) , . 
Wt a A {z,t) = -A±A^ aAM + ^a L (z,t), 

(43) 

where 7 is the total decay rate of the excited state and 
we h ave omitted the the noise operators given in Eq. 
(105). The admixed noise is vacuum in both equations 



of motion, and since the equations only couple annihila- 
tion operators to annihilation operators we can ignore the 
noise for the calculation of any normally ordered prod- 
ucts (see below). In the above equations the imaginary 
part of the first terms on the right hand side of each line 
refers to the phase shift caused by the index of refraction 
of the medium and the AC-Stark shift of the atoms, and 
the real part of these terms represent damping by spon- 
taneous emission. The last term on each line represents 
the coupling of light and atoms which is our main interest 
here. 

Solving these equation for a resonant field A = with 
no classical field Q = we find that the intensity is re- 
duced by a factor of exp(— <f), with the optical depth d 
given by 



d = L 



Mgf 



3A 2 7 Q 

2-7T7 



nL = naoL, 



(44) 



a>L,out(t) 



dz m(Q(t');t, L — z)a,A,m(z) + 



(45) 



where the operators are identified with hats for clarity, 
and the argument Q(t') indicates that the solution de- 
pends on the driving field at all times. The remaining 
terms denoted by dots in Eq. ( 45 ) are similar linear com- 



binations of the input light and noise operators. If we for 
instance solve the equations of motion with a complex 
function a,A,m(z) as the initial condition, we will due to 
the linearity of the equations get the same solution only 
without the hats 



&L,out(^) 



dz m(Cl(t),t, L — z)a,A,m(z). (46) 



We can thus obtain most of the solution (45) by simply 
inserting hats in the solution. Because we have ignored 
some input operators, however, the resulting operators 
will not necessarily have the right commutation relation. 
If there is no incident light all of these other modes will 
be in vacuum, and one can obtain the right commutation 
relation by adding a su itable amount of vacuum noise 
(Gorshkov et al. 2007b). Another way to see why the 
complex number equations are sufficient to obtain full 
information about the dynamics is to note that if all in- 
put modes are in classical coherent states we can take 
expectation values of the equations of motion (43) and 



obtain the same equations of motion for the mean val- 
ues. These classical equations of motion are thus iden 
tical to t he quantum equations of motion (Raymer and 
iMostowskil 1981). The equations of motion correspond 



to a general beam splitter relation, so that with coher- 
ent states as input states the output quantum states will 
also be a set of coherent states with amplitudes given by 
the mean values. Since any initial state can be expanded 
on the set of coherent states, the knowledge about the 
evolution of the coherent states obtained by solving the 
equations of motion for the mean values gives us com- 
plete information about the evolution for any quantum 
state. 

The equations of motion ca n in fact be solved analyt- 
ically (Gorshkov et a/., 2007c). If we consider the situa- 



tion where the only non vanishing initial value is a a we 
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find the solution for the output light in Eq. ( 46 ) with 



/^w /x x rid (—i)Q(z,t) J 2a+h h(o,t) djz 



L4(A-<J) 



x/ -ie l4> Jh(0,t) 



L 



(47) 



where /o(^) denotes the 0-th order Bessel function of the 
first kind. Here we have assumed that the coupling con- 
stant is real, which can always be done by absorbing any 
phase into the definition of f2, and we have replaced the 
coupling constant by the optical depth through the rela- 
tion (44). The function h{t,t') is defined by 



h(t,t') 



dt' 



dimnf 



4A 2 



r 



(48) 



and as we shall see below, this is also the function charac- 
terizing the strength of the Faraday interaction. Finally 
the phase <\> is defined by 



tan( 



_7_ 
2A* 



(49) 



Similarly we find for the output atomic variables 



«A,out(^) 



JO 



dtm(n*(T-i)]T-t,z)a Liin (t). (50) 



The fact that the kernel m(Q*(T — t)\T — t,z) is similar 
to the kernel in Eq. (|46|) i s a direct cons equence of time 
reversal symmetry ( Gorsh kov et al.\ |2QQ7c| ) . 

While the above solutions for the equations of motion 
are exact , they are sufficiently complicated and it is hard 
to gain any physical intuition from them. More insight 
can be gained by applying to the equations the Laplace 
transform in space using 



This equation now has a simple interpretation. Let us 
suppose that the incident light field is in vacuum so that 
we can ignore the last line in Eq. (52). The initial atomic 



state can then decay through two different mechanisms: 
either through spontaneous emission, or through coher- 
ent interaction with the forward light mode. Consider 
now the fraction in the first line of Eq. (52). The imagi- 



nary part of the first term, proportional to the detuning 
A will give rise to an unimportant phase, while the real 
part describes decay of the atomic excitation at an effec- 
tive rate (proportional to) 7(1 + d/2u). Here the decay 
rate 7 is due to spontaneous emission whereas the sec- 
ond term ^d/2u is due to the coherent interaction with 
the forward light mode. The optical depth d has been 
introduced in this expression through its relation with 
the coupling constant g in Eq. (44), and thus charac- 



terizes the strength of the coherent interaction with the 
forward light mode. For a sufficiently high optical depth 
d ^> 1 and sufficiently smooth atomic excitations u ~ 1 
the coherent interaction will dominate the spontaneous 
emission jd/u ^> 7 and we can obtain an efficient inter- 
face between atoms and light. 



Parametric gain interaction In case of the parametric gain 
interaction we need to use the plus sign in one of the sub- 
stitutions as well as include some terms associated with 
the time derivative of (j x (r,t)} as discussed in Appendix 
iDLDl We then find 

d , , .g*(z)Q(z,t) t 

a A\ z -> z )-> 



d -a L (z,t)=i 2A 



d_ 

m 



a A (z,t) 



\n(z,t)\ 2 A 



4A 2 



■r 



a A (z,t) 



(53) 



W n <*'*>«t( M ), 



2A 



■ ij 



a A (u,t) 



1 

7z 



where the z dependence of the Rabi-frequency is given 



/ dz e~ uz / L aA(z,t), (51) by the exp[z J dz' \g(z')\ 2 / '(A — ij /2)] dependence associ- 

JO Q-hfarl -unf h flno cY\n-ncrp> in flip T^rrm a ctq-H rirt r»f flip placcipQ 



where we have chosen the normalization so that u is a 
dimensionless number of order unity. We can then derive 
the equations of motion for the Laplace transformed vari- 
ables, which only couple operators with the same param- 
eter u. (Note, however that because the atomic operators 
only have support on z E [0, L] and not [0, oo[ different 
Laplace components are not orthogonal and care should 



ated with the change in the propagation of the classical 
field caused by the index of refraction and scattering of 
the medium. To arrive at these equations, we have as- 
sumed that the decay from the excited state goes into 
some states \a m ) different from |0) and |1), cf. Fig. [2jt>) . 
This model for the decay is not always best for real sys- 
tems used for parametric- typ e interaction, for example, 
in DLCZ-type applications (Duan et a/., 2001) a much 



be ta ken when applying these formulas ( [Gorshkov et al. 



2007c)). Since the Laplace transformed equation for the 
light field does not involve derivatives, we can then elim- 
inate the light field and obtain a single equation for the 
atomic operator 



d_ 
dt 



a A (u,t) 



W)l 2 



ig>/LSl*(t) 



a A {u,t) 



stronger decay to auxiliary states would decrease the ef- 
fective optical depth, and a large optical depth is crucial 
for the interface. The model used here is, however, the 
simplest possible model, and we therefore restrict ourself 
to this situation. The interpretation of the above equa- 
tion is similar to the interpretation of the equation of mo- 
tion for the beam splitter interaction (43), only in these 



(52) 



2ix(A 



2u)) 



0<L : m(t)- 



equations there is no coupling of the light to itself be- 
cause the light talks to an almost empty transition. The 
reason why there is no decoherence term in the atomic 
equation is that spontaneous emission drives atoms from 
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state |0) into some other state \a), where they are lost 
from the system. Collective states like the one in Eq. 
( p78J ) are, however, immune to removing one of the atoms 
in state |0) and this source of decoherence has therefore 
no effect within the approximations we are using here. 
See the work by Mewes and Fleischhauer ( 2QQ5[ ) for a 
more detailed discussion of the robustness of collective 
atomic states. 

These equations of motion can also be solved ana- 



lytica lly (Carman et al. 1970 Raymer and Mostowski 
1981). The solution is similar to the solution for the 



beam splitter interaction and we shall not go further 
into it here. Again it can be shown that it is possible 
to obtain a large coherent coupling with a large optical 
depth, by Laplace transforming the equations. By doing 
so one finds a strong gain term for the Laplace compo- 
nent with argument u, which is ~ d/u times the single 
atom scattering rate. For a large optical depth d we are 
thus dominated by the coherent interaction, regardless 
of the assumption about the final state after the decay, 
which we made above. Unlike the beamsplitter interac- 
tion, where one can work at any detuning, this strong 
coherent interaction only works in the far off-resonant 
regime, because the classical beam would be completely 
depleted if one works on resonance in a medium with 
large optical depth. 



Faraday interaction The Faraday interaction is a com- 
bination of the beam splitter interaction and the para- 
metric gain with equal weights, and the easiest way to 
obtain the equations of motion is therefore to just com- 
bine the results obtained in the previous two subsections. 
As discussed in Appendix |IX.D| the resulting equations 
of motion are then (assuming Q and g to be real) 



_9 

dz 



x Uz,t)= -^ 



2v^)Q«) wW _ ^SL,^ 



+ 7 a 



4A 2 + 7 2 ' 



Wz VhM = 4A2+V '- PA{Z) ~ <*■> ' -* PL{Z) 



4A 2 + 7 2J 

d XA(zt) - 2y/2Ag(z)Q(t) 7 » 2 () 

x A {z,t)- 4A2+7 2 PiW 2 (4A 2 +7 2 ) A( h 

7ft 2 



dt 



m PA{z ' t)= 4A 2 + 7 2 V ( ^ 



2(4A 2 +7 2 ) 



Pa(z), 



(54) 



where we have omitted the noise operators, which are 
given in Eq. (108). Again these equation are derived 



under the assumption that the decay goes to some aux- 
iliary state \a m ). For a treatment with decay back t o 
the interface levels, see for instance 



Duan et al. 



Hammerer (2006); Madsen and M0lmer (2004) 



(2000a); 



Let us now consider the solution of Eq. (1541) in the 
limit of a small damping. The Faraday interaction is 
only used with far off resonant light A ^> \fd^ since the 
classical field would be completely absorbed if we were 
working close to resonance. We will therefore ignore the 



first term in the evolution of the momentum operators 
Pl and pa, which is much smaller than the similar term 
in the evolution of xl and xa- Furthermore, a time de- 
pendent driving Cl(t) c an be accounted for by a simple 
rescaling (see Appendix |IX.E ), so we only consider a con- 
stant driving field Q. Because the quantities pl and pa 
appearing in the coupling to the operators xa and xl are 
conserved quantities apart from the small decays, the dy- 
namics effectively only involve the integrated operators 



X L = -±= f dt x L (t), 

vi Jo 



Pl= Vt'o dtpL{t) ' 

X A = —7= / dz xa(z), 
VL Jo 



(55) 



Pa 



-I 



dz pa(z), 



where the normalization is chosen as for single mode op- 
erators [X, P\ — i. In the limit of small scattering the 
resulting dynamics are then given by 



X L , out ^e-^/ 2 X L , in + ^P A , in , 

p ~ e -VL/2p . 

1 L,out ~ c J L,m? 

A A ,out « e-^/ 2 A A , in + kP lm , 

^A, out ~ e Pa, in- 



(56) 



Apart from the decay, this solution is the same as the 
results derived in Eq. fl40| with a minor modification of 
the coupling constant (|41[) which is now given by 



h(0,t) 



4A 2 



4A 2 



-T 



(57) 



where the function h(0,t) is the same as the function 
introduced for the beam splitter interaction in Eq. (48). 



In the limit of large detuning A the coupling constant 
here is the same as the one derived in Eq. (41). 



In the solutions above, the light and atomic operators 
are damped by factors of exp(— tjl/2) and exp(— 77^/2) 
respectively. The damping factor for the light is given by 



Vl 



i\g\ 2 L 



d 7 2 



i(a* + £) 



2 4A 2 + 7 2 " 



(58) 



Note that the optical depth used here is for linear polar- 
ization. The definition of d therefore differs by a factor of 
2 from (44) because we have taken g to be the coupling 
constant for circular polarized light and not for linear po- 
larization. From the expression above it is clear that we 
can ignore the decay of the light field if we use a suffi- 
ciently large detuning A ^> \fd^. A large detuning also 
reduces the coupling constant but this can be compen- 
sated by using stronger laser fields. 
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The atomic decoherence can be related to the coupling 
constant through 



drjA 



7dJdt\Q\ 2 

4A 2 + 7 2 



(59) 



This relation between <i, t)a-> and k is valid regardless of 
which assumption one makes about the final state after 
a decay. But different decay channels will have different 
effects in the equations of motion, and Eq. (54) is only 



directly applicable for the particular model considered 
here. Nevertheless, with a fixed interaction strength n 2 
we can always obtain negligible decoherence with a suf- 
ficiently large optical depth d. Conversely, for a given 
optical depth d there is always an optimal damping fac- 
tor t]l balancing the losses against the gains in some 
figure of merit (such as e.g. state transfer efficiency, 
light-matter entanglement etc.) whic h depends on the 



coupling strength ft, see for example (Hammerer et at. 
2004). 



Scaling with atom and photon number Further insight into 
the connection between the coupling constant, atomic de- 
coherence and dissipation of light can be obtained by ex- 
pressing these parameters through the number of atoms 
Na and photons Np of the classical field. Let us assume 
that the classical beam has a square profile with a cross 
section of area A. Using the definitions of the electric 



field (12) and Rabi frequency (31) we find 



/ 



dt\n\ 2 



3A 2 

''2ttA 



N P , 



(60) 



where j c (j q ) is the decay rate between the two states 
coupled by the classical (quantum) field with an emit- 
ted photon of the same polarization as the classical 
(quantum) field. The function /i(0, T) characterizing the 
strength of the interaction in the far off-resonant limit 
can be then expressed as 



«^ -(£)"«&** 



7 



A 2 4A 2 



T 



V 
N A N P , 



(61) 



where in the last line we have introduced the resonant 
cross sections for the scattering of the quantum and clas- 
sical fields (a m = 3A 2 7 m /27T7). The atomic and light 
decoherence can then be related by 



NaVa = N P T] Ll 



(62) 



which simply reflects the fact that the number of scat- 
tered photons is the same as the number of atoms which 
have scattered a photon. If the number of input clas- 
sical photons is much larger than the number of atoms 
Np ^> Na the decay of light is much smaller than the 
atomic decoherence and can be neglected. 



The function /i(0,£) = n 2 characterizes the solution 
of all three model systems considered here. Further in- 
sight into the reason for this can be gained by rewriting 
the equatio ns in t erms of rescaled dimensionless variables 
(Appendix |IX.E). 



D. Realistic Multilevel Atoms 

In the previous section we derived the main equations 
describing three types of light-matter interactions, the 
beam splitter, the parametric gain and the Faraday-QND 
type, Eqs. (30), (33) and (34) respectively, for simple 



few-level model atoms. In this section we will give ex- 
amples of how these interactions are commonly realized 
with real atoms. 



Faraday interaction Applying the Faraday interaction 
(34) to a full hyperfine level with many nearly degen- 



erate Zeeman states would at first sight seem to violate 
the simple two level approximation that we have used in 
the theoretical derivation. For alkali atoms, however, the 
full theory actually reduces to what we have derived in 
the previous sections if the detuning is much larger than 
the hyperfine splitting in the excited state. Consider the 
S\/2 ~^ -F3/2 transition in alkali atoms as indicated in 
Fig. [3J If atoms are optically pumped into one of their 
hyperfine ground state levels F, an off resonant probe 
will couple in general to three dipole allowed transitions 
F -> F' = F — 1, F, F + 1. Each of these transitions will 



contribute to the effective Hamiltonian in Eq. (25) de- 



scribing the interaction of a single atom with off-resonant 
light. It will be convenient to rewrite this Hamiltonian 

as 

H int = E(-\r)%E(+\r), 

where we introduced the atomic polarizability tensor op- 
erator (Deutsch and Jessen 1998| Happer, 1972 Happer 



and Mathur, 1967) 



EE E 



m' ,m" m" ,1 



\gm>)(gn 



and £^/ m // A D^J m denotes the dyadic vector product. 
The polarizability operator is a rank- 2 spherical tensor 

and can therefore be de- 



( Edmonds 1964 Zare 1988 



composed into irreducible tensor components, 



4ZJ 



^ (ao(A)T<°> + ai(A)T« + a 2 (A)T (2 ^ 



where the tensor operators T^ transform under rota- 
tions as a scalar, vector and matrix for k = 0,1,2 re- 
spectively. In this expression A = Apf = p+i is the laser 
detuning from the uppermost level and 2Dq = (J'\\D\\J) 



is the reduced dipole matrix element for the Si/ 2 —> Pi 



3/2 
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transition. It re lates to the spontaneous decay rate in- 
troduced in Sec. |lI.C| as 70 = 4uj 3 \D \ 2 /3c 3 . 

The real coefficients a/ c (A) follow from elementary cal 



culations ( Geremia et al. , 2QQ6| |Hammerer, 2006 Juls- 



|gaard , 2003 ; Kupriya nov et al] [2005 ) and are given in 



appendix |IX.F| The essential feature of these coefficients, 
which is proven in the appendix |IX.F[ is that for a laser 
detuning, which is large compared to the hyperflne split- 
ting of excited states, |A| ^> |A^ + i — A F /|, the rank-2 
tensor component vanishes, 02(A) — » 0. For the case of 
Cs the coefficients a&(A) are shown in Fig. pi 
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FIG. 3 Coefficients ao, ai, a2 for, respectively, scalar, vector 
and tensor polarizability versus (blue) detuning A of probe 
light driving the 6 Si/ 2 — >• 6P 3 / 2 transition in 133 Cs. The inset 
shows the relevant levels and energy scales. 

In the asymptotic limit the interaction Hamiltonian is 
thus given by 



D 2 - 



a 



g(-\r) • £ (+) (r) 



-CL\ 



£<"> 



if) • (. 



J x 



For an atomic ensemble, the first term will give rise to 
an index of refraction, while the second term accounts for 
the Faraday effect. For light propagating along z and the 
classical light polarized along x the Faraday Hamiltonian 
(34) can be easily derived. 

When the detuning is larger than the hyperflne split- 
ting, the interaction is thus essentially the same as for a 
spin 1/2 ground state, where CL2 vanishes exactly for all 
detunings. This can be understood by looking at how the 
energy levels of an alkali atom appear. The full Hamilto- 
nian can be written H = Hq + Hfs + Hhfs + ^int> 
where the four Hamiltonians represents the Coulomb, 
fine structure, hyperflne structure, and interaction with 
light. Normally one just considers the first three terms as 
an atomic Hamiltonian and does perturbation theory in 
the interaction Hamiltonian. When the detuning is larger 
than the hyperflne structure it is, however, more appro- 
priate to do the adiabatic elimination before treating the 
hyperflne interaction. Without the hyperflne interaction 



the optical fields only talk to the electron spin, where 
there cannot be a rank two tensor, since the spin cannot 
be changed by two, and a^ therefore only appears as a 
perturbation in Hhfs I 'A. 

The Faraday interaction of light with a true spin 1/2 
ground state atom, which therefore has only scalar and 
vector polarizability, can be achiev ed with the Ytterbium 
isotope 171 Yb, as was explored in (Takeuchi et al. 
20061). 



2007 



A-Systems Although for the beam splitter and the para- 
metric gain interaction, which require a A configuration, 
magnetic sublevels are sometimes used (iNovikova et all 



2007b), the most common implementation makes use of 
Sij2(F — I ± 1/2) hyperflne levels as ground states |0) 
and 11) and one of the P states as an excited state le), see 



(Chaneliere et al. 



Choi et al, 2008 



2005 2007 Chen et al. 



sukevich and Kuzmich 



Chou et al.\ |2004[ I2007J |2005| |Mat 



200812007b; 



2004). This approach gives ex- 



cellent results in zero magnetic field even if atoms are 
not optically pumped initially to one Zeeman substate . 
It has also been suggested by de Echaniz et al. (2008); 



Kupriyanov et al. ( 2005[ ); |Mishina et al. (2007) to make 
use of the second rank tensor polarizability to engineer 
effective A-schemes for the beam splitter and the para- 
metric gain interactions using two degenerate Zeeman 
states and polarized light. 



E. Ensemble in Magnetic Field 

So far, our theory for light matter interaction assumed 
degenerate ground state levels. For pure three level A 
schemes non-degenerate ground states do not make any 
difference, as the level splitting can be compensated for 
by choosing appropriate frequencies for light. In this sec- 
tion we will mainly deal with the Faraday interaction 
for atoms in an external magnetic field along the axis 
of atomic polarization (Fig. E|. The Zeeman splitting 
caused by this field can be advantageous in several re- 
spects. On the one hand, combined with the homodyne 
detection typically used in connection with the Faraday 
interaction, this results in low-noise AC signals, as will 
be detailed in sec. 



II.F 



On the other hand, it can also 
simplify and enhance protocols aiming for an e ffici ent 
creation of entanglement of two ensembles (Sec. IV) or 



between an ensemble and light (Sec. |VTJ ). 

The free Hamiltonian for an ensemble of atoms in a 
uniform magnetic field (Eq. ([24])) oriented along the x- 



Ho 



2 



dz (x 2 A (z) + p 2 A {z)) , 



(63) 



where ujl is the Larmor frequency. This Hamiltonian 
generates Larmor precession of the transverse spin den- 
sity components, xa(z) andp^(^), about the x-axis. The 
full Hamiltonian describing Faraday interaction in a mag- 
netic field is H = Hq + Hp, where Hp is given in Eq. 
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FIG. 4 A magnetic field along the axis of polarization causes 
Zeeman splitting ujl of ground state levels. Photons will be 
scattered from the classical light, driving the tv transitions, 
to the quantum field, coupling to the cross a transitions, at 
sideband frequencies cjq ± ojl from the carrier frequency cjq • 



(34). In an interaction picture with respect to Hq this 



Hamiltonian is, 



H-,-j 



g*(z)Q(z,t) 

dZ v^A PL{Z) 



x (cos(oj L t)p A (z) +sin(uj L t)x A (z)). (64) 



Operators xa and pa refer now to spin components in 
a frame rotating at ujl about the x axis. For the sake 
of simplicity we use in the following the same symbols 
for canonical operators in both frames, as it will be clear 
from the context, to which one we are referring. In the 
rotating frame the canonical operators for transverse spin 
components are related to the spin components in the lab 
frame via 



Jy( z )/V n ( z ) = cos(uj L t)x A - sm(u L t)pA, 
Jz(z)/V n ( z ) = cos(u L t)p A + s'm(u L t)x A . 

with the number density of atoms n(z). 

The Maxwell- B loch equations in the rotating frame are 
accordingly, 



x L (z,i) = - 

p L (z,t) = 0, 



V2A 



d_ 
dz 
_9 
dz 
d , N g*(z)to(t) , x , x 

—x A {z,t) = cos(u L t)p L {z,t), 

d , N g*(z)Q(i) , N , N 
—p A (z,t) = — sm{u L t)p L {z,t). 



(cos(uj L t)p A (z) +sm(oj L t)x A (z)), 



F. Quantum measurement and feedback 

In this section we will deal with measurements which 
can be done on light and then fed back onto atoms. 
We will focus on homodyne detection of light, which is 
of importance for experiments using Faraday interaction 
and EIT but also briefly describe photon counting, which 
is used in combination with parametric gain and beam 
splitter interactions. 



Homodyne detection of light The discussion of the quan- 
tum interface in the language of canonical variables for 
light is most fruitful because these variables can be mea- 
sured with almost perfect efficiency by the balanced ho- 
modyne technique. We will concentrate here on the po- 
larization homodyne version which is relevant for several 
protocols described in the article. In particular in the 
context of the Faraday interaction the polarimetric mea- 
surement of light is an important tool. Balanced homo- 
dyning employs overlapping the quantum field of interest 
with a strong coherent field, a local oscillator, on a 50/50 
beamsplitter and measurement of the difference of the 
power in the two outputs. In its polarization version as 
shown in Fig. |5j the local oscillator field and the quan- 
tum field are overlapped on a polarizing cube so that they 
have orthogonal polarizations and the role of the beam- 
splitter is played by a polarizing beamsplitter which splits 
the light into 45° and —45° modes. The measurement of 
the differential power of these two modes corresponds 
to the measurement of the S y Stokes operator, whereas 
with an extra A/4 plate in front of the beamslitter the 
S z Stokes operator is measured: 

Sx(t) = (a\ x a Lx - a\ iV a Liy )/2, 

by\t) ~ ( a L,+45° a L,+45° — a L,-45° a L,-45° )/* 

= ( a L, X a L,y+ a L,y a L, X )/ 2 i 
S z(t) = ( aJ L,a+ a L J a+- a L J a- a L,a-)/ 2 

= -K a L,x a L,y ~ a L,y a L,x)/ 2 - 

The third Stokes operator S x is equal to the total pho- 
ton number of the strong field for the case of the strong 
coherent field \a) in linear x polarization, such that 



(Sx(t)) 



(65) 



Now the atomic momentum operator, i.e. the spin pro- 
jection along the axis of light propagation, is not con- 
served anymore and the overall interaction is not of QND 
character. Integration of these equations becomes some- 
what more involved than before. We will resort to this 



2 /2. In this case the measurement of the 



problem in sections IV and VI 



two other Stokes operators amounts to a homodyne de- 
tection of y polarized light with the coherent field in x 
serving as the local oscillator, 



The homodyne detection can also provide an excellent 
suppression of the technical (classical) noise if the fre- 
quency of the local oscillator and the quantum field differ 
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FIG. 5 Polarimetric measurement of light 

by ool lying in the radio frequency domain, see Fig. [4] 
In this case the relevant canonical variables are encoded 
in sideband modulation modes of ^/-polarized light which 
are read in the cos(cjl£) and sin(o;x,t) components of the 
photodetector output: 



x Ls = 

Pl = 



I dtcos(u)Lt)xL(t), 
/ dt cos(lj L t)p L (t), 
/ dtsin(ujLt)xL(t), 
/ dt sm(u> L t)p L (t). 



(66) 



These components of the photocurrent can be measured 
by lock-in amplifiers. The bandwidth of this measure- 
ment can be adjusted to BW « r _1 where r is the optical 
pulse duration. In this way fluctuations at all frequencies 
outside this bandwidth are effectively irrelevant. In the 
case where atomic ground state levels are non-degenerate, 
e.g., are split by the Larmor frequency ujl as discussed 
in Sec. |II.E| and shown in Fig. [4| the atoms couple to 
the sidebands of light and the entire measurement and 
interaction can be encoded at sideband frequencies ±co>l, 
as in several experiments described later in the article. 



Feedback Another important tool in many quantum in- 
formation protocols is the feedback of results of measure- 
ment of light onto atoms. The theory of quantum feed- 
back is a wide field on its own, especially in the case of 



continuous measurement and feedback (Thomsen et al. 



2002a). Here we deal with a relatively simple measure- 



ment and feedback scheme, where light observables of 
the type (66) are measured by integrating a photocur- 



rent over the whole pulse duration and the measurement 
result, a single number, is fed onto the atoms. The oper- 
ations which need to be done on the collective spins are 
small rotations about the y or z axis, i.e. small tilts of 



the collective spin. In the language of canonical opera- 
tors Xa, Pa this amounts to displacements in the phase 



plane (Arecchi et al. 1972). 



In this case - feedback of integrated measurement re- 
sults via displacement operations - a simple rule can be 
applied for describing the overall effect on a state of the 
atoms. Assume that the state of the system is described 
by certain input-output relations of the type (40). If a 



quadrature of light, say Xl, is measured and the corre- 
sponding measurement result £ is used to displace the 
atomic state, i.e. to tilt the collective atomic spin, in 
such a way that the mean of, say Pa, is transformed as 
(Pa) — >• (Pa) + g£, then the statistics of Pa after the 
feedback operation can be calculated from 



^U,finai = Pa + gX L , 



(67) 



that is, one simply needs to add the measured observable 
multiplied by the gain to the operator which is subject to 
the feedback. This rule for describing the feedback holds 
strictly as an operator identity, irrespectively of the state 
of the system being Gaussian or Non-Gaussian. 

The proof is most easy in the Schrodinger picture. As- 
sume the state of some bipartite system is Pal, where 
the indices refer, e.g., to atoms and light respectively. 
A measurement of Xl gives a result £ with probability 
Pt = (^a{pal}\0^ where X L \£) = f |f), and the state 
of the system A collapses to 

The feedback affecting the desired displacement is de- 
scribed by a unitary transformation of the state of the 
system A, 

which gives in the ensemble average over all possible mea- 
surement results the final state 



pf = J d£ Pi pT = J diie 



ig£,x A 



(ZlPALlOe 



-ig£X A 



I 



dm 



igX L x, 



PALe 



-igX L X A 
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= tr L { 



igX L XA 



PALe 



-igX L XA 



}• 



From this equation all the moments of Pa can be calcu- 
lated as 

{P!) = ^a{p1p { 1 ) } 

= tY A L{(PA+gX L ) n pAL}, 

where the cyclic property of the trace was used in the 
second equality. This justifies the rule given above. 

Note that for Gaussian states, for which it is enough 
to keep track only of the first and the second moments 
in order to have the full knowledge of the state, the sim- 
ple linear transformation of operators as above is exactly 
equivalent to a full description in the Schrodinger picture 
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e.g. on the basis of the Wigner function. The measure- 
ment of and the feedback on more than one mode can 



be described by an immediate generalization of (67) as 
shown in (Hammerer et a/.l|2005a| Sherson et al. 2006c), 



provided the measurement involves commuting observ- 
ables only. 



Photon Counting In some cases measurements in the 
Fock state basis, i.e., photon counting are convenient for 
characterization of the interface performance. As dis- 
cussed in Sec. [Vj storage and retrieval of some nonclas- 
sical states can be characterized by the measurement of 
the second order correlation function g 2 (l,2) which is a 
normalized pr obability of ph oton counts at two points 
in space-time (Loudon, 2004). Atoms typically used for 



the interface are Rubidium and Caesium and the corre- 
sponding spectral lines are around 780nm and 850nm re- 
spectively. In this spectral domain commercial avalanche 
photodiodes (APDs) typically have quantum efficiency 
around 40 — 50 % and the dark count rate of a few hun- 
dred per second. Such parameters are sufficient to de- 
termine non-classical correlations via g 2 (l,2) which are 
insensitive to losses if dark counts are neglected. 



G. Other Strategies 

In the discussion so far we have focussed on what we 
consider to be the main protocols in this field. There are, 
however, numerous variations of all of these protocols, 
some of which we will briefly discuss here. 



Non-copropagating beams. In the derivation above we 
have only considered the situation, where the quantum 
and classical light fields are co-propagating. For many 
applications of the beam splitter and parametric gain 
interaction t his assumption is, however, not necessary 
(|Balic et a/.[ |2005| [Braje et al] |2004| |Chaneliere et~aT > 
2005). For instance, an incoming photon which is ab- 



sorbed in an atomic ensemble by using a co-propagating 
classical field generates an excitation of the form of Eq. 
(18). The reason why one can later retrieve this quan- 



tum state is constructive interference. During readout all 
atoms will radiate in phase in the direction of the classi- 
cal laser and this is the effect, which allows for efficient 
interfaces between atoms and light in the limit of a large 
number of atoms (high d). If the photon which is ab- 
sorbed has a different direction than the classical drive 
field, the generated atomic excitation will still have the 
form of Eq. (18). The only difference is that the state 



will have a phase factor exp(zA/c • fi) on the component, 
where the zth atom is in state |1), with Afc being the dif- 
ference between the k- vectors for the two fields. In order 
to have constructive interference in the readout process 
the difference in the k- vectors of the outgoing photon and 
the classical field in the read out process should cancel 



the phase factor imprinted on the atoms in the first step 
of the protocol (see also Sec. |V.D for a discussion of the 
effect of difference in the k- vectors). 

Expressed in different terms, the initial atomic state 
|000 — 0) has a homogeneous phase corresponding to a 
zero momentum state. The initial process involving the 
absorption of a photon from one beam and the emission 
of a photon into a different beam imprints the difference 
momentum Afc onto the atomic spin wave. Construc- 
tive interference is achieved for processes returning the 
atoms to the initial zero-momentum spin wave |000....0) 
and the differential momentum in the read out process 
must therefore carry away the momentum in the spin 
wave. To achieve momentum conservation (or equiva- 
lently phase matching) the total momentum of all the 
absorbed photons should thus match the total momen- 
tum of all the emitted photons. 

A disadvantage of using non co-propagating beams is 
that the storage into non-symmetric modes limits the 
storage to the time it takes atoms to move a distance 
~ l/ |Afc|. In particular for room temperature gasses (see 
Sec. III. A) this, as well as the differential Doppler-shift, 



makes it undesirable to use geometries, where the beams 
are not nearly co-propagating. Even for cold atomic en- 
sembles this effect limits the stor age time, as was ob- 
served experimentally by Zhao et al. (2009a). On the 



other hand it is often a major experimental advantage not 
to have the classical and quantum beams co-propagating 
since this makes it much easier to count photons in the 
quantum beam. 



Interaction based on phase shift It is instructive to have 
a look at the Faraday interaction from another perspec- 
tive. Consider first Fig. J6^a) where, as before, the quan- 
tization axis is taken along x, the direction of polariza- 
tion of both atoms and light. Selection rules for dipole 
transitions dictate that strong classical, x-polarized light 
with amplitude a x ~ a drives the \±) x — » \±) x transi- 
tions, while the \±) x — >> |=f) x cross transitions are cou- 
pled to the weak quantum field a y in ^/-polarization. In 
this picture it is evident that the Faraday interaction is 
the sum of the beam splitter and the parametric gain 
interaction, cf. Fig. [2] The same level configuration 
can also be looked at by taking the axis of quantiza- 
tion along the z direction, as shown in Fig. Rub), where 
\±)z = (\+)x ± \-)x)/V%- Light propagating along this 
direction naturally couples with its circular components 
a a ± = (a x ± ia y )/V2 ~ (a ± ia y )/y/2 to the cross tran- 
sitions \±) z — >> |=f)s only. The off-resonant coupling will 
thus give rise to AC Stark shifts on atomic levels \±} z 
depending on light intensities of a 1 * 1 -polarization com- 
ponents. Vice versa, light polarization is thus rotated 
according to the difference in level populations |=b) 2 . 

From this observation it is clear, that all that is re- 
quired for a Faraday interaction is a mechanism of non- 
destructive measurement of level populations via phase 
shifts of light. Making use of the vector polarizability for 
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(b) 



X 



potential application for four wave mixing. For a review 
of this see Fleischhauer et al. (2005). 



|0> = |+)x h)x=|l) \-) z \+) z 

FIG. 6 Level scheme for Faraday interaction: (a) For the 
quantization axis along x: Atoms are polarized to \+) x , laser 
light is linear polarized along x and drives the up-transit ions, 
while the quantum field in y polarization couples to the cross- 
transitions, (b) The same interaction with the axis of quan- 
tization taken along z: Circular light components of equal 
intensity cause AC Stark shifts of equally populated states 



probing the collective hyperfine spin-angular momentum 
as described above, is therefore just one way to achieve a 
Faraday or QND interaction. In addition, several propos- 
als and experiments pursue the idea to probe coherences 
of the pseudo-spin c onsisting of the two Si/2JF = Jzbl/2) 
hyperfine levels, see (Chaudhury et al] 2006 Oblak et al] 



2005||Petrov et al] [2007} | Windpassinger et a/.||2008| ). As 
shown in Fig. [7| at a certain detuning the phase shift of 
the probe light due to F — I ± 1/2 —> F r transitions ex- 
actly cancels for equal populations of F = / ± 1/2 levels. 
Any imbalance of populations will yield an interferomet- 
rically detectable phase shift. 




FIG. 7 QND interaction for measurement of pseudospin com- 
posed of hyperfine ground states of Cs: (a) Level scheme 
for 133 Cs with probe light tuned in between hyperfine lev- 
els F — 3,4 of the 6 2 Si/2 ground state, (b) Differential 
phase shift $ due to the F = 3 -► F' = 2,3,4 and the 
F = 4— ^F' = 3,4,5. At magic frequencies the phase shift 
vanishes. 



Other Hamiltonians and level structures. There are sev- 
eral possibilities involving more complicated atomic level 
structures than the ones shown in Fig. |2J A particular 
example is the so called double A-systems, with two ex- 
cited states. An interesting feature of this system is its 



In the derivation of the theory we adiabatically elimi- 
nated the excited state to arrive at the effective ground 
state Hamiltonian. As discussed in Subsec. III.DI this adi- 
abatic elimination in general leads to a Hamiltonian in- 
volving spherical tensors of rank zero, one, and two with 
strength characterized by the three coefficients ao, cl\ and 
G&2- The three protocols that we have mainly considered 
thus correspond to suitable initial states and particular 
combinations of these spherical tensors. By adjusting the 
detuning as well as laser polarizations and atomic initial 
state there is, however, a lot of freedom in varying the rel- 
ative strength and effect of the different tensors, which 
allow for a richer dynamics. Kupriyanov et al. (2005) 
and Mishina et al.\ (|2007|) considered how the higher or- 



der tensor operators modify the equations of motion and 
in particular how the Faraday inter action is influenced by 
the rank two tensor. For instance, Mishina et al. (2007) 



it was shown that a particular choice of detuning removes 
the AC- Stark for a detuned beam splitter interaction and 
thus removes the need to adjust the frequency of the clas- 
sical driving field in order to keep the field in two photon 
resonance with the AC-Stark shifted transition. 

To arrive at the Faraday interaction we just combined 
the beam splitter and parametric gain interaction with 
the same strength, but in Mishina et al. (2006) it was 
shown that an arbitrary combination of the paramet- 
ric gain and beamsplitter interaction can be obtained by 
choosing suitable initial conditions, detunings and com- 
binations of the elements of the spherical tensor. 

An example of a protocol where the light-atom inter- 
face involves excited states without adiabatic elimination 
is shown in Fig. [IJa) flHald et al] |1999 



Kuzmich et al. 



1997). Another protocol of this type is considered in Sec. 



V.E[ where we discuss spin echo techniques. A disadvan- 
tage of such protocols is, however, that the storage time 
is limited by the coherence time of the optically excited 
state, which is often shorter than the coherence time of 
ground states. 



Optical cavities. The key parameter in characterizing the 
applicability of an atomic ensemble for a light matter 
quantum interface is its optical depth which for a free 
space ensemble is limited by the size and atomic interac- 
tion. An alternative strategy is to use multiple passes of 
the light through the atomic ensemble by enclosing the 



ensemble in an optical cavity (|Black et al. 2005 Dan 



tan et aJ.||2QQ5||Josse e£ a/.||2004||Simon et al. ||2007a|b| 



Thompson et al, 2006). In this case the parameter char- 



acterizing the usefulness of the system is the cooperativ- 
ity parameter C = iV a9c/ ' n cli where g c is the coupling 
constant for a single atom to the cavity mode, and n c 
is the cavity decay rate. The cooper ativity can also be 
expressed as C ~ Td, where T is the finesse of the cavity 
which roughly equals the number of passes that the pho- 
ton makes through the cavity (|Gorshkov et al] |2007b[) . 
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The gain achieved by using a cavity thus equals the num- 
ber of round trips. 



Non Gaussian operations. In this review we consider the 
quantized light fields which are much weaker than clas- 
sical control and driving light, and the quantum fluctua- 
tions of the atomic ensemble which are much smaller than 
the mean spin. In this limit we only include the lowest 
order terms in the atomic and light field operators a a and 
cll- Because there are no first order terms the effective 
Hamiltonian will be quadratic in the harmonic oscillator 
operators. The operations which may be performed thus 
fall into the class of Gaussian operations and the solution 
of the equations will in general be a Bogoliubov transfor- 



interactions of atoms in optical lattices (Muschik et al. 



2008) or the so called Rydberg blockade, where the ex- 
citation of a single atom to a Rydberg level blocks the 
excitation of other atoms, and therefore creates a uni- 



form long distance interaction (Lukin et al. 2001). Al 



ternatively one can exploit the fact that atomic ensembles 
are particularly well suited for "catching" traveling pho- 
tons and then afterwards transfer the excitation to some 
other system for processing the information. A proposal 
along these lines is presented by Rabl et al] ( |2006 ) based 
on a transfer of excitations from an ensemble of dipolar 
molecules to a solid state system. A review of techniques 
for achieving other types of operations , e.g., Kerr inter- 
actions, is given by Fleischhauer et al. (2005). 



mation of t he incident mode operators (Braunstein and 

van Loock| |2005[). For any input state with a Gaussian H. Summary of the theory 



Wigner function the output Wigner function will also be 
Gaussian. The main advantage of using atomic ensem- 
bles is that the dynamics resulting from these Gaussian 
operation are collectively enhanced so that a perfect op- 
eration is achieved in the limit of large optical depth. 
While the resulting dynamics allow for a variety of quan- 
tum information protocols to be performed, such as quan- 
tum teleportation and quantum memory (sec [V] and VI), 
the fact that higher order terms are not collectively en- 
hanced limits the applications for quantum information 
processing. In particular it is known that Gaussian op- 
erations alone do not allow for distillation of entangle- 



ment from Gaussian states (Eisert et al, 2002 Fiurasek 



2002 Giedke and Cirac, 2002) and that algorithms for 



efficient classical simulation of any evolution involving 
onl y Gaussian op eration and Gaussian initial states ex - 
ist flBartlett et~al\ |2002| |Lloyd and Braunstein) |1999). 



These limitations may, however, be avoided by combin- 
ing the Gaussian operation with photon counting ( Genes 
and Berman[ |2006| |Neergaard-Nielsen et al] |2006| |Our 



joumtsev et al. 2006|). In the pioneering paper by |Duan 



et al] ( |2001| ) such photon counting techniques were pro- 
posed as a means for quantum communication over long 
distances using the probabilistic entanglement protocols 
discussed in Sec |IV.C[ Such techniques could in principle 
also allow for even more advanced quantum information 
protocols to be implemented. 

The fundamental obstacle for directly achieving non 
Gaussian operations is that they rely on an interaction 
between individual excitations. Such interactions be- 
tween the excitations can for instance be achieved if two 
photons interact with the same atom. But since this is 
essentially a single atom effect it is not enhanced by a 
large optical depth. An approach which allows for an en- 
hancement of this nonlinear effect by optically imprint- 
ing a Bragg mirror which localizes excitation similar to 
an optical cavity is explored in Bajcsy et al. (2003) and 



Andre and Lukin ( 2002 ) . An alternative approach to non 



Gaussian operations is to engineer a strong interaction 
between excitations stored in different atoms. Interest- 
ing proposals in this direction are to use the collisional 



We have presented a detailed theory for the quantum 
interfaces between light and atomic ensembles. In par- 
ticular we have presented a unified theory for the three 
model systems presented in Fig. [2ja)-(c), the beam split- 
ter, the parametric gain, and the QND (Faraday) inter- 
action. The three systems have distinct features, but 
are also interconnected, e.g., the Faraday interaction is 
just a combination of the beam splitter and parametric 
gain interaction with equal weights. Most importantly 
all three systems achieve ideal operation in the limit of 
high optical depth d. This feature can be understood 
as constructive interference or collective enhancement of 
the coupling: the coupling between the state where all 
atoms are in the ground state and the collective state 
(18) scales as \/Na, whereas spontaneous emission is a 



single atom effect, which is independent of the atom num- 
ber. We emphasize, however, that the theory is derived 
without the optical broadening present in many experi- 
mental realization. One therefore cannot directly replace 
the optical depth d appearing in the equations of this 
section with the actual measured optical depth in the 
presence of inhomogeneous broadening. Nevertheless the 
optical depth in one version or the other remains the key 
parameter for characterizing the usefulness of an atomic 
ensemble (see Sec. VII for a discussion of inhomogeneous 
broadening). 

For all three types of interaction, in the far off- 
resonant limit A ^> c?7, the strength of the coupling is 
parametr ized by e xactly the sa me function h(0,T) « ft 2 , 
cf. Eqs. ( 46|47|48 ) and ( 56|57 ), which is most easily seen 
by rewriting the equations in the dimensionless form as in 
Appendix |IX.E In this limit the decay and phase shift 
of the light can be ignored and the coupling constant 
can be related to tja, the spontaneous emission proba- 
bility per atom, through ft 2 = t]a^^ which we explicitly 
derived for the Faraday interaction, but which also ap- 
plies to the other systems. From this expression one thus 
directly sees that the spontaneous emission can be elimi- 
nated for high optical depth. In a special case of the beam 
splitter interaction, the resonant EIT setting, these argu- 
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ments are not directly applicable, but the ratio between 
the desired evolution and spontaneous emission, i.e. the 
constructive interference discussed above, is completely 
independent of detuning ( [Gorshkov et al. 2QQ7a|c| ) . 

It is instructive to discuss the bandwidth of the quan- 
tum interface, i.e., how fast an operation such as a stor- 
age or a read process can be performed. In the theory 
we have adiabatically eliminated the excited state, which 
means that the shortest time r on which an operation can 
be performed is limited by the low saturation parameter 
condition s <C 1 and the condition on the value of the 
coupling constant n necessary for a particular process. 
Bearing in mind the relations tja = T^s, r = (£?W) _1 , 
we can draw some general conclusions on the Fourier lim- 
ited bandwidth BW of the process. For protocols where 
k 2 ~ 1 the low saturation condition yields the limitation 
on the bandwidth of the light pulse, BW <C jd (again 
the arguments given here only directly apply in the far off 
resonant limit, but the conclusion is also valid for EIT). 
Going beyond the low saturation regime s ~ \ allows 



to increase the bandwidth somewhat (Gorshkov et al. 



2008), but it remains limited by BW < jd. For a typical 
alkali atomic ensemble the bandwidth of the order of 10 
MHz can be achieved. This scaling of the BW provides 
an upper limit, B W <C 7V^, e.g., for an atomic EPR en- 
tanglement protocol discussed later in the article where 
the conditions n 2 ^> 1,7? « 1/y/d has to be met. Note, 
however, that the term bandwidth is also often used in 
connection with the number of different modes which can, 
e.g., be stored in an atomic ensemble. This is a different 



question, which we will return to in Sec. V.E 



III. ATOMIC MEDIA FOR QUANTUM INTERFACE 

A few common requirements can be formulated for all 
ensemble based interfaces described in the previous sec- 
tions. A long lived (ground) state of atoms is commonly 
used. This could, e.g., be Zeeman levels or hyperfine 
levels. The ensemble should be initialized to a polarized 
state (coherent spin state), that is, one of the ground sub- 
states should be populated by optical pumping or other 
means. Most importantly, the sample should have a large 
resonant optical depth d. Up to date, experimental re- 
alizations of the ensemble based interfaces utilize alkali 
atom gases at room temperature, alkali atoms cooled and 
trapped at temperatures of a few tens or hundreds of mi- 
crokelvin, or impurity centers in solid state. Below we 
describe these and other media used for quantum inter- 
faces. 



A. Room temperature gases 

A gas sample of alkali atoms is one of the simplest 
atomic ensembles to have in the lab. Surprisingly enough 
such an object can also work very well as a quantum 
memory, if proper care of decoherence is taken. The ther- 



mal motion and associated with it Doppler broadening 
are not necessarily a problem. For the Faraday interac- 
tion, the Doppler broadening plays little role if the de- 
tuning is much greater than the Doppler width (200-300 
MHz for Cesium or Rubidium). For other protocols such 
as the beam splitter interaction the Doppler broadening 
has a detrimental but still tolerable effect as discussed in 
Sec. Ivm 

In addition to the Doppler broadening the atomic mo- 
tion also leads to changes in the atomic positions. A 
quick glance at the solution to the beam splitter interac- 
tion (47) and (50) reveals that the atomic operators with 



different longitudinal coordinates experience different dy- 
namics. The atomic motion in the process of interaction 
leads to washing out of these spatial modes which is a 
much more pronounced problem for the beam splitter in- 
teraction, as compared to the Faraday interaction. 

In order to reduce the deleterious effect of atomic 
motion a buffer gas is usually used in the experiments 
which utilize the beam splitter interaction in gas cells at 



room temperature (Eisaman et al. 2005 Novikova et al. 



2007a). Adding a few torr of a noble gas allows to suffi- 



ciently localize the diffusive motion of alkali atoms. An 
extra benefit of this approach is that it prevents alkali 
atoms from depolarizing collisions with the walls of the 
cell which otherwise could lead to a rapid decoherence. 
The lifetime of the atomic memory in cells with a buffer 



gas can reach milliseconds (Novikova et al.\ 2007a). Note 



that for some protocols, collisions of atoms in the excited 
state with the buffer gas lead to an energy redistribution 
of scattered photons which may lead to larg e errors in th e 
absence of careful spectral filtering ( |Manz et al., 2007). 
Despite the difficulties, alkali atom cells with a buffer gas 
has been successfull y used for experimen ts on quantum 
memory using EIT ( Eisaman et al. , 2005[ ) . 

The effect of atomic motion for the Faraday interac- 
tion can be almost completely eliminated. As follows 
from the propagations equations (40) the Faraday inter- 



action couples light to a symmetric atomic mode defined 
in Eq. ( [16] ) . In this case atoms with different coordi- 
nates z along the direction of light propagation couple to 
light in the same way. Hence the atomic motion along 
z does not affect the interaction. The transverse motion 
of atoms along x and y axes will affect the performance 
if the spatial profile of the light beam is inhomogeneous 
which is almost always the case. This effect can be re- 
duced in two extreme cases, either for times short com- 
pared to the motion time, or in case when the duration 
of the light pulse is so long that atoms have a chance to 
cross the beam many times during the interaction and 
the effect of motion averages out. The latter was the 
case of the experiments by Julsgaard et al. (2001, 2004a) 



and Sherson et al. (2006a), where pulses of about 1 msec 
duration have been used. 

The possibility to eliminate the effect of atomic motion 
on the efficiency of the interface based on the Faraday 
interaction has allowed to conduct high-fidelity experi- 
ments with room temperature Cesium atoms (Julsgaard| 



21 




(2001 2004a); Sherson et al. (2006a) the bandwidth of 



FIG. 8 Paraffin coated Caesium cell. 



et al] |2001[ |2004a| [Sherson et al] |2006a| ). Atoms were 
contained in cells with a paraffin coating of the internal 
walls (Fig. [8]). Such coating has been used in precision 



magnetometers for the past three decades (Alexandrov 



2003| Groeger et al, 2006), and ground state coherence 
times of up to a second have been demonstrated. In 
paraffin coated cells atoms can withstand tens of thou- 
sands collisions with the cell walls before significant spin 
depolarization occurs. Since it is the number of collisions 
with walls that matters, the larger is the cell, the longer 
is the quantum memory lifetime. As discussed in details 
in Sections IV, V, and VI, quantum memory time of the 
order of several milliseconds has been achieved in cells 
with dimensions 25 x 25 x 25 mm 3 . 

Room temperature ensembles of Cesium atoms of a few 
cubic cm size contain about 10 11 — 10 12 atoms. For off- 
resonant Faraday-type interfaces the Doppler broaden- 
ing does not af fect t he ensemble effective resonant optical 
depth (see sec. VII) which is the same as for atoms at rest 
reaching the values of the order 50 or even higher. The 
experimental challenge lies with the fact that the quan- 
tum spin noise of such an ensemble is \/Na ~ 10 -6 . 
In order to reach the level of the spin quantum noise, 
all types of technical spin fluctuations, such as driven by 
stray magnetic fields or fluctuations of the lasers used for 
optical pumping, have to be reduced below this level. 



The solution to this problem used in Julsgaard et al] 
(|2001[ |2004a[); |Sherson et al] (|2006aD has been to apply 



a bias magnetic field along the direction of the collective 
atomic spin. A field of the order of one Gauss provides 
the Zeeman splitting of the ground state uol of a few hun- 
dred kHz. As discussed in Sec. III. El this means that the 
collective transverse components of the spin which cor- 
respond to the atomic canonical variables rotate at the 
Zeeman frequency. At the same time, as described in 



Sec. II. F[ canonical variables of light which couple to the 
rotating atomic spin can be measured via homodyne de- 
tection also at the Zeeman frequency. Thus all relevant 
variables for light and atoms are now encoded at a fre- 
quency of a few hundred kHz. At these frequencies tech- 
nical noise can be reduced below the 10 -6 level, so that 
both spin and light fluctuations are dominated by quan- 
tum noise. In practice the photocurrent detected in the 
homodyning process is measured with lock-in amplifiers 
which allow access to the light variables ( 66 ) encoded at 



the memory has been reduced to around 1 kHz for the 
reasons discussed in Sec. III. HI 



B. Cold and trapped atoms 

Cold and trapped ensembles of alkali atoms have been 
among the first atomic objects to be used for quantum in- 
terfaces with light. The first experiment mapping quan- 
tum properties of light onto atoms was performed with 
Caesium ato ms in a MOT, a magneto-optical trap (Hald 
et 



1999). A MOT provides a relatively simple way 



to achieve a cold atom sample suitable for the quantum 
interface, however it also has its limitations. A typical 
resonant optical depth in a MOT lies in the range be- 
tween 2 and 10 which is not very high. Another consid- 
eration concerns the transverse crossection of the light 
beam which couples to the atoms. In most of the exper- 
iments on interfaces which use a MOT, light is focused 
down to a few tens of microns which is much less than 



a typical MOT crossection of a millimeter (Chaneliere 



~eTal\ [2005] jChen et a/.||2008l|Chou et a 1. 1 [2004 1| Pant an 



et al] 2005 Simon et al. 2007b). Such geometry lim- 
its the atomic memory lifetime to the transient time it 
takes the atoms to leave the probe volume. For a typ- 
ical MOT temperature of 100 /iK this time is around 
hundreds of microseconds, as was demonstrated by Zhao 



et . 



(2009a). If, on the other hand, the beam crossection 



is such that the light couples to the entire MOT the tran- 
sient effects become irrelevant. However, in this case the 
number of photons in the strong driving field Np grows 
proportionally to the cross section A, as evident from 
Eq. (61). When a photon number of the strong field is 



too high it becomes more difficult to implement protocols 
based on separating and photon counting of the quantum 
mode. Protocols based on homodyne detection also place 
limits on the maximal number of photons in the driving 
field. In most cases the driving field is also used as the lo- 
cal oscillator for homodyne detection. This implies that 
in order for the detection to be shot noise limited the 
classical fluctuations of the strong field should be sup- 

— 1/2 

pressed to better than N p . In practice this places the 
limit Np < 10 10 . For higher values of Np modulation 
techniques similar to those used with thermal ensembles 
allow to get down to shot noise limited detection. 

Another difficulty of working with a MOT is due to 
the presence of gradient magnetic fields and the associ- 
ated difficulty of optical pumping and magnetic states 
decoherence. This problem can be overcome by switch- 
ing off the MOT fields which in turn limits the lifetime 
of the interface, and by using a clock transition as was 



done in (Zhao et al. 2009a) where the memory lifetime 



up to 1 msec was achieved. An advantage of using a small 
sub-ensemble of a large MOT is that a single MOT can 
then serve as a source of two or more atomic ensembles 
Choi et al. 



the frequency ujl- In the experiments by Julsgaard et 



(Chen et al. 2007a 



Kuzmichl 2004). 



2008 Matsukevich and 
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Using a far detuned dipole trap allows to overcome a 
number of problems associated with the MOT. A dipole 
trap forms an atomic ensemble with a typical transverse 
size of around 10-50 jam which is a good size for the 
interface. The resonant optical depth can reach 20 or 
more. Dipole trapped atomic ense mbles demonstrate co- 
herence times exceeding 10 msec (IWindpassinger et al. 



2008) a nd storage times for single excitations exce eding 



stantial optical depth can be created, although its value 
is usually limited by the ion-ion interaction at high den- 
sity of doping. Note, however, that one cannot directly 
compare the measured optical depth in the presence of 
inhomogeneous broadening with the optical depth intro- 
duced in the theoretical derivation in Sec. [ill see sec - 



VII EIT- and Raman- based memory (Longdell et al. 



2005) has been explored in REIC materials as well as 



100 /as (Chuu et al. 2008) or even a few msec (Zhao a photon echo approach based on Controlled Reversible 



et al.< 2009b). Dipole trapping can also be insensitive Inhomogeneous Broadening (CRIB) (Hetet et al. , 2008a 



to the magnetic quantum number and, in some cases, 
to the hyperfine quantum number. A detailed investiga- 
tion of dipole trapped atoms as the medium for the spin 



squeezing and the int erface is given in ( |Oblak et a/.||2005 
Windpassinger et 



L) [2008j ). 

A Bose-Einstein condensate (BEC) is an attractive 
medium for the quantum interface due to its very high 
resonant optical depth. Indeed, BEC has been the 
medium where classical coherent storage of light, the 
so called stopped light has been first demonstrated (Liu 



et al., 2001). It is important to note that it is not only 



the on-axis optical depth that defines the strength of the 
interface coupling. For example, a "one dimensional" 
sample will not couple efficiently to a focused Gaussian 
beam because the diffraction of the beam means that it 
will not be overlapping with the ensemble if the Fresnel 



number is less than unity ( Miiller et al. t 2005 ) 



One problem with a BEC-based quantum interface is 
the low rate at which experiments can be performed. A 
typical BEC requires tens of seconds to be created. Then 
the sample can be used for interface experiments a few 
times, after which a new sample should be created. BEC 



on a chip ([Hansel et al. 2001[ Schneider et al. 2003 ) of- 



fers an attractive alternative where much faster loading 
times can be combined with very efficient optical cou- 
pling. 



C. Solid state 

An optically dense collection of atom-like impurities in 
a solid state host is an excellent candidate for the quan- 
tum interface. Both A schemes as well as the photon 
echo-based memory (see Sec. V.E) have been investi- 
gated. The absence of motion in solid state means that 



complex spatial structures can be generated by light and 
stored. Recently substantial progress has been achieved 
with crystals or glasses doped with rare-earth elements, 
such as erbium (Er), thulium (Tm), praseodymium (Pr), 
neodymium (Nd), and europium (Eu). The rare-earth 
ions doped into glass and crystal materials (REIC) dis- 
play up to a second coherence times of the ground state at 
liquid Helium temperatures. The ions experience strong 
inhomogeneous broadening up to 100 GHz due to local 
lattice fields. The technique of spectral hole burning 



( Nilsson et al. , 2004 ) followed by spectrally selective anti- 



hole populating allows to create a sub-ensemble of ions 
with nearly natural optical transition bandwidth. A sub- 



Kraus et at] [20061 [Moiseev and Kroll||2001||de Riedmat-| 
ten et al.[ |2008||Staudt et al.[ |2007| ) 



The optical Raman coupling to the nuclear spin coher- 
ence has been investigated for the past decades in REIC, 
mostly in praseodymium- or europium-doped crystals. 
These materials seem particularly suitable for the EIT- 
and Raman-based quantum memory protocols since they 
exhibit a hyperfine structure where a A-system can be 
found, together with long optical coherence lifetimes, and 
also long hyperfine coherence lifetimes (15ms for Eu:YSO 
and 550 /as, that can be extended up to 30s by dynamic 
decoherence control techniques, for Pr:YSO ( [Fraval et al] 
2005)). The absorption wavelength of these materials is 



in the domain of dye lasers (606 nm for Pr and 580nm for 
Eu) . In order to take advantage of their long optical co- 
herence lifetimes, the dye laser sources must be stabilized 
down to less than 1 kHz. The absorption wavelength of 
Tm lies in the convenient range of diode lasers (793 nm). 
It also exhibits long optical coherence lifetimes, similar to 
that of Pr. It has been recently shown that it is possible 
to build a A-system in thulium b y applying a magnetic 
field in a very specific orientation dLouchet et al. , 2008). 



In rare-earth ion-doped crystals, the transitions are not 
polarization-selective, so the only way to address them 
separately is to use a source whose bandwidth is smaller 
than the ground state sublevel splitting. In Pr and Eu 
the splittings are fixed (10 and 17 MHz for Pr:YSO, 75 
and 102 MHz for Eu:YSO) whereas in Tm they can be ad- 
justed with magnetic field (36MHz/T in Tm;YAG). The 
hyperfine coherence lifetime of up to 300/isec has been 
measured in Tm:YAG. 

Er-do ped materials are st udied with photon echo tech- 
niques (Staudt et al. 2007). The prime interest to this 



ion is due to the optical wavelength in the telecom band 
1.5 jam. The optical coherence lifetime of this material is 
very short (a few /isec at most), but can be dramatically 
increased by applying a very intense magnetic field. The 
most promising results up to date have been achieved in 



Pr-doped crystals (Hetet et al. 
20051). 



2008a Longdell et al. 



Materials containing a high concentration of quantum 
dots may be interesting candidates for the ensemble- 
based interface. Experiments with spin polarized dots 
show sufficient ground state coherence times and possibil- 
ity of optical pumping and quantum nondemolition cou- 
pling (Atature et al] 2007). However, up to now the work 



with ensembles of dots has not reached quantum limits 
probably due to an insufficient optical depth. A different 
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solid state medium which can be used for a quantum in- 
terface is Nitrogen vacancies in diamond, where EIT has 
been observed by 



ori entanglement. This observation lies at the heart of the 



Hemmer et 



(2001) 



original quantum repeater protocol (Duan et al. 2001), 



D. Other possible media 

Optical lattices have attracted a lot of attention lately 
due to the exciting possibilities for generation o f entan- 
glement by controlled atom-atom interaction ( |Mandel| 
The lattices can also display high opti- 



ca al. 2003) 



which is based on entanglement of atomic ensembles her- 
alded by detection of single photons, as will be described 
in Sec. HV.Cl Our main focus in Sec. IIV. 51 will be on the 
a priori entanglement achieved via a QND-Bell measure- 
ment and feedback on two ensembles. To introduce this 
method, we first explain how a single atomic ensemble 
can be prepared in a spin squeezed state by means of a 
QND interaction, homodyne detection of light and feed- 
back on atoms. Note that in Sec. IV. Dl we describe the 



cal depth since structures of up to 100 x 100 x 100 atoms 
spaced by hal f a micron can be created. In the recent 
experiment by Zhao et al. ( 2009b ) long memory lifetimes 



memory experiment (Choi et all 2008) which involves a 



exceeding 6 ms were achieved in a one dimensional opti- 
cal lattice. A quantum interface with such a lattice would 
offer an exciting possibility to transfer entanglement from 
atoms to light and to combine the quantum information 
processing capabilities of lattices with the quantum net- 
working provided by the quantum interface. Theoretical 
studies of quantum interfacing of light with lattices have 



heralded entanglement as an intermediate step. 



A. Spin-Squeezing in a Single Ensemble 

Spin sque ezed states of atomi c ensem bles were in- 
troduced by Kitagawa and Ueda (1993) in analogy to 



recen tly appeared ([Eckert et al. 



2008 



Muschik et al. 



squeezed states of the radiation field and suggested by 
Wineland et al. ( 1994 1992[ ) to be of use for enhancing 
the sensitivity in atomic spectroscopy, Ramsey interfer- 



20081) and the first EIT-based stoTaiTof classical light ,22^' and atomic clocks ' Accordingly, |Kitagawa and| 



with the coherence time of 240 msec has been demon- 



strated ( Schnorrberger et al., 2009). 



Another system where a collection of atoms can be 
efficiently interfaced with light is a large ion crystal in 
an ion trap. Very clean and large ion crystals have been 
created and first attempts towards achieving quantum 



Ueda| ([1993]) define the state of a collective spin J to 
be squeezed if the variance of one spin component J±_ 
transverse to its mean polarization is smaller than the 
transverse varian ce corresponding to an atomic coherent 
(Bloch) state (Arecchi et al. 1972[ ), that is a product 



coupl ing to light have been undertaken (Herskind et al. 

2008). 



IV. ENTANGLEMENT OF ATOMIC ENSEMBLES 

In this chapter we describe generation of entangled 
states of two distant macroscopic objects. The first 
method, which is based on QND interaction, mea- 
surement and feedback, generates an EPR (two-mode 
squeezed) state of atomic spins. The second method, 
which relies on parametric gain and beam splitter inter- 
actions and single photon counting, creates Bell states in 
two collective spins. 



van Enk et al. (2007) gave a useful classification of 



the various types of entanglement, which are generated 
in experiments. They distinguish a priori entanglement, 
which can be deterministically generated, a posteriori en- 
tanglement, which is generated probabilistically and de- 
stroyed when measured, and finally heralded entangle- 
ment, which is as well probabilistically generated, but 
success can be testified by measuring an auxiliary sys- 
tem, such that the entangled state is still available for 
use. Using post-selection of successful cases, all types 
of entanglement are in principle equally useful. When 
combining a large number of entangled states, e.g. via 
entanglement swapping in a quantum network, the over- 
all success probability will be dramatically different for a 
posteriori entanglement as compared to heralded or a pri- 



state of fully polarized atoms. With this definition, a 
state is squeezed if t;s = AJ j/ \J J/2 < 1 and necessarily 
consists of correlated atoms. Wineland et al. (1994, 1992) 
on the other hand show that the figure of merit for the 
suppression of quantum fluctuations, which ultimately 
limit the sensitivity of atomic Ramsey interferometry, is 
£r = (2J) 1 / 2 AJ±/\(J )| < 1 and provides an alternative, 
stronger definition of spin squeezing. 

We will here follow yet another definition and refer to 
a spin state as squeezed if £ = AJ±/(\(J )|/2) 1 / 2 < 1, 
which is stronger than the definition due to| Kitagaw a~and| 
Ueda] but weaker than the one due to | Wineland et al.\ be- 
cause £r > £ > £5 as can be easily seen. The conditions 
are the same for nearly fully polarized states (J) ~ J. If 
we take the mean polarization along x, and assume that 
the transverse component with minimal variance is along 
z, the definition of spin squeezing adopted here is 



APi< 



2' 



(68) 



where we use the Gaussian approximation (16). We use 



this definition, because it immediately translates into the 
entanglement criterion for a bipartite state of two ensem- 
bles in Sec. IrvTBl 

Various ways to create squeezed states in ensembles of 
two level systems were proposed . They involve direct in- 
teraction of spins (|Andre et al. , 2002 Pu and Meystre 



2000| |S0rensen et a/.[|2OOl||S0rensen and M0lmer| 



mapping of squeezed light onto atoms (|Appel et al. 
Dantan et al.| |2006a||Hald et a/.[[T999l|Honda et ~ 



1999) 



2008 



2008 
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Kuzmich et al. t 1997) 



atoms (Hammerer et al. 2004 Takeuchi et al. 



multiple passes of light throu gh 
2005k or 



a projective, Faraday interaction based QND measure- 
ment 



dBrj 

(19' 



Kuzmich 



Braginsky and Khalili 
et al.\ (119981 120001). The main idea in the last method is 



1996), as used by 



that light correlated with a collective atomic spin via a 
Faraday interaction can be used as a meter system, read- 
ing out one of the spin components. Homodyne measure- 
ment of light, as discussed in section |lLF} then provides 
information about this spin component, projecting the 
collective spin into a state with reduced fluctuations in 
this component. 

Mean values and variances of transverse spin compo- 
nents conditioned on a homodyne measurement of light 
can be easily evaluated by means of the following classical 
formulas for the mean value and variance of a Gaussian 
random variable £ conditioned on the measurement of 
another (possibly correlated) Gaussian random variable 
C with an outcome z, 



(0 



c= 



(0 



(JO 

<C 2 > 



A£ 2 



k= 



A£ 2 



(jcl 

(C 2 ) 



If we assume that both light and atoms are initially pre- 
pared in their vacuum states, i.e., atoms are completely 
polarized along x, then the states of atoms and light af- 
ter the Faraday interaction still have Gaussian statistics, 
as given by (40), and the above formulas apply. Suppose 



that Xl is measured on the state after the interaction 
with an outcome xl- According to the formulas above, 
the conditional variances of atomic spin components are 
then given by 



AX 



1 



A \x L 



AP 



1 



A \x L 



1 
?2' 



(69) 



and exhibit spin squeezing. 

Conditioned on the measurement outcome xl. Pa is 
squeezed with the mean value given by 



(Pa)\ 



Xl—xl 



1 



x L . 



If we ignore the measurement outcome the evolution is 
given by Eq. (40) and there is no squeezing since Pa 
is conserved. If a feedback operation is applied to the 
atoms, e.g., by applying a pulse of magnetic field, the 
atomic spin can be tilted such that Pa is displaced by 
—kxl/(1 + ft 2 ), and the mean value is zero (Pa) — 0. 



The variances (69) for the (anti) squeezed variances then 
hold also in the ensemble average. 

In Sec. III.FI we introduced another method for describ- 
ing the linear feedback and it is instructive to apply it 
here. The result xl is fed back on atoms by displacing 
Pa by an amount gxL, where g is a suitable gain factor. 
By Eq. (67) and (40), the Pa component after feedback 



is in the ensemble average given by, 

^A,final = PA,out + fl^X,out = (1 + ^)^A,in + 0^L,ii 



(70) 



Minimizing the variance of P4,fmai with respect to the 
gain g, yields an optimal feedback gain g = — k/(1 + n 2 ) 
and the reduced variance of Eq. ( 69 ) , in agreement with 
the discussion abov e. An exp eriment alo ng these lines 
has been reported in Kuzmich et al. ( |2000| ). Spin squeez- 
ing of 1.8 dB in the collective spin of a cold ensemble o f 
Ytterbium atoms was reported by|Takano et al. (2009). 



Recently spin squeezing on atomic clock transitions ha s 
been demonstrated for Cesium in Appel et al. (2009), 



wi th £ = — 4.5dB and £r = — 3.4dB, and for Rubidium 
([20091), with £ R 



Schleier-Smith et 



-3.2dB. 

The discussion so far ignores the impairing effects of 
spontaneous emission and light absorption. In order to 



take it into account we have to resort to Eq. (56). For 



the relevant case of small atomic decay, t\a <C 1, and 
dominant light losses due to reflection at glass cells and 
detector inefficiency parametrized by e (1 ^> e ^> tjl) 
these equations read 



XA,out — v 1 
PA,out — V 1 

^L,out = VT- 
P L ,out = VT- 



Va(X a , 



VAfx A 



VaPa 



<*l, 



eP L A 



n + ftP. A, in J 
+ VVAfp A , 

VtfPL, 



where we explicitly included Langevin noise operators for 
atoms fx A (P A ) an d light fx L (P L )- For both systems one 
can to a good approximation assume vacuum properties 
(fafp) — 5<x,p/2- Using these expressions and minimizing 
the variance with respect to the gain g yields a minimal 
variance 



AP 2 



1 + 7^(1-6)^ 

l + (l-e)^ 2 



1 > VA 

2 ~ 2 ' 



(71) 



The bound on the achievable squeezing is not so surpris- 
ing, given that the state of atoms suffered essentially a 
decay by a fraction tja- Due to the relation n 2 = dr]A, 
cf. Eq. (59), there is always an optimal choice for the de- 



cay t] a given a certain optical depth d (Hammerer et al. 



2004 ) . For the decoherence model adopted here the limit 



to spin squeezing by QND measurement and feedback 
is AP 2 final > (1 + Vl + d)- 1 - cT 1 / 2 , for large op- 
tical density. For a detailed discussion on the limits of 
spin squeezing by means of Faraday interaction and QND 
measurement we refer to Bouchoule and M0hner| (2OO2). 
The limits strongly depend on the particulars of the QND 
sche me. For example, for t he, so -called, two-color prob- 
ing ( Windpassinger et al, 2008) the 1/d scaling of the 
squeezing limit is achievable. 

Our description here covers only feedback where the in- 
tegrated photocurrent of the homodyne detection is taken 
as the measurement result and used to correct the atomic 
state after the probe pulse has passed. It is of course 
possible to perform a continuous feedback of the pho- 
tocurrent while the probe pulse is still on. An exhaustive 
theory for this procedure giving a description in terms 
of the stochastic Schrodinger equation can be found in 
Thomsen efal] ( |2002a|b ). 
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It is interesting to note that, beyond atomic interfer- 
ometry, spin squeezed states received renewed interest in 
the theory of many particle entanglement. It was shown 
S0rensen et at. (|2001) that spectroscopic squeezing 



£r < 1 is a sufficient condition for bipartite entangle- 
ment within each pair of spins in the atomic ensemble, 
see also |Wang and Sanders ( 2QQ3| ). More general spin 
squeezing inequalities were fruitfully studied in the con- 
text of experimental verification of multipartite entan- 
glement^ Ko£biczet_a/. 2QQ5a|b S0rensen and M0hner 



2001 Toth et all 2007) 



B. Deterministic entanglement 

The procedure discussed in the previous section can be 
applied to deterministically create entanglement between 
tw o, spatially separated atomic ensembles, as proposed 
by Duan et al (2000a). A precursor to this proposal 



inv olving entangled light resou rce has been put forward 
by Kuzmich and Polzik (2000). Each of the ensembles 
is described by a collective spin Ji (i = 1,2), or, taking 
polarizations along x for both systems and adopting the 
Gaussian approximation, by a pair of canonical opera- 
tors [xauPa ■} — i5i,j- Consider a probe pulse undergoing 
Faraday interaction with both ensembles, first with en- 
semble 1 and then, after propagating some distance, with 
ensemble 2. By linearity, the state of light is desribed by 



Pl, out 






■«(-FUi,in +^A 2 ,in), 



c.f. Eq. (40). Note that for both ensembles P Ai out 



-Pa*, in is a conserved quantity in the Faraday interaction. 
Just as for the single ensemble, a measurement of X^ 50nt 
will then give a reduced variance of the non-local observ- 
able P Al + P A2 , 



A(P Ai +Pa 2 ) 2 



1 



<1, 



where the bound corresponds to uncorrelated ensembles 
in coherent states. Using feedback this non-local squeez- 
ing can be achieved unconditionally. In a second step, 
the spins of both ensembles are rotated by an angle of 
7r/2 about the x-axis in such a way that 



Pa x — >• X Al , 
Pa 2 -> —Xa 2 , 



Xa x — » — Pah 

X A2 ^Pa 2 . 



A second light pulse interacting with both ensembles as 
before, will then read out the observable XA lj i n — ^A 2 ,in> 



XL,out — ^L,in + K>(Xa 1 ,v 



Xa 2 ,iyi), 



such that a measurement and feedback procedure will 
produce a squeezed variance of Xa x — Xa 2 , as before. 
Note that simultaneous squeezing of these two observ- 
ables is possible only, because we are dealing here with 



commuting observables, [P Al +Pa 2 ,Xa 1 —Xa 2 ] = 0. The 
counterwise rotation of the two spins about x is therefore 
crucial. Overall, this will produce a state which fulfills 
the inequality, 

A(P Al + P A2 ) 2 + A(X Al - X A2 ) 2 < 2. (72) 

Losses and decoherence will affect this result similarly 
as the single ensembles spin squeezing in (71). The sig- 
nificance of this inequality is that it constitutes a neces- 
sary and sufficient entanglement condition for symmetric 
(with respect to 1 j-> 2) Gaussian states of two systems, 
flDuan et a^2000b||Simon[|2000| ). 

In the limit of large squeezing, where the variances 
of both non-local observables vanish, the corresponding 
state approaches the (unphysical) ideally correlated state 
with Wigner function W ~ S(X Al - X A2 )S(P Al + Pa 2 ), 
which was considered by Einstein, Podolsky and Rosen 
(EPR ) in their famous Gedankenexperiment (Einstein| 



et al> 1935), speculating about the incompleteness of 
quantum mechanics. See Keyl et al. (2003) for com- 



ments on whether and how this limit can be understood 
in a more rigorous mathematical sense. Because of this 
connection, the quantity on the left-hand side of Eq. 
( 72 ) is sometimes termed EPR- Variance and denoted by 



Aepr- Its importance is supported by the fact that for 
symmetric Gaussian states the quantity r = — - In Ae 2 pr 
provides an entanglement measure and uniquel y deter- 
mines the entangl ement of formation of the st a te ([Gi edke 
et all 120031 iWolf et all 120041 iWoottersl 119981) via 
E p = cosh (r)log 2 (cosh r) — sinh (r)log 2 (sinh r). As 
follows from the discussion of losses and decoherences, 
we have to assume a lower limit on the EPR variance 
Aepr ^ 2d~ 1 / 2 and thus an upper bound on the bi- 
partite entanglement between the two ensembles of e.g. 
^oF ^ 1.15ebits (£?oF ^ 2.77ebits) for an optical density 
of d = 10 (d — 100). A useful and comprehensive review 
of the theory of entanglement in systems of continuous 



variables was recently given by Adesso and Illuminati 



(2007), for a concise introduction to the basic facts on 



the same topic see Eisert and Plenio (2003). 



Protocol with counter-rotating spins As explained in 
Sec. II. E[ for ensembles at room temperature containing 
a very large number of atoms a constant magnetic field 
helps to reduce technical noise, because scattered light 
can be detected at sideband Zeeman frequencies. In the 
following we will show, that application of an external 
magnetic field provides in fact an elegant and efficient 
way to achieve entanglement of two atomic ensembles 



with a single prob e pulse, as was demonstrated in Juls- 
gaard et a/.|(|200lD. 



In order to show this, we have to resort to the Maxwell- 
Bloch equations ( [65] ). It is straightforward to generalize 
these equations to the case of two atomic ensembles and 
to include Larmor precession of the two spins. We as- 
sume that the two collective spins precess in opposite 
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directions, which can be achieved by either using op- 
positely oriented fields or by using parallel fields with 
atoms prepared in opposite Zeeman substates. Replacing 
ujl — >■ —ml for the second ensemble when using Eqs. (65 ), 
the equations of motion for light quadratures are, 



a_ 

dz 



x L (z,t) 






s(u L t) (p Al (z, t) + p Al (z, t)) 



+ sm(u L t) (x Al (z, t) - x A2 0, t)) 



_9 
dz 



PL(z,t) = 0. 



(73) 



Obviously light reads out sums of momenta and differ- 
ences of positions, which are commuting observables for 
oppositely rotating spins. It will be instructive to study 
directly the evolution of these global observables. One 
finds, 

9 / , x , \\ 9*(z)Q(t) / x / x 

— {x Al (z,t)+X A2 (z,t)) = C08(u> L t)p L (z,t), 

d_ 
dt 
d_ 
dt 

d_ 

dt 



(p Al (z,t)+p A2 (z,t)) =0, 
(x Al (z,t) -x A2 (z,t)) =0, 

(pAt(z,t) -pA 2 (z,t)) 



g*(z)n(t) 

V"2A 



s'm(uj L t)p L (z,t). 

(74) 



The solutions fall into two groups, 



XL c ,out — -^Lcjin + ttPA+jinj -^L s ,out — X 



-* L c ,out — *L C 



Pl s , out — -Pl. 



L s ,in ^^-A_,in? 
in 7 



^A + ,out — ^A + ,in + ^^L c 
■* A+,out — -TA + ,iri) 



-^A_,out — X A _ 
Pa_, out — Pa_ 



"«ft- 



(75) 



which involve non-local atomic variables 
1 



^4± 
^4 ± 



2L 
1 



^(x Al (^) ±xa 2 (^)), 
dz(p Al ( z ) ± PA 2 (z)) 



and cosine and sine modulation modes, Xl c , Pe,,, and 



Xl s , Pl s , which were introduced in (66) in section II. F 



From the discussion of squeezing in a single ensemble it is 
evident that a measurement of the sine and cosine com- 
ponent of the Xl quadrature will produce a two mode 
squeezed state with reduced EPR variance (72). 



Implementation Experimental demonstration of deter- 
ministic entanglement of two atomic ensembles has been 
first reported by Julsgaard et al. (2QQ1[), with further 



developments reported by Polzik et al. (2003); Sherson 



with Cesium vapor at temperatures in the range 15° - 
50° C. Atoms are contained in glass cells coated from 



inside with a transparent layer of paraffin (Alexandrov 



2003), as discussed in Sec. III. A The temperature stabi 



lized cells are placed inside cylindrical magnetic shields, 
see Fig. [9j Windows made in the shields allow for optical 
axis in two directions - one along the axis of the shield 
used for optical pumping and another in the radial di- 
rection used for the probe light. A solenoid produces 
a homogeneous axial magnetic field inside each shield. 
Typical cells have the near-cubic shape with the size of 
25 — 30 mm. The magnetic field inhomogeneity is of the 
order of 10 -3 . This rather modest homogeneity is suffi- 
cient since the duration of the light-atoms interaction of 
1 msec is sufficiently long so that the atomic motion leads 
to the effective time averaging of the spatially dependent 
Zeeman shifts. As a result the magnetic field inhomo- 
geneity has only a quadratic effect on decoherence. 

The experimental sequence begins with a few msec 
pulse of optical pumping along the direction of the axial 
magnetic field, see Fig. [9] The level scheme and fre- 
quencies of light pulses are shown in Figs. [3] and |4j The 
two cells are pumped with the same lasers with opposite 
circular polarization of optical pumping in the first and 
second cell. 99 % or more of the atoms in F = 4 state 
are pumped into the mp = 4 magnetic substate in one 
cell and into mp = —4 in the other cell, as verified by 



the magneto-optical resonance method (Julsgaard et al. 



et al. (2006b). The experiments have been performed 



2004b Sherson et al. , 2006b[ ) . The total angular momen- 
tum J x of F = 4 state is calibrated by measuring the 
Faraday rotation angle of a weak linearly polarized light 
pulse propagating in the direction of the optical pump- 
ing. After optical pumping a probe pulse linearly polar- 
ized in x — y plane is fired and its polarization rotation, 
i.e., the value of the operator S y , is measured by two 
detectors via a balanced polarization measurement, see 
Fig. [5] The detected photocurrent pulse is sent into the 
lock-in amplifier which detects the cos(u)Lt) and sin(cj^^) 
components, Xl c and Xl s introduced in Sec. II. F| 

The critical condition for the implementation of the 
deterministic interface based on homodyne mesaure- 
ments is quantum noise limited measurement of light 
and atoms. This is made possible by employing light 
and atomic detection at high frequency, typically around 
ujl = 320 kHz. Homodyne detectors utilize silicon pho- 
todiodes with quantum efficiency more than 98 — 99 % 
and low noise photo-amplifiers with the response peaked 
around uo l- The photodetectors have dark noise equal to 
shot noise of light for the light power as low as 100 /iW. 
This means that with a few mW probe pulse, the detec- 
tion can be almost perfectly shot noise limited. Light 
losses have been dominated by reflection off the inner 
surfaces of the cell windows (the outside surfaces are an- 
tireflection coated) amounting to 15% and propagation 
losses from cells to detectors of about 8%. 

The duration of the interaction is chosen to fulfill the 
condition of optimal entanglement which, analogous to 
the case of spin squeezing, is k 1 t = \[d obtained with 
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FIG. 9 Deterministic entanglement of two atomic ensembles via QND measurement, a) Experimental setup. Dashed lines 
- optical pumping, solid arrow - entangling light direction, b) Pulse sequence and the layout of the experiment, (l)-optical 
pumping pulse, (2)-entangling pulse, (3)-verifying pulse, c) Projection noise of atoms d) EPR variance of the entangled state 
normalized to the projection noise level, i.e., to the variance for a separable coherent state of two ensembles. 



rjA = l/\d (cf. Eq. (59)). In the experiment the detun- 
ing has been chosen within the range 800 — 1000 MHz 
to be larger than the Doppler width and the hyperfine 
splitting of the excited state. Together with the optimal 
value of the optical power set by the detectors around a 
few mW, the above conditions lead to the minimal pulse 
duration of the order of a msec and the corresponding 
bandwidth in the kHz range. This value can be in prin- 
ciple increased by reducing the transverse size of the sam- 
ple and/or using different detectors. 

The first experimental run in the presence of atoms 
aims at the establishment of the projection noise level of 
the atomic ensembles. As seen from Eq. ( [75] ) the sum 
variables Xa_ , Pa+ can be measured in a QND way us- 
ing a single probe pulse. In the experiment the sequence 
of the optical pumping followed by the QND measure- 
ment is repeated several thousand times and the variance 
of the measured photocurrent pulses is calculated. The 
variance is then plotted as a function of the macroscopic 
spin of the sample, see Fig. [9Fc). The linear dependence 
along with the almost perfect spin polarization proves 
that the spin noise is at the projection noise level. The 
projection noise level has been also independently calcu- 
lated from the macroscopic collective angular momentum 
of the sample measured via Faraday rotation of the aux- 



iliary probe pulse (iSherson et al. , 2006b). The calculated 



value agreed with the measured projection noise to within 
10 % which is well within the uncertainty of this calcu- 
lation. The projection noise level defines the right hand 
side in Eq. (72). When normalized to the shot noise of 



the probe the projection noise value is equal to the total 
ft 2 of the two samples according to ( 75 ) . 



After the projection noise level is established the ex- 
periment proceeds with generation and verification of the 



entangled state. In the original paper (Julsgaard et al. 



2001) no feedback was applied to atoms and hence a 



conditional entangled state was demonstrated, that is a 
non-local state with reduced variance but with a non- 
deterministic mean value. For possible applications, e.g., 
teleportation this entanglement is as good as the uncon- 
ditional one because the knowledge gained with the mea- 
surement on the first entangling pulse can be applied to 
achieve teleportation. The creation of a deterministically 
and uncondit ional entang l ed sta te has been achieved sub- 
sequently ( Polzik et al.\ 2003). The ex perimental se- 
quence which realizes such entanglement (|Sherson et al. 



2006b) involves a feedback applied to the atoms in be- 
tween the two probes (Fig. ^p). The feedback pulse of 
320 kHz magnetic field with the cosine and sine compo- 
nents proportional to XL c ,out an d ^L c ,out respectively is 
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applied to the rf magnetic coils surrounding the cells. An 
appropriate electronic gain must be chosen so that the 
feedback pulse rotates the atomic collective spins such as 
to generate the minimal EPR variance that is the mini- 
mal variance of the angle between the spins. The choice 
of the gain g = — k/(1 + ft 2 ) minimizes the EPR vari- 
ance in the absence of decoherence. In the experiment 
the optimal gain has been cho sen oper ationally by mini- 
mizing the EPR variance (see Sherson et al. ( 2006b ) for 



details). The results for the EPR entangled state of two 
atomic ensembles are shown in Fig. [9j The variance of 
the entangled state obtained after applying the feedback 
pulse is measured by the verifying pulse. Fig. [9] shows 
this variance normalized to the projection noise variance 
as a function of the number of atoms. A higher num- 
ber of atoms leads to a higher value of n and hence to a 
higher degree of entanglement. The mimimal EPR vari- 
ance observed in these experiments was Aepr = 1.3. 
This variance corresponds to an entanglement of forma- 
tion of E q f = 0.28 ebits. 



C. Probabilistic entanglement 

One of the main motivations for studies of atom light 
quantum interfaces is an application for quantum re- 
peaters, which would enable long distance quantum com- 



munication (Briegel et al. 1998). A protocol (also known 



as the DLCZ-protocol) for such a repeater based on 
atomic ensembles and linear optics was first presented by 
Duan et al. (2001). Several improvements of the proto- 



col have been suggested ( Chen et al. |2007b Jiang et al 
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Zhao et all 2007). We leave the detailed discussion of the 



DLCZ-protocol and quantum repeaters for a dedicated 
review. In this section we provide a basic discussion of 
the probabilistic entanglement generation. 

The entanglement generation in the D LCZ -protocol 
uses the parametric gain interaction Eq. (112) between 
the light and the atoms as shown in Fig. lOfr . In prin- 
ciple this interaction could be used in the strong cou- 
pling regime (k, ~ 1 in the notation below) to generate 
continuous variable entanglement along the lines of the 
entanglement generation protocol used for quantum tele- 



portation in Sec. VI Instead the DLCZ-protocol works 
in the weak coupling limit (k<1) and gene rates proba- 
bilistic entanglement. The first term in Eq. (112) for a a 
is just the AC-stark shift of the ground state, which can 
be removed by a simple rescaling, and the phase <\> van- 
ishes for a large detuning. In the limit of weak coupling 
the dynamics only involves collective operators analogous 
to the ones defined in Eq. (55) and is equivalent to the 



evolution with the ideal two mode parametric gain evo- 
lution operator exp(m(a^a^ + H.C.)/2). For k<1 the 
joint state of the collective atomic and light harmonic 
oscillator degrees of freedom is then 




FIG. 10 DLCZ-protocol. a) Entanglement between light and 
atoms is generated by a parametric gain interaction in two dis- 
tant ensembles. Detection of a photon by one of the two pho- 
todetectors after the optical beam splitter probabilistically 
generates entanglement between the two ensembles, b) The 
atomic excitation in one half of an entangled pair of atomic 
ensembles (dashed line) is read out onto light with the light- 
atoms beam splitter interaction and mixed on an optical beam 
splitter with its counterpart from another entangled pair of 
ensembles. Photo detection swaps the entanglement so that 
now the outmost ensembles become entangled. 



where \mn) al describes the state with m (n) excitations 
in the atomic (light) harmonic oscillator. This result can 
be understood rather intuitively from the level scheme 
in Fig. 



10 1, which shows that the interaction generates 



simultaneous excitations of the atoms and the light as 
described by the expression above. 

The state in Eq. (76) is in itself an entangled state, 



but in the DLCZ protocol it is used to probabilistically 
generate an entangled stated of two atomic ensembles. 
The outgoing light modes from two different ensembles 
are combined on a beamsplitter as shown in Fig. [To] 
Conditioned on a click in one of the two photodetectors 
the ensembles are prepared in the state (|01) ± |10))/-\/2, 
because one cannot know from which ensemble the pho- 
ton was emitted (the sign in the superposition depends 
on which detector detected the photon) . If two pairs of 
ensembles are independently entangled in this fashion, as 
shown in Fig. [lOfr , the entanglement can be extended to a 
larger distance. Towards this end the atomic states from 
the two closest atomic ensembles, each belonging to a dif- 
ferent entangled pair, must be read out onto light modes 
which are mixed on another beam splitter. Detection of a 
photon after the beamsplitter extends the entanglement 
to twice the distance by entanglement swapping. This 
read out process can for instance be done using the beam 
splitter interaction as discussed in V.D Note, however, 
that the atomic mode functions ([29| and ( 32 ) suitable for 
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the parametric gain and beam splitter interactions are 
different leading to a mode mismatch if the excitations 
are read out in the same direction as the entanglement 
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generation. This problem can be avoide d by reading out 
in the backward direction ( Andre , [2005 ) (this analysis as- 
sumes that the atoms retain their positions throughout 
the experiment; see also |Duan et al. (p002) for a discus- 
sion of the limit, where the motion of the atoms leads to 
an averaging over the atomic position). 

Following the initial DLCZ-protocol several experi- 
ments have been performed, which demonstrate a num- 
ber of important ingredients of the repeater, of which 
here we only mention a few. The first experiments 



Chou et all 2004 Kuzmich et al p 2003 Matsukevich and 



Kuzmichl 20041 van der Wal et al, 2003) demonstrated 



quantum correlations in pairs of photons generated via 
the creation of pairs of atomic excitations and photons 
in Raman scattering, followed by a beam-splitter interac- 
tion converting the atomic excitations into another pho- 
ton after a programmable delay. Chaneliere et al ( 2005 ) 



and Eisaman et al (2005) reported the storage and con- 



trolled release of single photons, which are deter minis- 
tically created from a first atomic ensemble, in a sec- 
ond one via elec tromagnetica lly induced transparency. 
Chaneliere et ql.|fl2007D;|Chou et a/.| ( |2005| ); |Minto et al 
(2006); Yuan et al (2007) demonstrated the quantum 



interference of photons emitted from two independent 



atomic ensembles in Raman scattering, with ( Chou et al 



2005) proving entanglement of the two atomic ensembles. 
Chou et al ( |2007[ ) implemented an elementary link for 
a DLCZ-type quantum repeater as depicted in Fig. [10] 
Other steps t owards the DLCZ quantum repeater were 
performed by Chen et al (2008), where the teleportation 



of photonic to atomic qubits he ld in an atomic en semble 
has been demonstrated, and by|Choi et al (2008) where 



storage and release of photonic entanglement from two 
atomic ensembles has been reported. 



V. QUANTUM MEMORY FOR LIGHT 

Atomic quantum memory for light is an important in- 
gredient for a number of quantum information routines. 
It is implicit in many quantum communication proto- 
cols, in particular in those which require local operations 
on more than just a single optical pulse, so that stor- 
age is necessary. It is required for linear optics quantum 
computing and for scalable cluster state quantum com- 
puting with photons. Quantum memory is a necessary 
ingredient of a quantum repeater. Different applications 
demand quantum memories fulfilling different require- 
ments. Some applications, such as those which include 
local operations and classical communication (LOCC), 
benefit from a high-fidelity deterministic write-in opera- 
tion into the memory. By deterministic we mean a proto- 
col which works with probability one, so that the fidelity 
is calculated for every try. Others can tolerate lower 
probability of write-in and read-out, provided the suc- 
cess is heralded. Nonetheless, the latter protocols, such 
as a quantum repeater, would also benefit from deter- 
ministic write-in which would lead to a higher efficiency, 



higher rate, and eventually to longer distances. 

In this chapter we review several main approaches to 
the quantum memory for light. We first discuss a fig- 
ure of merit and a classical benchmark for determining 
the quality of a quantum memory. We present the pro- 
tocol based on a QND interaction and feedback, which 
demonstrates a quantum memory channel with the fi- 
delity higher than the classical benchmark. We then dis- 
cuss the memory based on the A scheme, concentrating 
on the recent achievement of the EIT-based memory ex- 
periments. We conclude with a discussion of memories 
based on various types of the photon echo. 



A. Figure of Merit 

From a fundamental perspective a quantum memory 
can be analyzed as a quantum channel, acting in time. 
A perfect quantum memory is nothing but the identity 
map taking arbitrary states as input and returning them 
unchanged, some time later. A realistic memory will be 
imperfect and the question arises what figure of merit to 
use in order to evaluate its performance. One sensible 
measure characterizing the performance of a memory is 
the fidelity, i.e. the average state overlap (Nielsen and 



Chuang 2000), which can be achieved between an in- 



put state drawn from a predefined set of states accord- 
ing to a predefined probability distribution and the state 
which is finally read out of the memory. The fidelity is 
of fundamental relevance if it exceeds the best classical 
fidelity, in which case the channel is thus outperform- 
ing the best classical channel. The classical fidelity relies 
on the simple strategy of measuring the given quantum 
state, storing the resulting classical data and on demand 
reconstructing the quantum state as good as possible. 

For example, in a special case where both input 
and output states are Gaussian states with amplitudes 

V^in(out))j Vin(out) 

the fidelity is given by 
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(77) 



In experiments the average of this fidelity can be taken 
with respect to a Gaussian distribution of coherent states 
centered at vacuum with a mean occupation number, n. 
If the average fidelity for a flat (infinitely broad) distribu- 
tion of coherent input states equals unity, the memory is 
ideal and would store also any Non-Gaussian state per- 
fectly. This follows from the fact that coherent states 
provide an (overcomplete) basis for the Hilbert space. 

It is worth emphasizing that the performance of a 
quantum memory can be, in principle, tested with coher- 
ent states only. Knowing the performance of the mem- 
ory (a quantum channel) for all coherent states, that is 
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performing a quantum tomography of the memory pro- 
cess with coherent states, one can predict the fidelity of 
the memory channel for an arbitrary class belonging to 
a single mode. In this sense the often used division be- 
tween continuous variable memory and discrete variable 
memory is not quite justifiable. It is more appropriate 
to speak about the protocols which are based on contin- 
uous variable measurements (homodyning) and discrete 
variable measurements (photon counting). 

The question then remains what exactly a measured 
average fidelity smaller than unity guarantees. The 
benchmark maximum fidelity of a classical channel F c i ass 
is known for a limited number of quantum states includ- 
ing qubit states and coherent states with a Gaussian dis- 
tribution in phase space of width n. For the latter case 



* 



Braunstein et a/.| ( |2000[ ) conjectured and |Hammerer et ■ 
(|2005b|) proved that 



class 



1 



-n 1 



(78) 



In the former case, for a class of arbitrary qubit states 
the maximum clas sical fidelity is F c i ass = 2/3 ( |Massar| 



and Popescu 1995). Very recently a classical benchmark 



fidelity for a third class of states, namely the displaced 



squeezed states, has been found (Adesso and Chiribella 



2008 Owari et al. , 2008 ) . For a class of pure squeezed 



states with the variance s in vacuum units, arbitrary ori- 
entation and arbitrary displacements the classical bench- 
mark approaches zero for large s as F c i ass = ^JL. Quan- 
tum memory which exceeds these classical benchmark 
fidelities thus shows performance which is classically im- 
possible. For similar fidelity benchmarks for finite dimen- 
sional systems see the work by Keyl and Werner (1999). 



The fidelity is not necessarily the one and only figure of 
merit for a quantum memory protocol. A relatively high 
fidelity may still be compatible with errors which are hard 
to correct. On the other hand, a lower fidelity protocol 
with particular kinds of errors may be more appropri- 
at e for a specific appli c ation. For instance the analysis 
in Brask and S0rensen (2008) shows that different types 



of errors can have a very different effect on the repeater 
protocol of Duan et al. (2001), and Surmacz et al. (2006) 



argue that in certain applications it might be more im- 
portant to preserve entanglement, when one partner of an 
entangled pair is stored, than to conserve the quantum 
state itself. 

For an important class of protocols discussed in Sees. 
|V.D| and |V.E[ the performance is theoretically described 
by a simple beam splitter relation between the input and 
output operators a out = y^^in + -\A — ^ where n is the 
efficiency, e.g, for mapping the input light intensity to the 
output light intensity, and v is a vacuum operator. The 
memory performance is then completely characterized by 
the single parameter 77, and quantities such as the fidelity 
for a given distribution of states may later be derived 
from it. The performance is therefore often discussed 
in terms of the single parameter 77 and we will use this 
characterization in Sees. IV.DI and lV.EI In an assessment 



of a given experiment one should, however, verify that 
the simple beam splitter relation is indeed applicable for 
this experiment. 

It is possible to define other figures of merit and also 
to consider different benchmarks than the one given by a 
classical measure and prepare strategy, in order to quan- 
tify the quality of storage - in quantum memories - but 
also of transmission of quantum states - as in quantum 
teleportation, cf. Sec. VI For experiments with single 
photons, it is common to consider the conditional fidelity, 
which characterizes the fidelity conditioned on the detec- 
tion of a photon after the interface. Since this condi- 
tioning suppresses the effect of losses, the conditional fi- 
delity is often much higher than the unconditional fidelity 
discussed above. When a conditional fidelity is used, 
another parameter, often called efficiency is introduced 
which describes the probability of success of the proto- 
col. In the context of teleportation of coherent states, 
the question of the "right" figure of merit was subject 



to cons i derable debate in the lite r ature, see|Bowen et al. 

f )03b|) ; |Braunstein et al. ( 2001 ); Grosshans and Grang- 
(2001); Ralph and Lam (|1998 ) and references therein. 
In particular, Grosshans and Grangier ( 2001| ) emphasize 
the importance of a benchmark Fi_^2 given by the max- 
imal fidelity achievable in a 1 to 2 cloning machine. For 
F > F\^2 the memory output is guaranteed to be the 
best possible copy of the input state. For c ohere nt in- 
put Fi^ 2 — 0.68, as shown by (Cerf et al. 2005), and 



is thus more demanding than the classical fidelity bench- 
mark Fdass , which can also be interpreted as the maximal 



fidelity of a 1 to 00 cloning machine (Hammerer et al. 



|2005b[) . Figures of merit different from fidelity were used 
Bowen efoT] d2003a|b| >; |Hetet et al\ ( |2008bD ; [Ralph and 



Lam 



(1998) to characterize both quantum storage and 
teleportation. There is, however, a consensus that the 
most sensible figure of merit ultimately depends on the 
specific application of the quantum memory or telepor- 
tation link within a quantum network for e.g. quantum 
cryptography or optical quantum computation. 



B. QND & Feedback Protocol 



The first demonstration of a quantum memory flJuls- 



gaard et al.\ 2004a) beating a classical benchmark (Ham- 



merer et al. t 2005b ) was based on the QND-Faraday in- 
teraction of a pulse of light, carrying the quantum state 
to be stored, with two collective spins counter-rotating in 
an external magnetic field. The basic input-output equa- 
tions describing the interaction are given by ( 75 ). Each of 



these constitutes a realization of the simpler input-output 



relations (40) for a pulse interacting with a single atomic 
ensemble without the magnetic field. For simplicity, we 
will base the theoretical discussion of the main idea on 
this single ensemble setup, and will return to the actual 
implementation based on the setup involving two atomic 
ensembles in the experimental part. 

The input light is described by canonical operators 
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Xk,in and Pl, in while the collective spin is prepared in 
the fully polarized state (XA,in) = (Pa, in) = and 
1/2. With a choice of ft = 1 in 



AXi >in 



APl, in 



(40) the entangled state of atoms and light after the in- 



teraction is described by 



^L,out 
XA,out 



Xl,\ 
Xaa 



Pa, 
Pl, 



Pl, out 
-Pa, out 



Pl;, 
Pa, 



Following the interaction the light quadrature Xc, >0 ut is 
measured and the corresponding measurement result £ is 
fed back displacing the atomic state such that Pa, out —> 



PA,out — £■ As shown in section II. F , the final transforma- 



tion of the collective atomic spin in the ensemble average 
is given by 



-^A,out 

Pa, out 



Xa,u 
PAAr 



-Pl^ 

^L,out 



-*la 



(79) 



This concludes the mapping of the quantum state of light 
onto atoms: the mean values are transmitted faithfully 
(apart from an unimportant phase change) as (X^ out ) = 
(PL,in) and (PA,out) = -(X L ; m ). The operator X L;in is 
mapped perfectly onto the atomic collective spin, while 
the operator Pl, in is mapped with the addition of one 
unit of vacuum operator which comes from the initial 
coherent state of the atomic ensemble. This latter im- 
perfection can be remedied if the initial atomic spin 
state is squeezed before the memory operation, such that 
AJ^ in — >• 0. If such squeezing operation is performed 
by, for example, an additional QND measurement, the 
fidelity of the quantum memory operation can, in princi- 
ple, approach 100%. 

In the experiment the quantum memory performance 
has been tested with a set of coherent states of light taken 
from a Gaussian distribution of coherent states centered 
at vacuum with a mean photon number, n. Given the 
measured gains and the measured variances, AX^ )OUt 
and APA jOUt of the state of the memory, the fidelity can 
be calculated as: 



P=(n(l-ft) 2 + l/2 + AXl ( 



-1/2 



x(n(l-£) 2 + l/2 + APi 50Ut ) 



-1/2 



As discussed above, if the protocol starts with the atomic 
ensemble in a coherent state A J^ out = 1, APj[ out = 

1/2, and hence F — y/2/3 « 82% with the choice of n = 
g = 1 which is optimal for the class of arbitrary coherent 
states. For an unknown qubit state (a|0) + f3\l)/V2 the 
same protocol yields the fidelity of 80% for optimal values 
of g and n. 

The experimental setup and the sequence of operation 
for writing into the quantum memory are similar to the 
sequence described in the chapter on deterministic en- 
tanglement. The quantum light mode which corresponds 
to the ujl sidebands in the polarization orthogonal to 
the strong field now carries the quantum state of light 
to be mapped. The state is generated by an electro- 
optical modulator (EOM) as shown in Fig. 11 The 



quantum sidebands together with the strong pulse prop- 
agate through the atomic cells and are analyzed with the 
polarization homodyning technique (Fig. [5]). The strong 
pulse thus serves a dual purpose, first as a strong driv- 
ing pulse for the interaction with atoms, and second as a 
local oscillator for the homodyne measurement. 

In the experiment two cells in a magnetic field play the 
role of one quantum memory unit. As discussed above 
this approach allows to achieve quantum limited noise for 
ensembles of trillions of atoms because quantum infor- 
mation is encoded and processed at ujl = 320 kHz side- 
bands where classical noise can be strongly suppressed. 
For two cells with magnetic field according to Eq. (75) 

the light variable Xl c - 



|t J dtcos(uLt)xL(t) should 



be measured and the result fed back to the atomic vari- 
able Pa + , out- The method for measuring Xl c by the ho- 
modyne measurement of the Stokes parameter S z of the 
light and the subsequent processing of the photocurrent 



by the loc k-in a mplifier indicated in Fig. [TT] is discussed 
" Eq. 



in section 



II.F 



Note that at the same time 



the variable Xl s — J ^ J dtsm(u)Lt)xL(t) can be mea- 
sured and fed back into the XA_,out variable of atoms. 
The memory could hence be used as a two mode memory, 
although this direction has not been pursued. 

After the projection noise level of atoms is established, 
as described in the section on entanglement, the optimal 
feedback gain must be determined. The gain is chosen to 
optimize the fidelity for the class of states to be stored 
in the memory. An example corresponding to the class 
of coherent states distributed around vacuum with n = 8 



is shown in Fig. 12 The optimal values for this class of 
states are n = 0.8, g = 0.8. After the memory sequence, 
cf. 



Fig. 
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is over the atomic variable Pa+, out is mea- 
sured witlTa strong QND verifying pulse. The quantum 
sideband modes of this pulse are initially in the vacuum 
state. After propagating through the memory the quan- 
tum sideband modes contain the atomic memory vari- 
able Pa + , out according to (75). 



Fig. 



12 



shows that the 



mean values of this variable for variousnght input states 
are very satisfactory proportional to the mean values of 
the input light canonical variables. As also shown in 
the figure the same is true for the other canonical vari- 
able of light Pl c stored in the memory variable XA + ,out- 
-^A + ,out has been measured in another experimental se- 
quence, where an RF pulse rotating the collective atomic 
spins by 7r/2 and thus converting Pa + , out into Xa , , ou t 
has been applied just before the verifying pulse (Fig.jnJ. 
The memory is thus shown to work very well as a classi- 
cal memory for light since the mean amplitude and phase 
of the input light pulse and the retrieved light pulse are 
equal to within a chosen factor. Note that as a clas- 
sical coherent memory this Faraday+feedback memory 
can have unity retrieval efficiency because the gain is ad- 
justable, whereas for the Raman, EIT, and photon echo- 
based memories discussed below the efficiency is less than 
unity in the presence of losses. 

The decisive demonstration of the quantum charac- 
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FIG. 11 QND+feedback memory experimental setup and pulse sequence ( Julsgaard et al. 2004a). The state of light is encoded 
by the EOM in the sidebands of the strong pulse yjn(t). Two cells serve as the quantum memory unit. The feedback pulse 
proportional to the photodetector signal is applied to the RF magnetic coils. Inset. (1) - optical pumping, (2) - input pulse, 
(3) - feedback RF pulse, (4) - RF pulse rotating atomic P into X used for half of verification pulses, (5) - verifying pulse 
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FIG. 1 2 Atomic coherent memory results (Julsg aard et al] 
|2004a| ). The mean values for both quadratures of the input 
light, atomic memory and the output light are identical to 
within a chosen factor (here 0.8). 



FIG. 13 Variance of the atomic memory state for input co- 
herent states with n < 8. The variance of 1/2 corresponds to 
the memory variance equal to the input light variance. Above 
the 3/2 level is the classical memory performance for arbitrary 
coherent states. The dashed line is the best classical perfor- 
mance for the states within the n < 8 class. Diamonds and 
squares are experimental results from Julsgaard et al. (2004a). 



C. Multipass Approaches 



ter of the memory follows from the analysis of the vari- 
ances of the stored quantum state. Fig. p~3] shows the ex- 
perimental variances of the state of the atomic memory 
for input light states with the mean photon number be- 
tween zero and eight. From these values the experimen- 
tal fidelity of 64% has been calculated, which is higher 
than the benchmark classical fidelity 52% for this class 
of states. 



The Faraday+feedback quantum memory and the pro- 
tocols for entanglement and spin squeezing in Sec. |lV|rely 



Duan 



on a pulse of light (or two pulses in the case of 
et al. (2000a)) interacting with the medium once. It is of 
course possible to have a pulse of light interacting several 
times with one (or more) atomic ensembles, as suggested 
in several theoretical studies. This section is intended to 
provide an overview over proposals relying on multiple 



33 



passes of light through atoms. 

Common to all these proposals is that they take ad- 
vantage of the possibility to perform phase shifts on light 
and to rotate atomic spins in between the passes. A rel- 
ative phase shift <p between the classical field and the 
quantum field, in x- and ^/-polarization respectively, will 
give rise to a rotation of field quadratures, 

X L ^coscj)X L +sm(j)P L , P L ^ cos (j)P L - sin 4>X L . (80) 

Rotation of atomic spins about the axis of polarization, 
via e.g. fast RF pulses, allows for an analogous rotation 
of Xa and Pa- Alternatively light can be sent through 
atoms from a different direction. As light is sensitive 
to the projection of the collective atomic spin along the 
axis of propagation, this is equivalent to a rotation of 
the atoms. Especially for room-temperature atoms in a 
cell, optical access from two orthogonal directions can 
be afforded trivially, while light impinging from different 
sides still talks to the same symmetric mode of atoms 
due to thermal averaging. 



The first proposal along this line is due to Kuzmich and 



Polzik (2003) and presents a protocol for atomic state 



read-out, i.e. mapping of the atomic spin state onto the 
polarization of light. As shown in Fig. [l4j in a first pass 
a pulse of light propagating along x interacts with atoms 
in a QND fashion, generating a state described by (40). 
A phase shift of <j) = tt/2 then changes X' L in = P^out 
and P' Lin = — ^L,outo where primed variables refer to 
the second pass of light. The pulse is redirected to the 
ensemble along the negative y direction, such that the 
input output relations become, taking k — 1, 



X L,out : 
P' 



zX L,in 
- P' 



- x 'a,\ 
--X 



L.in 



Paa 



•pOl,m + -PL,in) = ~ ^A,m, 



(81) 



Aside from an unimportant phase change, the state of 
atoms is mapped on light. The spurious effect of light 
noise Xl^ in the second line can be removed by using 
squeezed light. The overall input-output relations are 
similar to the ones for quantum memory for light in Eq. 
(79), but require neither measurement nor feedback. 



The previous protocol naturally raises the question 
whether in principle a perfect state transfer, or state 
swapping, could be achieved in several passes and spe- 
cific r otations on atoms a nd light without using squeezed 
light. Kraus et al. (2003) addressed this question in full 



generality and gave necessary and sufficient conditions 
for what types of quadratic Hamiltonians can be achieved 
in two modes, given a specific interaction - such as e.g. 
the Faraday interaction ~ PlPa ~ when combined with 



'local' operations of the type (80). It was found there, 
that the Faraday, or QND, interaction is most capable 
for simulating, in this sense, other interactions, while the 
beam splitter and parametric gain interaction have no 
potential to emulate any other interaction. Furthermore, 
optimal strategies for generation of squeezing and entan- 
glement were devised in the same paper. Fiurasek (2003) 



extended these results asking what types of unitary trans- 
formations (instead of Hamiltonian interactions) can be 
achieved and showed in particular that a perfect state 
swap can be performed with three passages of light, and 
not for less, see also the work by Takano et al. (2008). 



Hammerer et al. (2004) applied these general results to 



the specific situation of a light-matter quantum interface, 
showing that atomic decay and light losses can be toler- 
ated. See also Kurucz and Fleischhauer (2008) for a dis- 



cussion of memory conditions and a comparison of these 
multipass approaches to Raman and EIT based memo- 
ries. 

These protocols all assume a QND interaction in light 
and matter, which is practical for large cells with room- 
temperature atoms only when two cells and counter- 



rotating spins are used, as discussed in Sec. IV. B This 
makes multiple passes with different directions of propa- 
gation a difficult issue. Another complications is due to 
the fact that room-temperature ensembles require msec 
pulses, implying an unreasonable long delay line in the 
loop of Fig.[l4]in order to prevent the pulse meeting itself 
in the atomic medium. These problems were overcome by 



Sher son et ^ ( |2006a| ), |Fiurasek eta l. (2006) and |Muschik 



et al. (2006 



who showed that protocols for quantum 
memory and entanglement generation involving multiple 
passes of one or more pulses can be matched to Larmor 
precessing ensembles and that light traversing atoms si- 
multaneously from different directions can in fact be ad- 



vantageous. In particular, iMuschik et al. (2006) consider 



a single cell Larmor precessing in a magnetic field ori- 
ented along x in a setup as shown in Fig. [l4j assuming a 
loop length much smaller than the pulse length, such that 
the pulse "meets" itself in the medium. For atoms rotat- 
ing in the sense of the light propagating along the loop, 
solution of the corresponding Maxwell-Bloch equations, 
taking care of propagation effects, yields input-output 
relations for atomic variables, 



Xi.out = e-* 2 ' 2 X AM + Vl - e-* 2 X+ in , 
-PA.out = e 



- K2 /2p 4 . + ^^ ? P+ 



L,in> 

f- 
L,in> 



where X~£ in , P^ in refer to a light mode centered at the 
upper side band frequency Uq + uol (^o is the carrier fre- 
quency of the classical driving pulse and ujl the Larmor 
frequency) with a weakly exponentially decaying mode 
function, n is again given by (41). An analogous input- 



output relation holds for light, such that this scheme re- 
alizes an exponentially efficient state exchange of atoms 
and light. 

The Faraday interaction was achieved as a sum of the 
beam splitter and the parametric gain interaction, c.f., 
Fig. [2] What we effectively achieve by having multi- 
ple passes is that we add two Faraday interactions with 
different relative phases between the beam splitter and 
parametric gain interactions in Fig. [2] (c). A suitable 
choice of geometry and phase shifts leads to cancelation 
of the parametric gain interactions after the two passes, 
and we are left with the beam splitter interaction making 
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an ideal memory transformation. With other configura- 
tions it is the beam splitter interaction which cancels 
and the parametric gain interaction takes effect, gen- 
erating an entangled state of light and atoms, whose 
EPR variance, cf. Sec. |IV.B[ sca l es asy mptotically as 
Aepr ~ exp(-^ 2 ) flMuschik e^"Z||2QQ6|. 



t Qo a 




FIG. 14 (a) Setup for two pass protocol for read-out of atomic 
states as suggested by Kuzmich and Polzik (2003) and, with 
an additional magnetic field applied along x, for full state ex- 
change or creation of entanglement between light and atoms 
as suggested by Muschik et al. ( 2006 ) . (b) State swap between 
two ensembles coupled to a common ca vity mode based on 
adiabatic passage as demonstrated by Simon et al. (2007b). 



The classical pulses Qi are applied in a counterintuitive se- 
quence. 



To see what happens in this situation let us now con- 



sider the A -> limit of the solution in Eqs. (46), (47), 
and ( 50 ) . If we assume a sufficiently large incident clas- 



sical driving field so that h(0,i)z/L > 1 we may use 
the asymptotic form of the Bessel function and write the 
integral kernel (47) as 



m(Q;t,z) 




L2y/2ittfvgtz 



d{yf^i-^) 2 /L 



(82) 

where we have for simplicity assumed that the classical 
driving field is time independent and have introduced the 
group velocity 



n 2 L 



(83) 



For a large optical depth this kernel is centered around 
z = v g t. If the spin wave mode is slowly varying on 
a length scale L/y/d, or if the incoming field is slowly 
varying on a time scale L/\fdv g , i.e., the input field is 



inside the EIT-transparency window ( Fleischhauer et al. 



2005J), the expression above can be approximated by a 

Vgt). The solutions of Eqs. 



delta function m = ^/Vg~5(z 
(46) and (50) then become 



D. Raman and EIT approach 

The beam splitter interaction is a quite natural choice 
for a quantum memory since the interaction maps exci- 
tations from light to atoms and back. The beam split- 
ter interaction was pro posed for a quantum memory by 
Kozhekin et al. (2000) who considered the far-detuned 
Raman limit A ^> dj. The same Raman limit was also 
considered by Nunn et al. (2007). Most of the atten- 
tion to the beam splitter interaction, however, came with 
the realization that on resonance with the atomic tran- 
sition, Electromagnetically Induced Transparency (EIT) 
can be used for quantum memory for light. The prin- 
ciples of EIT and its applications for coherent memory 



for light have been reviewed earlier ([Fleischhauer et al. 



2005 Lukin , 2003 ) , hence we will concentrate here mostly 



on the recent advances on quantum memory. The EIT 
is achieved when a strong control optical pulse renders 
the A system transparent for the signal pulse, see Fig. 
[ijb). This transparency is accompanied by a strong re- 
duction in the group velocity of the signal pulse. The re- 
sult is that the signal pulse entering the atomic ensemble 
is spatially compressed. If the compression is sufficient 
to make the signal pulse fit inside the sample, the control 
field can be turned off at this point and the signal pulse 
is "frozen" into the atomic ground state coherence - the 
dark state polariton wave. The process can be inverted 
by turning on the control field after some delay, which 
leads to the generation of a signal pulse which "remem- 
bers" the classical and quantum properties of the input 
signal pulse. 



«A,out(^) 
0>L,out(t) 



za L ; m (T - Z/Vg) 



/Vg~a A ,m(L-v g t). 



(84) 



The first line describes the writing of the input field into 
the atomic memory, whereas the second one describes 
the readout of the atomic memory back into light. To 
accomplish the writing (storage) process the control field 
must be turned off when the entire input pulse is in- 
side the medium, which requires T s i gna i < L/v g . This 
prompts the reduction of v g to zero and hence " stopping" 
of light. The process of the readout or the retrieval is ac- 
complished according to the second line by turning the 
control field on again which leads to the mapping of the 
atomic memory operator back onto the light operator. 

Recently several experimental implementations of the 
EIT atomic memory ( Chaneliere et al. t 2005| Choi et al. 



2008 Eisaman et al. 2005 Kuzmich et al. 



2003) have 



demonstrated that such quantum features of light as vio- 
lation of a Cauchy-Schwarz inequality and entanglement 
can be preserved by the memory. Very recently the 
EIT based memory for quantum fluctuations has been 
also demonstrated using squeezed vacuum light in (|Ap-| 



pel et al. 2008 Honda et al.\ 12008b. As mentioned in 
(Honda et al. 2008) a fidelity calculation along the lines 
of the discussion in Sec. |V.A| would be misleading for the 
case studied in these two papers since only one particular 
state has been used and, in addition, the overlap of this 
weakly squeezed state with vacuum is not far from unity. 
An overall efficiency of storage and retrieval of around 
10 — 15% has been achieved, and weakly squeezed light 
(around — 0.2dB) has been retrieved from the memory. 
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The experiments ( Chaneliere et al. t 2005 Choi et al. 



2008 Eisaman et al. p 2005 IKuzmich et al, 2003) have 



had a relatively low overall efficiency of the storage- 
retrieval process, however it did not preclude observation 
of the storage of non-classical light. The reason is that 
the violation of the Cauchy-Schwarz inequality is based 
on the measurement of the normally ordered second order 
correlation function g 2 (l,2) = (: fi\fi2 '•) / {fh\) {fb2) where 
ni 5 2 denote photon number operators measured at space- 
time points 1,2. g 2 (l,2) which describes the normalized 
probability of detecting a photon at point 2 conditioned 
on the detection at point 1 is insensitive to losses be- 
cause they affect the numerator and the denominator in 
the same proportion. For an ideal single photon state 
g 2 (l,2) = 0, whereas any g 2 (l,2) < 1 is a signature of 
the non-classical character of the field which for the case 
of a stationary photon flux is referred to as photon anti- 
bunching. 



Here et al\ 2005 ; Eisaman et al. , 2005 ; Matsukevich et al. 



The source of the non-classical fie ld used in flChane- 



2006||Yuan et al.||2007D, and flChoi et al.||2008) has been 



an atomic ensemble prepared in an approximate atomic 
single excitation state by a weak parametric gain interac- 
tion (see Sec. IV. C) which was then retrieved onto light 
by a beam splitter process. Non-classical states pr epare d 
in this way have been re ported by Balic et al. (2005) 
IChaneliere et al.\ fl2005[ ); JChou et al.\ Q2004D; |Du et al 



Matsukevich and Kuzmich (2004); van der Wal et 
( 20031) . The sequence is similar to that described in Sec 



( 200 8P; |Eisaman et al\ Q2005P; |Kuzmich et a/.|"|2 003 ) 



IV. C except for only one atomic ensemble is involved 



at this stage. The lay out of the experimental setup of 
Eisaman et al. (12005) is shown in Fig. 15 The weak 



parametric gain interaction driven by the "write" clas- 
sical field so that k«1 creates a single atomic excita- 
tion in the "source" ensemble conditioned on the detec- 
tion of a photon in a particular spatial "Stokes" mode. 
The requirement of low gain n <C 1 is necessary so that 
multi-photon pulses are suppressed, which is a typical 
situation for the parametric- type interaction. The "re- 
trieve" strong pulse converts the atomic excitation into 
an "anti-Stokes" light pulse. The non-classical charac- 
ter of this field is manifested by the correlation func- 
tion conditional on the detection of one Stokes photon 
g 2 (AS || ns = 1) < 1. This condition means that if a 
Stokes photon has been detected and if there is a photon 
in the anti-Stokes pulse, then the probability of having a 
second photon in the anti-Stokes pulse is less than that 
for a random process. In the ideal case this probability is 
zero and the anti-Stokes pulse contains either no photons 
or just a single one. 

The non-classical light pulses produced by the " source" 
ensemble are then directed towards the memory ensem- 
ble. There they are stored and retrieved using the EIT 
pulse sequence as shown in Fig. [16] from |Choi et 
(2008). The figure shows the probability of detecting 



the quantum input pulse appears at r = 0. Then the 
driving field is turned off leading to the storage of the 
light in the medium. For ideal storage there should be 
no counts corresponding to the input pulse. The counts 
around r = hence correspond to the "leakage" of light 
through the memory ensemble. After a delay time of 
1 /isec the strong "read" control field is applied and the 
atomic excitation is read out generating a light pulse as 
prescribed by Eq. ( 50 ) . An overa ll storage-retrieva l prob- 
ability of 17% has been sho wn in|Choi et a/.| ( |2008[) . This 
experiment and the one by Chaneliere et a/.| ( |2005 ) used, 
respectively, Cs and Rb atoms cooled and trapped in 
a mag neto-optical trap (MOT), whereas Eisaman et al. 
(2005) used Rb atoms in a neon buffer gas cell at room 



temperature. In the two former experiments a small sub- 
ensemble of atoms contained within the volume of the fo- 
cused optical beam inside the sample served as the mem- 
ory. The storage time in these experiment was restricted 
by the motion of atoms and/or by magnetic field inho- 
mogeneity to about 1 /isec. The two lower levels of the 
A system were the two ground state Si/ 2 hyperflne levels 
and the upper level was the P 3 / 2 level. 

The results for the conditional correlation function 
g 2 (out || ns = 1) < 1 for the me mory output field ob- 
tained by Chaneliere et al. ( 2005[ ) are shown in Fig. 17 
The results clearly show the nonclassical character of 
the retrieved light. The overall efficiency of the stor- 
age/retrieval process of a photon was 6.4 %. Similar re- 



sults have been obtained by Eisaman et al. (2005). 



The experiment by Choi et al. ( 2008 ) have taken the 



EIT approach one step further by demonstrating it for 
an entangled state, more precisely for a superposition 
state of a photon being in one of two possible pathways 
(Fig. [l8|). In this experiment a conditional non-classical 



a photon after the storage medium. The strong "stor- 
age" driving EIT field has a constant amplitude until 



state with g (AS \\ ns = 1) < 1, an approximate sin- 
gle photon state, was split on a beam splitter in two 
parts. Conditioned on the registration of the Stokes pho- 
ton in the source ensemble (not shown in Fig. 18) the 
input consisted of a single photon with 15% probability. 
The rest was mostly vacuum with a small addition of 
a two-photon component with the probability of 9% of 
that for a Poisson source with the same average photon 
number. The input state with these properties has the 
concurrence of 0.10 corresponding to an entanglement of 
formation of E Q p = 0.025 ebits. (A concurrence of unity 
corresponds to E q f = 1 ebits, i.e. to a maximally entan- 
gled Bell state.) The input state was linearly polarized at 
45° with respect to the polarizing beam splitter (Fig.[l8|). 
The photonic state of the two outputs of the beam split- 
ter conditioned on the successful preparation of the single 
photon state is * in = (\0 L )\1 R ) + e^|l L )|0 H ))/V2. The 
two components of this state have been directed into two 
atomic ensembles. The two ensembles were two groups of 
atoms within a Magneto-optical trap (MOT). Rb atoms 
were optically pumped into a particular magnetic sub- 
level of the ground state F = 4, nip = which con- 
tributed to a better performance of the memory. 

After the EIT storage and retrieval steps similar to 
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FIG. 15 EIT-based memory setup after |Eisaman et al. ( 2005 ). a) Level scheme and sequence of control pulses. Left - parametric 
interaction for a probabilistic generation of an excitation in the source ensemble, center and right - EIT-based beam splitter 
interactions for storage and retrieval, b) Setup showing the source ensemble of the non-classical photon flux and the memory 
ensemble (target atoms) . A count in detector SI or S2 is a condition for starting the process. Detectors AS1 and AS2 analyze 
the statistics of the input and output memory photons. 



those described earlier in this section the two retrieved 
components were combined on a polarizing beam splitter. 
With a suitable choice of the phase the state recreated 
from an ideal memory would make the single linearly po- 
larized photon, just like the input one. Changing the 
phase between the two components of the state stored 
in the two ensem bles by adjusting the A/2 wave plate 
Choi et al. (2008) performed tomography of the entan- 
gled state and obtained the density matrix of it. The 
results of this procedure are shown in Fig. 19 From 



by |Gorshkov et al.\ J2QQ7a|b|c|cH |2008|> (see also related 
work by Dantan et al. (2005, 2006b|);|Dantan and Pinard 



these results the concurrence of the retrieved state of 
0.017 (E q y = 0.001 ebits) has been inferred demonstrat- 
ing that the retrieved state has retained entanglement 
after the storage process. 

The overall efficiency of the storage/retrieval of entan- 
glement measured by the ratio of the concurrence of the 
output to the concurrence of the input is 20 % after the 
storage time of 1.1 /xsec. It is limited by the finite optical 
depth of the sample and can be improved by optimization 
of the pulse shapes as discussed below. 

These experiments represent exciting experimental 
progress, however a higher efficiency is still desirable for 
future applications. Higher efficiency can be achieved 
with an increased optical depth, but also by optimizing 
the shape of the classical driving field Q(t). This op- 
timization problem is the subject of a series of papers 



d2004| ) and by |Nunn et al.\ p007| ). These studies show 
that for any slowly varying pulse of duration T such that 
Tdj ^> 1 (bandwidth BW <C dhf) the optimal storage 
and retrieval efficiency is the same and is independent 
of the detuning from the excited atomic state. Further- 
more, the optimized inefficiency only depends on the op- 
tical depth d and scales as 1/d. These results can be 
understood from a new "universal" physical picture of 
the storage and retrieval process: First of all for a given 
stored spin wave, the retrieval process is essentially a 
constructive interference effect similar to super radiance, 
where the radiated fields from all atoms interfere con- 
structively in a certain direction. As a result there is a 
fixed branching ratio between the decay into the desired 
quantum field mode and the decay into all other modes, 
which is independent of the control field shape as long 
as sufficient optical power is used. Secondly, the opti- 
mal storage procedure is the time reversal of the retrieval 
procedure, and, by time reversal symmetry, the optimal 
storage efficiency is identical to the retrieval efficiency. 

The retrieval efficiency does, however, depend on the 
spin wave mode, and for optimal storage the classical 
drive field Cl(t) has to be chosen so that it maps an in- 
coming field mode into the optimal spin wave. This opti- 
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FIG. 16 EIT-based memory flChoi et a/.|fl2008| )). Probability 
of detecting a photon downstream the memory ensemble (left 
axis). The points around r — 0, the moment when the input 
pulse is launched are due to the imperfect memory leading to 
the "leakage" of photons through it. The strong field (right 
axis) is turned off around r = and turned on again at r = 
1/xsec leading to the retrieved pulses with the overall storage 
retrieval probability of 17% 




0.03 



FIG. 17 The correlation function for the inp ut (b) and output 
(d) lig ht for the EIT-based memory from |Chaneliere et al. 
(2005). g 2 (0) < 1 is a signature of non-classical character 



corresponding to the single photon in the limit of g (0) = 0. 



mal field shape £l(t) can be found by a direct calculation 
but an alternative experimental procedure for finding op - 
timal shapes is demonstrated by |Novikova et al (2007a). 
In this experiment which worked in the classical regime 
with many photons in the signal beam, the shape of the 
pulse to be stored was optimize d for a given classical drive 
field. Novikova et al. ( |2007a ) first stored a given pulse 
and recorded the shape of the retrieved pulse. Then the 
time reverse of the recorded pulse was used as an input 
for the next round of storage, retrieval and measurement. 
This procedure rapidly converges and yields the optimal 



efficiency (Gorshkov et al. 2007c), which in the experi- 



ment was in the range 42-45%. In Novikova et al. ( 2008 ) 



on the other hand the full theoretical optimization of 
the classical field shape Q(t) was used to store arbitrary 
field shapes and retrieve them into a possibly differently 
shaped mode with a similar efficiency. 

Because the optimal strategy for storage and retrieval 



is based on time reversal, higher efficiency can actually 
be achieved if the excitation is read out in the backward 



direction compared to the direction of storage ( Gorshkov 
\et a/.|[2007c| ). This change of direction, however, requires 
a redefinition of the atomic operators (29). The mode 
functions u m (z;f_\_) in Eq. (29) are solutions to the 3- 



dimensional Maxwell equation in the forward direction 
but it may not be so in the backward direction. This will 
complicate the dynamics unless the mode function can 
be chosen real, which requires a Fresnel number much 
bigger than unity (for a discussion of a related problem 
see Andre (2005)). Furthermore a finite energy difference 
ojqi between the two ground states introduces a momen- 
tum difference Afc = ujoi/c which red uces the achiev- 
able memory efficiency unless AkL < 1 ( [Gorshkov et al.\ 
2007c). A way to cope with this problem is presented by 



Surmacz et al. (2008). 



For applications in a quantum network, stored excita- 
tions will have to be processed in the quantum memory. 
A first step in this direction was performed in recent 

~ ( |2007cD and |Simon et al. 



experiments by Simon et al. 



(2007b). The latter experiment involved two atomic en- 
sembles in a medium Finesse (J 7 = 240) cavity. It demon- 
strated adiabatic transfer of a single excitation stored 
in one ensemble to the other ensemble with the cavity 
mode serving as a quantum bus, as well entanglement 
of the two ensembles by a partial transfer. First, a sin- 
gle excitation is generated in one ensemble by driving a 
weak parametric gain interaction with subsequent detec- 
tion of a Stokes photon, as described above. Conditioned 
on the successful generation of a single excitation in en- 
semble 1, beam splitter interactions are switched on and 
couple both ensembles to a common cavity mode, as in- 
dicated in Fig, p^b). This is done adiabatically, first for 
the "empty" ensemble 2 and then for ensemble 1, in a 
counterintuitive sequence generating an adiabatic dark- 
state passage 1 1) 1 10) 2 — >• |0) 1 1 1)2 - After the transfer, the 
single excitation was read out from ensemble 2, demon- 
strating a transfer efficiency between 10% and 25%, de- 
pending on the optical depth. Simon et al. (2007b) also 
demonstrated a partial swap of the excitation, generat- 
ing in the ideal case an entangled state of the ensembles, 
Jl)i|0) 2 -> cos0|l)i|O) 2 +exp(^)sin<9|0)i|l) 2 , where 
is controlled via the intensities and <\> by the relative 
phase of the laser fields Qi in Fig. [l4|b). Reading out 
the collective excitations and measuring photon correla- 
tion functions, similar to what was done by |Chou et al\ 
(2005), a lower bound on the entanglement of the ensem- 



bles was determined, giving a concurrence larger than 
0.0046 corresponding to an entanglement of formation of 
E oi > 0.0001 ebits ( |Wootters[|1998l ). 



E. Photon Echo 



It has long been known that the photon echo technique 



(Kurnit et al. 1964) could be used to store classical light 



pulses. Recently it has been realized that one can ex- 
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FIG. 18 EIT-based memory for an entangled state (Choi et al. 2008). Two memory ensembles store two components of an 
entangled state which are later retrieved and analyzed by two pairs of coincidence detectors D\^- 
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FIG. 19 Characterization of the retrieved entangled state 
(Choi et al. 2008). The reconstructed density matrices of 



the input (a) and the output (b) states of the light. 



tend these techniques into the quantum interface domain 
flKraus et al\ |2006| |Moiseev and Kroll[ |2001[ ). If a light 
pulse is absorbed by an ensemble of two-level atoms, the 
quantum properties of light will be stored in atoms as 
shown already in the early experiment by |Hald et~ 
(1999). The problem is that for a meaningful memory the 



coherence time of the optical transition should be longer 
than the duration of the pulse. However, this means that 
the bandwidth of the interaction which is set by the in- 
verse coherence time is narrower than the bandwidth of 
the light set by its inverse duration, so that the entire 
pulse cannot be stored. A way out of this problem is 
to use a medium with inhomogeneous broadening. Then 
different frequency components of the light pulse will be 
effectively stored in different sub-groups of atoms. To 
avoid the dephasing of the stored state caused by the 
broadening a photon echo technique is used which also 
allows to control the release of the stored excitation. 

The essence of the photon echo approach is to have an 
inhomogeneously broadened line which is then reversed. 
In the original approach an incoming light field is ab- 
sorbed by a two level system with a broadened optical 
transition. Due to the inhomogeneous broadening the 
optical coherence from each atom (i) precesses at differ- 
ent frequencies exp(— ioj{t). These different precession 
frequencies dephase the optical polarization such that 
it does not radiate because the radiation from different 



atoms interfere destructively. In the simplest version of 
the photon echo, a strong 7r-pulse, which interchanges 
the ground and excited states, is applied after a time 
T/2. This strong 7r-pulse effectively reverses the phase 
acquired by each atom such that the subsequent time 
evolution causes a rephasing of the optical coherence. At 
time T all the atomic polarizations are again in phase 
causing an echo signal to be emitted. 

A modification of the photon echo technique 
called Controlled Reversible Inhomogeneous Broadening 
(CRIB) which in principle allows an ideal memory effi- 
ciency was introduced by Nilsson and Kroll (2005) and 



Moiseev and Kroll (2001 ) and further developed by Kraus 



et al.|Q2006D . |Moiseev arid Kroll| ( |200l[ ) considered a A- 
system similar to Fig. [lib). The photon is absorbed 
while an applied electric field broadens the |0)-|e) ab- 
sorption line. A ir pulse traveling in the same (for- 
ward) direction is applied, which takes the population 
from the excited state to the initially empty state |1) for 
long term storage. Later on another 7r-pulse traveling in 
the opposite (backward) direction is applied, which re- 
leases the excitation. During the retrieval process the 
broadening with the applied electric field is reversed, 
which reverses the time evolution resulting in an ideal 
retrieved signal in the backward direction, provided that 
the excited state decoherence is negligible. Experimen - 
tally this scheme was realized by Alexander et al. ( 2006 ) . 



Staudt et al. (2007) demonstrates that similar photon 



echo techniques can preserve the coherence of time bin 
qubits conditioned on having an outgoing photon. There 
are by now many versions of this CRIB approach, dif- 
fering, e.g., by whether the broadening arises from dif- 
ferences in the response of the atoms to a homogeneous 
field (transverse broadening) or from spatially dependent 
broadening (longitudinal broadening). A full account of 
all differe nt types is beyon d the scope of this review (see 
e.g. |Longdell et al. (2008) for a theoretical discussion of 
several different schemes). 



In a recent experiment (Hetet et al. 2008a) CRIB has 



been implemented in a particularly elegant and simple 
way. A simil ar technique was also considered theoreti- 
cally by Sangouard et al. (2007a). In contrast to most 
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other techniques considered in this paper, no strong clas- 



sical control pulse is required in this scheme. |Hetet etal\ 
(2008a) applied a gradient electric field to the crystal 
used for the memory. Light interacted with Pr 3+ ions 
which were doped into a Y 2 Si05 host. Due to the gradi- 
ent field the resonance frequencies of the ions at different 
parts of the crystal were distributed by the Stark effect 
within the 2 MHz bandwidth. The homogeneous line 
width of the ions is 100 kHz. Hence the ions should have 
been in principle capable of storing a pulse of light for 
this time which is 20 times longer than the pulse dura- 
tion. The stored light has been released by reversing the 
sign of the electric field gradient. In the experiment a 
classical pulse of 2 /isec duration was stored for 2 /xsec. 
An overall efficiency of the storage and retrieval of 15 % 
has been achieved. The results for the echo memory for 
classical pulses are shown in Fig. [20] The figure shows 
the transmitted part of the input pulses (centered around 
zero) and stored and retrieved pulses for different storage 
times. All pulses are normalized to the amplitude of the 
input pulse. 
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FIG. 20 Transmitted and retrieved pulses for different storage 
times for an echo based memory (private communication from 
Hetet et al. ( 2008a[ )). The small peak at t ~ represents a 
small leakage (non stored component) of the incoming field, 
whereas the later peaks are the retrieved field. The emission 
time of these light pulses is controlled by reversing the sign 
of an applied electric field halfway during the storage period. 
See comments in the text. 

An important question is to which extent the addi- 
tion and reversal of broadening improves the memory 
performance compared to the other approaches consid- 
ered h^re^^or_th^approaches of Hetet et al. (2008a) 
and Sangouard et al. ( |2007a ) a major practical advan- 
tage is that it does not rely on any classical laser fields. 
This, however, limits the attainable storage time to the 
coherence time of the optical transition, which is typ- 
ically shorter than for the ground states. The advan- 
tage of reversible broadening for the storage of a single 
mode in the ground state cohe rence was investigated by 
Gorshkov et al. (2007d 2008) where it was found that 



little was gained. The situation is, however, different 
for the st orage of mult i ple mod es. I n particular, it was 
fl2007bp 



found by Simon et 



and 



Nunn et 



(2008) 



modes 



that CRIB changes the scaling of the number of 
which can be stored from \fd to <i, thus significantly in- 
creasing the multi mode capacity. Furthermore a novel 
type photon echo approach using atomic frequency combs 
was recently proposed ( [Afzelius et al\ |2009| ) . This pro- 
posal takes advantage of a large inhomogeneous broad- 
ening in a solid state system. By exploiting atoms with 
different resonance frequencies one can increase the effec- 
tive number of atoms participating in the memory and 
thereby achieve a very efficient quantum memory with 



high mult imo de cap acity. Using this approach |de Ried- 
|matten et al.\ ( |2008[ ) demonstrated the coherent mapping 
of light at the single photon level onto a neodymium ion 
doped crystal, and collective release of the stored light 
at a pre-determined time. Moreover, this experiment 
demonstrated the storage of pulses in different temporal 
modes, proving the multimode capacity of this approach. 



VI. QUANTUM TELEPORTATION BETWEEN LIGHT 
AND ATOMS 

A. Quantum Teleportation 

Quantum teleportation is a means for sending quan- 
tum states from A to B in a disembodied fashion using 
two separate channels, a quantum channel connecting A 
and B and a classical one. It makes use of entanglement 
shared via the quantum channel and classical commu- 
nication. Apart from being one of the most surprising 
and mind-boggling discoveries in quantum information 
theory, quantum teleportation has become an essential 
primitive in quantum computation and quantum com- 
munication. 

Shortly after the first theoretical layout for quan- 



tum teleportation of states of a qubit ( [Bennett et al. 



1993[ ) the pro tocol was extended to s t ates of continu- 

Vaidmanl 



1998 



ous v ariables (Braunstein and Kimble 

1994). Both sorts of protocols were first demonstrated 

with light, utilizing either probabilistically generated Bell 



states of photon pairs ([Bouwmeester et al.\ 1997 ) or de- 



terministically generated EPR beams (Furusawa et al. 



119981), both obtained from parametric down conversion. 
The first teleportation involving massive particles was 
performed recently in th e ion trap experim ents at Inns- 
bruck (iRiebe et a/.l|2004|) and NIST {Barrett et a/.l|2004). 



For a recent review on teleportation of states of continu- 



ous variables see (Furusawa and Takei 2007). 



Our focus here is on teleportation protocols involving 
both matter and light. In such a scenario, quantum states 
carried by traveling pulses of light are teleported onto a 
stationary quantum memory, for example, an ensemble 
of neutral atoms, using the quantum interface. The first 
realization of teleportation involving matter and light uti- 
lized the entangled state created via the quantum Fara- 
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day - QND interaction of a pulse of light with the collec- 



tive spin of an atomic ensemble (Sherson et al. 



2006c). 



In this experiment deterministic teleportation with the fi- 
delity higher than any classical state transfer can achieve 
has been demonstrated. Very recently probabilistic tele- 
portation between light and matter has been also demon- 
strated using parametric gain- and A-type interactions 
( [Chen et al.\ |2008[ ) and has been extended to entangle- 
ment swapping, that is the teleportation of entanglement 



(Yuan et al. 2008). 



In a nutshell, a teleportation protocol involves three 
steps: First, an entangled state is created and shared as 
a resource between two stations (usually termed "Alice" 
and "Bob"). In the cases considered below, the entan- 
glement between atoms (Bob) and light is generated by 
sending a strong driving pulse through atoms. As a re- 
sult of this interaction the forward scattered photons of 
the orthogonal polarization sent to Alice become entan- 
gled with atoms kept by Bob as shown in Fig. [2l] The 
ideal continuous variable EPR entanglement corresponds 
to A(X A + X L ) 2 = A(P A -P L ) 2 -)► 0. Alice also receives 
another, unknown quantum state described by canonical 
variables Y, Q sent by a hypothetical sender " Victor" , 
which is to be teleported to Bob. For this, Alice per- 
forms a joint measurement of Xl + Y and Pl~Q, called 
a Bell measurement, on the photonic part of the entan- 
gled state that she has received and the unknown quan- 
tum state of light to be teleported. This measurement is 
performed by mixing two light pulses on a beamsplitter 
as in Fig. [21] The results of this measurement are com- 
municated to Bob. Bob uses this classical information to 
perform a correcting operation on his quantum system - 
atoms - by shifting Xa,uu = Xa + (Xl + Y) —> Y and 
Pam = Pa — (Pl — Q) —> Q, thereby recovering the 
original unknown state. 

We see that in the hypothetical case of vanishing EPR 
variances, first and second moments and thus any Gaus- 
sian state (i.e. a state with Gaussian wave function or 
Wigner function) is transmitted perfectly. This in turn 
implies that any (non-Gaussian) state, including a qubit 
of the form a\0) + /3|1) G L 2 (R), would be teleported 
faithfully, as the set of coherent states is a subset of all 
Gaussian states and provides a basis for the full Hilbert 
space L 2 (R) (for the proof of this in the Schrodinger pic- 
ture and for the Wigner function, see |Braunstein and van| 
|Loock| ( [2005[ )). In this sense the distinction between tele- 
portation protocols for qubits and continuous variables 
is superficial. This said, for realistic cases of imperfect 
teleportation where fidelity is not perfect, performance of 
each teleportation protocol should be evaluated in detail 
having in mind a particular application. 

If EPR variances do not vanish, as is necessarily the 
case due to energy restrictions, teleportation will not be 
perfect, and it is necessary to evaluate the performance 
of the teleportation. The teleportation is essentially a 
protocol mapping the state of one system to another. 
We can therefore use the fidelity as the figure of merit 
as discussed in Sec. V.A[ i.e., find out when the perfor- 



mance of the protocol becomes better than that of the 
best classical protocol for a given class of input states. 



B. Teleportation based on Faraday interaction in magnetic 
field 



In Sherson et al. (2006c) the entangled state of light 
and atoms for teleportation is obtained via the Faraday 
interaction of light with a single collective atomic spin, 
precessing in an external magnetic field. The relevant 
level scheme is shown in Fig. [4] and in the inset to Fig. 
|2lj The situation is described by the Hamiltonian and 



Maxwell-Bloch equations given in section II. E In this 
case, the Hamiltonian does not fulfil the QND criteria 



(Holland et al. 1990 Poizat et al. 1994) because the 



Faraday interaction Hp, c.f. Eq. (34), does not com- 



mute with the free Hamiltonian (63) describing Zeeman 



splitting of ground states. In fact, as seen from the inset 
in Fig. 21, the interaction resembles the Raman pro- 



cess. The Maxwell-B loch equation s (65) were integrated 
by Hammerer et al. (2005a |2006|) and the solutions are 



again conveniently expressed in terms of experimentally 
measurable cosine and sine modulation modes (66), 



XA,out — Xa,iu H ~7=PL c ,im 

Pa, out = Pa, in H j=Pl s , \n-> 

-* L c ,out = -* L c ,in? 
-^L c ,out — -^L c ,in H y^ Pa, in 



1 fK\ 2 



(2) ^-^u) ^ 



-*L s ,out — -* L s ,in? 
-^L s ,out — -^L s ,in T^Xa 



V2 

( K\ 2 p 1 ( n \ 2 JD 

- {2 ) PL - in ~^UJ PL °*«*>' m ' 

The terms proportional to n 2 are a new feature, specific 
for this setup, and represent atom-mediated back-action 
of light onto itself. This effect involves the previously de- 
fined cosine and sine modes as well as yet another pair of 
canonical independent "back-action modes" Pl c{s) bacfc ,in 
which can be treated as vacuum noise operators. To see 
that this interaction can be used for teleportation let us 
inspect the EPR type correlations between the atomic 
mode XA,out, Pa, out and the light mode of the upper 
sideband X + , P + with the frequency u + ul 



XL + ,out — —7=\Xl s , out — Pl c , out) 5 
-PL + ,out = ~^=(^L c ,out +-PL s ,out)- 
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FIG. 21 Experimental setup for the teleportation of light to atoms (Sherson et al. |2006c| . Strong y-polarized pulse interacts 
with atoms and generates an entangled x-polarized mode at the upper sideband frequency (inset). This entangled light is 
overlapped with the pulse to be teleported on a 50/50 beam splitter BS after which the Bell measurements are performed. The 
results of the Bell measurements are fed back onto atoms via RF magnetic coils. Detailed comments in the text. 



One easily finds an EPR variance of 

AePR = A(XA,out+X L 



out) +A(Pa 5 ou t - 



PL + ,out) 



where the lower bound is achieved for n ~ 1.48. Hence 
indeed the modes of light and atoms are in an entangled 
state which can be used for teleportation. 

The state to be teleported is encoded in the lower 
sideband uj — ujl with respect to the carrier frequency 
(Fig. 21), expressed in terms of measurable cosine and 
sine modulation modes as 



Y 



1 

71 



(y B + Q c ), q 



i 

71 



(Y C -Q B ), (86) 



where [Y, Q] = i. The Bell measurement of the commut- 
ing observables after combining the entangled and the 
to-be-teleported states yields 



X c 



1 
1 

71 



(^L C) out +Y C ) , X s 
(Pl c , out — Q c) •> Qs 



1 

71 

1 

71 



(^L s ,out + ^s) J 

(87) 

{Pl s , out — Qs) • 



Conditioned on these results the atomic state is then dis- 
placed in order to get in the ensemble average the final 



state 
-Pa, fin = Pa, out ■ 



\- X s — Q c — (Xa, out + Xl + ^ out ) + Y 
- X c — Q s = (Pa, out — Pl + , out) + Q- 



In the hypothetical case of vanishing EPR variances of 
PGi.out +^L_,out) and (Pa, out -^l_, out) atoms would 
correctly display the statistics of the Y, Q mode, repro- 
ducing any input state (coherent, Fock, etc) as desired. 
For the given minimal EPR variances of 0.66, that is for 
a variance of 0.33 in each EPR variable Xa, out +-^l + , out 
and Pa, out — Pl + , out, teleportation will not be perfect, 
but still below the classical limit corresponding to the 
total EPR variance of 2. This discussion ignores atomic 
decay and light absorption which can be included as de- 
scribed in section Hi. CI 

The expe rimental implementati on of this teleportation 
protocol by Sherson et al. ( 2006c ) was performed with a 



new generation of paraffin coated cells filled with Cae- 
sium similar to those shown in Fig. [8] Atoms are ini- 
tially prepared in a coherent spin state by a 4 msec 
circularly polarized optical pumping pulse propagating 
along the direction of the magnetic field, into the sublevel 
F = 4, rriF = 4 of the ground state (insert in Fig. 21). 
Then an entangled light-atoms state is generated by send- 
ing a strong pulse polarized along the y axis (Fig. 21). 
The initially vacuum state of the x polarization of this 
pulse is populated after the interaction with the field 
Xl, out 5 Pl, out which is entangled with the atomic vari- 
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ables X a, out? Pa, out according to ( |85| ). (Compared to the 
theoretical derivation of the Faraday interaction in Sec. 
|TI|the polarizations of the classical and quantum fields are 
interchanged. This is of minor importance since it merely 
swaps the vertical and diagonal transitions in Fig. [2j) 

The nearly optimal value of the coupling constant 
k, ~ 1 was achieved with 4 x 10 13 photons in the strong 
pulse with the duration of 1 msec and a crossection of 
4.4 cm 2 detuned by 825 MHz and a number of atoms 
on the order of 10 12 corresponding to the Cs tempera- 
ture of 25° C. At Alice's location the mode of light en- 
tangled with atoms is combined with the input pulse to 
be teleported on a 50/50 beamsplitter BS (Fig. [21). The 
strong y polarized pulse which travels along with the en- 
tangled quantum field conveniently serves as the local os- 
cillator for the polarization homodyne measurements of 
the Stokes operators S y and S z performed at two outputs 
of BS. The output cos(u Lt) and sm(u Lt) components of 
the photocurrent are processed by the lock-in amplifiers 
to produce the feedback signals Q c ,s,Yc,s- The two feed- 
back signals at ujl = 322 kHz phase shifted by tt/2 with 
respect to each other are fed into the RF magnetic coils 
surrounding the atoms with a variable electronic gain. 
The gain is chosen so that the atomic variables are shifted 
by one vacuum unit if the input light mode contains one 
vacuum unit of excitation. This condition corresponds to 
the "unity gain" teleportation. 

To prove the success of the teleportation protocol a 
strong verifying pulse reads out the atomic operators 
(collective spin projections). The same homodyne polar- 
ization measurement setup is used for that. An example 
of the atomic state readout is shown in Fig. 22 The 



gain of 0.95 is found from the slope of the linear fit to 
the data, whereas the variances of the final atomic state 
after teleportation AX^ )OUt = AP^out = 1-2 are found 
from the variance of the distribution of the data. These 
variances are below the classical teleportation limit cor- 
responding to three units of vacuum noise, that is to 3/2. 
Accordingly the fidelity of the teleportation calculated as 



F=(n(l-£) 2 + l/2 + AX 2 t ) 



-1/2 



(n(l-£) 2 + l/2 + APi 50Ut )- 



-1/2 



for a distribution of coherent states with a width of n is 
greater than the best classical fidelity. Indeed the fideli- 
ties of, e.g., F = 0.60 ± 0.02 and F = 0.58 ± 0.02 have 
been obtained for sets of input states with n = 5 and 
n = 20 respectively whereas best classical fidelities for 
these cases are 0.54 and 0.51. 



As sh own in the Supplementa ry Notes to Sherson et c 
(2006c) and in Hammerer et al. (2006) the knowledge of 



AJ^ )OUt , AP^out f° r the teleportation of the coherent 
states corresponds to the complete knowledge of the tele- 
portation map and allows to calculate the fidelity of the 
qubit teleportation as 



F q = (6 + 16s 2 + 24s 4 + 4(# - 1)(1 - 2s 2 ) + 

(£-l) 2 (l-6s 2 ))/6(l + 2s 2 ) 3 , 



(89) 
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FIG. 22 Teleportation results (Sherson et al. 2006c). a) and 
b) Two canonical operators of the verifying pulse plotted as a 
function of the corresponding canonical operators of the input 
pulse for 2 • 10 3 realizations (vacuum units). The slope which 
is close to 1/2 should be multiplied by 2 to account for the 
attenuation of the verifying pulse on the beam splitter. From 
the variances of the distributions along the vertical axis the 
atomic state variances and the fidelity are obtained, c) A set 
of input states with (n) = 5 and random phases. 



where s 2 = 2AA/| out — 1. The quality of mapping 
for the two canonical operators is assumed to be equal 
AX^ ou t = AP^out- Direct demonstration of the qubit 
teleportation under the conditions of |Sherson et al.\ 
(2006c) has not been possible due to the absence of a 
light qubit source wit h the pulse duration of 1 msec. As 
shown theoretically in Sherson et al. (2006c), a qubit fi- 



delity of F q = 0.74 is achievable for n = 1. 

The protocol can in principle be improved by properly 
taking back-action modes, treated here simply as noise 
terms, into account. This is to a large extent a question 
of detector bandw idth and improved post p rocessing of 
photocurrents. In Hammerer et al. ( 2005a| ) it is shown 
that in this way the fidelity can be increased up to 80% 
corresponding to half a unit of added vacuum noise. This 
noise stems from the initial vacuum fluctuations of light 
before light-atom interaction which can in turn be re- 
duced by using a squeezed light for entanglement with 
atoms. One ends up with a similar situation as in the 



standard p rotocol ( |Braunstein and Kimble 1998 Vaid- 



man 



1994|) : the quality of teleportation is in the end lim- 



ited by the amount of available squeezing, which emerges 



here once more to be an irreducible resource (Braunstein 



2005). 



A number of alternative proposals for deterministic 
teleportation inv olving light and atomic ens emble s have 
been sugges ted. Horoshko and Kilin ( |2000[ ); [Mista and| 
Filip (2005) study the application of the state of light 
and atoms created in a QND interaction, cf. (40), as 



a resource for teleportation and show that the use of 
squeezing of light and atoms as well as unbalanced beam 
splitters in the Bell measurement can improve the fidelity. 
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Teleporation of states of light to atoms, based on entan- 
glement between motional degrees of freedom of a Bose- 
Einstein condensate and light has been suggested by |Cola 



et 



(2004); 



Paris et al. 



(2003) 



VII. ERRORS AND FIDELITY FOR DIFFERENT 
INTERFACES 



instance Brask and S0rensen (2008) show that memo- 
ries based on the single pass QND - Faraday interaction 
are not very well suited for the DLCZ-repeater proto- 
col (Duan et al., 2001), because the repeater protocol is 
specifically designed to correct only for the photon loss 
errors. A full evaluation of the performance thus depends 
on the particular application one has in mind. 



Despite impressive successes the unconditional fidelity 
and/or efficiency of the quantum interfaces demonstrated 
so far does not exceed 70 % and in many cases is at the 
level of 20 %. Several factors, some more fundamental 
and some more technical contribute to this. 



Scaling with optical depth. In the theory section we have 
shown that spontaneous emission can be avoided for large 
optical depth d ^> 1. Here we discuss the errors due 
to a finite d for different approaches. In the entangle- 
ment section we have shown that the single ensemble 
squeezing APj[ or the two ensemble correlation Aepr 
is ~ 1/y/d. This limit, however, depends on the exact 
decay mechanism, and a better scaling of single ensem- 
ble squeezing 1/d can be achieved if the QND interaction 
through pha se shi ft measurements using two closed tran- 
sitions (Sec. |II.G| ), is used, since in this case spontaneous 
emission does not add noise to Pa- 

For a quantum memory the inefficiency of storage with 
the beam splitter interaction (Raman or EIT) scales as 
1/d if suitably shaped spatial mode functions are used. 
If storing or reading out the spatially symmetric mode 
is preferable, as in the repeater case where the prob- 
abilistic entanglement is generated in this mode (Sec. 



IV.C ), the scaling is again weak, 1/y/d ( Gorshkov et al. m 
2007a c). However, a much better performan ce can be 



achieved with ensembles in optical cavities (Gorshkov 



et al. 



2007b), where the scaling is 1/Td {T is the cavity 



finesse) , see in particular the e xperiments by Simon et 
fl2007b| ), |Simon et al\ ( |2007c[ ). The protocols based on 
the Faraday interaction naturally couple to the symmet- 
ric modes and may be better suited for this mode than 
the beam splitter interaction. These memory protocols 
have a coup ling constant k 2 ~ 1, w hich yields ~ 1/d or 
~ \og(d)/d ( |Sherson et al. 2006a) error due to sponta- 
neous emission . 

Different protocols are also characterized by different 
kinds of errors. The protocols based on the single pass 
QND - Faraday interaction have variable gain, faithfully 
reproduce mean values of the input states, and thus give 
good fidelity over a large phase space. The error in 
the beam splitter protocols corresponds to a loss on a 
beam splitter and thus gives bad fidelity for states with 
a large amplitude. On the other hand the low fidelity 
beam splitter (Raman and EIT) protocols usually add 
only vacuum noise, whereas single pass QND-based pro- 
tocols add multi-photon errors. Which kind of error is 
less harmful depends on a particular application. For 



Optical losses Losses of photons obviously affect the per- 
formance of memory, and different memories are affected 
to a different extent. Clearly, in the protocols insensi- 
tive to the vacuum component optical losses only lead to 
lower efficiency, that is to a lower probability of success. 
On the contrary, if the goal of a protocol is an uncon- 
ditionally high fidelity, as in the case of protocols based 
on homodyning, optical losses directly affect the fidelity. 
Beside usual losses due to reflection on windows of cells 
and chambers which can be reduced by coating, there are 
losses due to absorption of light by atoms of the mem- 
ory. These losses were included already in our theoretical 
discussion in Sec. |II.C[ where we found that the decoher- 
ence of the quantum fields and the atoms vanished in the 
limit of large d. In addition one should also account for 
the damping of the classical fields. In Sec. |II.C| we also 
discussed that the probability of photon absorption tjl 
is linked to the probability of spontaneous emission of 
an atom t]a by t]l = t]aNa/Nl where Nl,Na are the 
number of photons of the driving field and the number 
of atoms respectively. If Nl ^> Na can be satisfied then 
the photon absorption can be made very small. 



Inhomogeneous broadening Another source of possible er- 
rors is inhomogeneous broadening of the optical transi- 
tions. Solid state systems have a strong inhomogeneous 
broadening because each of the emitters sits in a slightly 
different environment, while inhomogeneous broadening 
in room temperature atomic ensembles is due to the 
Doppler broadening of the atomic lines. In the theoretical 
derivation we only assumed a homogeneous broadening 
of the optical trans itions (the photon echo approach dis- 
cussed in Sec. 



V.E 



is a notable exception). One therefore 
cannot just replace the optical depth d appearing in the 
formulas derived in the theory section by the measured 
optical depth in the presence of inhomogeneous broaden- 
ing. 

A detailed study of the effect of inhomogeneous broad- 
ening for a quantum memory using the beam splitter 
interaction is presented by |Gorshkov et al] ( [200 7d) who 
have shown that even far off-resonance the inhomoge- 
neous broadening of the line still plays a role. The rea- 
son is that the strong field leads to a considerable AC- 
Stark shift of the ground state |1) even far off resonance. 
The inhomogeneous broadening introduces energy shifts 
which are different for each atom and thus causes deco- 
herence of the collective states. As a result the effect 
of inhomogeneous broadening is as severe off resonance 
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as it is on resonance. Hence the total inefficiency of the 
memory protocol has two contributions. The first is the 
spontaneous emission which scales as l/^hom, where <ihom 
is the optical depth of the ensemble without the broad- 
ening. The second contribution is from the inhomoge- 
neous broadening and scales as l/d? nhom , where <ii n hom 
is the actual optical depth in the presence of broaden- 
ing (this scaling assumes that the broadened lines fall 
off sufficiently fast; for other profiles, e.g., Lorentzian, 
the scaling is 1/dinhom)- If one increases the length of 
the sample the system will eventually be dominated by 
the 1/dhom contribution and thus behaves as if it were 
homogeneously broadened. 

The imperfections induced by inhomogeneous broad- 
ening discussed above, can to some degree be reduced 
by taking a more active approach where one tries to en- 
gineer the broadening. In essence the photon echo ap- 



atomic volume and the interaction time is much longer 
than the atomic transient time of flight. This suppression 



proaches discussed in Sec. |V.E| provide examples of such 
engineering of inhomogeneous broadening, where the ex- 
ternally imposed inhomogeneous broadening becomes a 
useful resource. An interesting example of engineering of 



i nhom ogeneous broadening is presented by Afz elius et al. 
(2009), who considers engineering a frequency comb in 
the atomic line shape, e.g. using the hole-burning tech- 
niques discussed in Sec. III.C| Essentially this allows 



for exploiting atoms in a solid state medium which have 
their resonance frequency shifted far away by inhomoge- 
neous broadening (thus reducing the 1/dhom error dis- 
cussed above), while avoiding the detrimental effects of 
inhomogeneous broadening (thus reducing the l/d? nhom 
error) 

For the Faraday interaction the situation is different 
than for the beam splitter interaction. The Faraday in- 
teraction is only employed far off resonance, and unlike 
the case of the beam splitter interaction the AC-stark 
shift does not cause decoherence. The strong classical 
light couples to both ground states and shifts them by 
the same amount. As a result the spin dynamics is 
not affected by a difference in the level shifts, and the 
Faraday interaction becomes insensitive to the inhomo- 
geneous broadening for detunings much larger than the 
hyperfine structure of the excited state and the Doppler 
width. 

In addition to broadening of the optical transitions, in- 
homogeneous broadening also affects the ground states, 
e.g., due to magnetic field gradients. This broadening 
is particularly bad because it leads to decoherence dur- 
ing the period when information is stored in the ground 
states of the ensemble, and is a major limitation for the 
coherence time in many experiments. 



Atomic motion As discussed in Sec. IIIII the atomic mo- 
tion in and out of the optical beam can severely affect 
the performance of the interface. For the Faraday inter- 
action a strong suppression of this effect is achieved in 
the experiments Julsgaard et al. (2001, 2004a); Sherson 



et al. 



is, however, not perfect. As analyzed in detail in Sherson 
(|2006b) the fact that atoms move accross the in- 



teraction volume leads to extra spin noise which should 
be accounted for when the projection noise level is being 
established. 

For experiments with cold atoms the atomic motion 
can still be of a problem, in particular when spin waves 
with very short wavelength are created in the memory. 



Atomic collisions For atomic gasses it is also important 
to consider the effect of atomic collisions. For exper- 
iments with paraffin coated cells atom-atom collisions 
contribu te up to 20 — 40 Hz t o the ground state deco- 
herence (Sherson et al., 2006b). In experiments where a 



buffer gas is used to suppress atomic motion, the alkali 
atoms - buffer gas collisions often have little effect on 
the ground state coherence of the atoms and the mem- 
ory time. The buffer gas does, however, change the dy- 
namics during the interaction with the light. The ho- 
mogeneous broadening due to collisions with the buffer 
gas can be included in the theory by modifying the ho- 
mogenous line width 7 ( Erhard and Helmj|2001[ ). At the 
same time collisions change the velocity of the atoms and 
thereby their Doppler shift. The collisions thus change 
the phase spreading which causes decoherence of the col- 
lective states from being ballistic to diffusive, and this 
reduce the inefficiency from the inhomogeneous broaden- 
ing. 

For the write stage of the probabilistic entanglement 
protocol a different effect related to the collisional broad- 
ening has been observed by Manz et al. (2007) and has 



been discussed theoretically in a different context (Chil- 



dress et al.\ [2005 ) . The collisions with buffer gas atoms 
cause the alkali atoms to emit photons at the resonance 
frequency of the atoms rather than at the anti-Stokes fre- 
quency. To observe the entanglement it is thus necessary 
to filter out these incoherent photons with a frequency 
filter. This effect, however, has little consequences for 
other protocols. 



Geometry of the ensemble. In most cases the two ground 
states of the ensemble are nondegenerate, so there is an 
additional phase factor ex.p[i(k f — ko)z] associated with 
the difference of the k vectors for the classical and quan- 
tum fields. This phase can be absorbed into the defini- 
tion of the mode functions u m (r,t) in (29) and (32), and 



does not play a role for the beam splitter and paramet- 
ric gain interaction applied separately. However, if one 
reads out the memory in the backward direction, which is 
sometimes advantageous, or combines the two protocols, 
the atomic operators should b e redefined and this phase 
can have a detrimental effect ( Andre[ 2005 Duan et al. , 
2002| |Gorshkov et al] [2007c| [SuTmacz et a/.[ |2008| ), see 



et al. (2006c) because the light beams cover most of the also Sec. V.D For the Faraday interaction on the other 
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hand, the mode functions must have a constant phase 
and this means that the atomic ensemble must be much 
smaller than the wavelength corresponding to the ground 
state splitting. In practice this wavelength varies from 
~ 100 m in case of Zeeman splitting with B c± 1 Gauss 
to ~ 3cm in case of hyperfine splitting. In addition the 
Fresnel number corresponding to the shape of the atomic 
ensemble must be large in order to stay within a single 
transverse spatial mode approximation for the Faraday 
interaction. 



Deviation from a two-level ground state model. When the 
ground state level used for the interface is magnetically 
degenerate with more than two states and the detun- 
ing of the strong field is not sufficiently larger than the 
hyperfine splitting of the excited state, the in teraction of 
the Faraday type (a\ vector term in Sec. II. D) is modified 
with the Raman interaction (a,2 tensor term). The Hamil- 



tonian becomes H = Xbs^l^a + Xp^l^a + H.C.. Var- 



ious aspects of this effect have been considered in Juls- 



gaard J2003), Kupriyanov et al. ( 2005 ) , andlSherson et al. 
( 2006b]) , as well as in Geremia et al. (2Q06J). Whereas in 



early work this effect was considered as a source of imper- 
fections, lately it became a subject of intensive studies as 



2009). 



a new resource for interfaces (Wasilewski et al. 

The beam splitter interaction (EIT) experiments 
mostly use two ground hyperfine levels. Magnetic de- 
generacy leads to imperfections which can be reduc ed by 
careful optical pumping and polarization filtering (Choi 



et al. 2008) 



VIM. OUTLOOK 

The light-matter quantum interface, a term coined in 
the end of 1990s, is one of the pillars of the field of Quan- 
tum Information Pr ocessing and Communication (QIPC) 
( Zoller et al. 2005 ) . In less than a decade since the first 



demonstrations of a quantum interface between light and 
an atomic ensemble the ensemble approach has become 
one of the most active areas of research in the field. The 
interactions which seem most promising for interfaces at 
the moment are discussed in this review: the QND - Fara- 
day and Raman interactions, EIT, and photon echo. 

Both fundamental and application driven aspects are 
obvious within this approach. One of the interesting fun- 
damental issues concerns the multi-particle entanglement 
necessarily present in the ensemble-based approach. At 
the same time the interface is a kind of a quantum chan- 
nel, hence its relation to the theory of quantum channel 
capacity should be explored in the future. This issue 
is connected to the fidelity of the interface since it is 
known that a quantum channel with F > 2/ 3 for coher- 
ent states has a non-zero quantum ca pacity ( jGrosshans 
and Grangier, 2001 Wolf et al, 2007). Quantum inter- 



face also allows for storing opti mal quantum clones o f 
a state of light as proposed by Fiurasek et al. (2004). 



Another feature which is intrinsic for the ensembles of 
atoms is their multi-mode capacity. This multimode 
capacity comes both in form a different temporal light 
shapes which are mapped into differ ent longitudin al or 
spectral atomic modes flAfzelius et al. t 2009 Fleischhauer 



and Lukin, 2002 Nunn et al. t 2008 Simon et al. 2007a) 



as well as different transverse modes or "quantum holo- 
grams" ( Surmacz et a l, 2008; Vasilyev et al., 2008), see 
also |Tordrup et al. ( |2008| ). Along the lines of the lat- 
ter the first experiments demonstrating storage of classi- 
cal images via the EIT approach have recently appeared 



(Shuker et a/., 2008 Vudyasetu et a/., 2008) 



Long distance quantum communication is one of the 
most actively pursued applications of the interface at the 
moment. It is based on the combination of probabilistic 
entanglement generation and deterministic entanglement 
swapping - a qua ntum repeater with atomic ensembles 
Duan et al. (2001)- which may serve as the basis for a 



"quantum internet" ((KimbleJ [2008J) . 

New applications for quantum memories and inter- 
faces, such as, e.g., quantum voting and surveying 



( Hillery et al. , 2006 Vaccaro et al.\ |2007[ ) should be ex- 
plored. 

Advanced architectures for quantum computing may 
be enabled by highly efficient photon-based connections 



betwe en small scale atomic processing nodes (|Jiang et al. 



2007). An interesting direction in this respect is the com- 



bination of photon counting and QND-Faraday continu- 
ous variable measurement techniques. It allows to com- 
bine the best of the two worlds: the high efficiency of the 
homodyne measurement and the non-Gaussian states, 
such as Schrodinger cat state s which can be g enerated 
by photon counti ng ( |Genes and Berman [2006; Massar| 
and Polzik 2003). Progress along these lines depends 
critically on the development of highly efficient photon 



counters and photon number re solving detectors ( [Achilles 
et al] 120041 1 Waks et al] 120061). 



For new applications of quantum interfaces a major 
challenge for experimentalists will be to improve the fi- 
delity and efficiency of the interface, and for theorists to 
find protocols where atomic memories with the fidelity 
and efficiency at the level of 90 — 95% - the likely levels 
to be achieved within the next few years - can help to 
achieve goals impossible with classical interfaces. 

Besides the fiducial write- and read-processes and long 
storage times, a quantum interface between light and 
matter will likely have to show yet another key element: 
the possibility to process stored quantum information 
and to allow for quantum logical gate operations for ac- 
tive entanglement purification and error correction (Diir 
and Briegel, 2007). Theoretical studies of the require- 



ments on gate operations in quantum repeater architec- 
tures have been performed in great detail by Briegel et al. 



(1998); Diir et al. (1999); Hartmann et al. (2007), and 



Dorner et al. (2008); Klein et al. (2006) recently invest! 



gated the usage of decoherence free subspaces in quan- 
tum communication. Proposals for ensemble-based im- 
plementations unifying an efficient light matter interface, 
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stable quantum memory and reliable small scale quantum 
processors are rare. A number of theoretical studies sug- 
gest gate operations on stored collective excitations via a 



Rydberg blockade mechanism (Brion et al. 2007 



Lukin 



et al.\ |2001| |Pedersen and M0hner[ |2009| |Petrosyan and 



Fleischhauer[ 2008), via an EIT enhanced optical non- 



linearity (| Andre et a/.||200^|Lukin and Imamoglu, 2001 



Ottaviani et a/.| |2003| |Wang et al.| |2006| ) or"m hybrid 



systems, such as in ( |Rabl et a/.r 2006), where a Cooper 



pair box serves as a saturable, nonlinear element. Initial 
experiments along those lines are in progress. 

Efficient ensemble-based quantum memories and 
matter-light interfaces, small scale quantum processors 
for error correction, repe aters, possibly sa tellite based 
quantum communicat ion (Aspelmeyer et al. 2003 [Pfen- 
nigbauer et al.\ 2005), and even hybrid systems |Ham- 
2009 ) - are the goals of today. The research 



merer et 



performed towards this end is of both fundamental inter- 
est for our understanding of quantum physics and of tech- 
nological importance. Its highly interdisciplinary charac- 
ter encompasses a broad spectrum of fields in physics as 
well as in computer science and information theory. 
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IX. APPENDICES 

A. Adiabatic Elimination 

The effective ground state Hamiltonian can be ob- 
tained by adiabatic elimination. We start with the dipole 
Hamiltonian H[ nt = —E • D. For magnetic sub- levels 
{\g m )} all matrix elements of the ground state dipole op- 
erator vanish {g m \D\g m t) = 0, and similarly for the ex- 
cited state (e m \D\e m f) = 0. Introducing positively and 



negatively oscillating components D = D^ + D^ ) and 
E — E^ + E^ we obtain in the rotating wave approx- 
imation 

H int = -(£(") • #+> + D<"> • £ (+) ), (90) 

where the first (second) term describe down (up) transi- 
tions. Expanding the Hamiltonian we obtain 

H int = -J2 £ (_) • D { +l\g m ){e m ,\ + H.C. (91) 



Using the Hamiltonian Ha + #int we obtain 

-n\gm)(e m '\ = -iA m /|# m )(e m /| 

+ t#+) J2 (&J, m »\9rn){grn"\ " ^Kj<W'><<W|) 

m" 

(92) 

For weak excitation (far below saturation) we can neglect 
the excited state operator |e m ")(e m /|. Next, assuming 
the dynamics is slow compared to the detuning A, we 
ignore the left hand side compared to the first term on the 
right hand side. These approximations are valid provided 
that £ (+) • i3f~L < A m /. In this limit we obtain 



\9m)(e>n 



' 2^ a m,m \9m)(9n 



A m / 



(93) 



To obtain an effective ground state Hamiltonian we 
substitute the expression ( 93 ) back into the Hamiltonian 
Ha + H{ nt • Note that we have to choose normal ordering 
of the operators in Eq. (91) when we insert Eq. ([93]) . A 



detailed discussion of this issue may be found in |Barnett 
and Radmore (1997). Secondly, since the terms in Ha 



involve the excited state population one could be tempted 
to ignore it since the population is proportional to 1/A^. 
On the other hand there is also a factor A m in front of 
this term, so that in total this is on the order of 1/A m 
and we have to include it. To treat this term we introduce 
any intermediate state \go) 



>){e m >\ = |e m /)(0 O | • \9o)(en 



(94) 



and substitute the result ( |93| ) for \go)(e rn '\. By doing this 
we arrive at the effective Hamiltonian in Eq. (25). 



B. Three-dimensional Hamiltonians 



In the main text we present the general expression 



(27) for the full three-dimensional Hamiltonian. In this 



appendix we discuss the 3D Hamiltonians for the three 
model system, beams splitter, parametric gain, and Fara- 
day interaction. 

First we consider the atomic beam splitter interaction 
inEfa). We assume that the state 1 1) is coupled to an 
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excited state by a classical field, so that E(r) = (E(r)). 
Inserting the expression for the quantized electric field in 
Eq. (12) we obtain 



H BS = / d 3 r 



-Mr,*) I s 
4A 



a ] A {r)a A {r) 



\m\ 



^— J2 \u m (rj_; z)\ 2 a} Lm {z)a Ltm {z) 



■E 



^^C(^) 

xa[(z)a A (r)+H.C. 



(95) 



where the coupling constant g(f) and resonant Rabi fre- 
quency are defined as in Eq. (31) except with z replaced 
by r. The first term in this Hamiltonian is the AC- Stark 
of the ground state caused by the classical field, the sec- 
ond term is the change in the index of refraction that 
the quantum field experiences, and the last term repre- 
sents the exchange of excitation between the light and 
the ground state coherence. We have here ignored the 
small AC-Stark shift caused by the weak quantum field. 
In case of the parametric gain interaction shown in 
Fig. [2Vb), the fields interact to the transitions which are 
flipped compared to the case of the beam splitter interac- 
tion. The emission of a photon is therefore coupled to an 
atomic transition in the opposite direction and we obtain 
the parametric Hamiltonian by making the replacement 
o>a{t) *+ &a(?) m the coupling between atoms and light. 
In addition, the effective AC-Stark shift has the opposite 
sign. The Hamiltonian thus reads 



H G 



-I 



d 3 r 



£ 



ift(r,*)i 2 



u* 



a\{r)a A {r) 



g*(fMr,t) 



2A 



u* m (r±;z) 



(96) 



at ro (*)4(0+H.C. 



There is no index of refraction here since now the quan- 
tum field couples to the almost empty transition. On the 
other hand the classical field now sees the atomic popula- 
tion, and we should include the index of refraction in the 
propagation of the classical light, which enters through 
the phase of Q. 

For the Faraday interaction we assume a strong x- 
polarized classical field, and consider the quantum field 
in the ^-polarization as shown in Fig. [2|c). There are 
now two paths which lead to creation of a photon: pho- 
ton creation (or absorption) can appear both through 
the creation and annihilation of an atomic excitation. 
As discussed in the main text the Faraday interaction 
is a combination of the two Hamiltonians in Eqs. (95) 



and (96) with equal weights. For the 1/2-1/2 transition 
depicted in Fig. [2| however, the two paths leading to 
the creation of a photon have opposite signs due to the 
Clebsch-Gordon coefficients. The coupling constant g is 
defined for circularly polarized light, so we should in- 
clude a factor of y/2 coming from the expansion of the y— 
polarized quantum field in the circular polarization basis. 
The Hamiltonian is then given by Hp = (Hbs—Hg)/V2- 
For the spin 1/2 system classical and quantum fields 
experience the same index of refraction, so we can remove 
the index of refraction from Hbs if we a ls° ignore it for 
the classical field. An important result for this Hamilto- 
nian is that the AC-Stark shift of the ground state dis- 
appears, because the classical field shifts the two atomic 
states by the same amount. 



C. Propagation equations for light 

The light Hamiltonian (for a single transverse mode) 
is given by 



H L = Y]\k\cala k , 



(97) 



where the sum is over all the longitudinal wave vectors. 
We now derive the time evolution of the operator cll(z) 
defined in Eq. (10). If we assume that the modes con- 



tributing to the slowly varying operator cll(z) are cen- 
tered around a large positive value fco, we can ignore the 
absolute value of the wave vector. Combining the time 
derivative caused by the explicit time dependence intro- 
duced in (10) with the time evolution caused by Hl we 
find the Heisenberg equation of motion 

— — = iuJoa L -i[a L ,(H L + H int )] = -c—a L -i[a L , H int ]. 

(98) 
To simplify the equations we introduce a rescaled time 
variable by defining new operators cll(z,t) = ap{z^t = 
r + z/c) and cla(z,t) = aA(z,t = r + z/c). The prop- 
agation equation for the light can then be simplified by 
using 



d_ 

dt 



d\ d 

c— )a L (z,t = r + z/c) = c—a L (z,r). (99) 



dz 



dz 



We also take the classical field to be moving in the pos- 
itive z-direction, in which case the Rabi- frequency is 
Q(z, t) = Q(t — z/c) = Q(t). The equations of motion in 
the main text always use this rescaled time and we omit 
the tilde on the operators and simply denote the rescaled 
time t by t. 

For the parametric gain the z dependence of the 
Rabi- frequency does not disappear completely from 
the equations of motion since we should include the 
ex.p(i j dz f \g(z f )\ 2 /A) dependence associated with the 
index of refraction for the classical field. For the Faraday 
interaction the same factor appears in the propagation of 
the quantum field and both factors can be omitted. 
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D. Inclusion of spontaneous emission 



In Sec. II. B we derived equations of motion ignoring 
spontaneous emission. To include spontaneous emission 
we should modify Eq. (92) such that it contains a de- 



cay as well as the Langevin noise operators associated 
with the decay. The general treatment of the Langevin 
noise operators is quite complicated because they depend 
on details of the decay mechanism. In particular, the 
equations for the parametric gain and Faraday interac- 
tion depend on whether the atoms end up in the states 
|0) and |1) or some auxiliary states \a m ) (Fig.|2|. For the 
beam splitter interaction, on the other hand, the precise 
state that the atoms decay to is less important ( |Gorshkov 
et al 2QQ7b|c ). Here we will for simplicity only consider 
the decay to some auxiliary states \a m ). 

To describe spontaneous emission we assume that each 
of the excited states \e m )j of the jth atom couples to 
some state \a m )j via a continuum of modes described 
by the annihilation operators bj(uu). We assume that 
each atom couples to its own continuum. By doing 
so we ignore collective scattering effects such as su- 
perradiance and Bragg scattering, and we also ignore 
dipole-dipole interactions mediated by other than for- 
ward modes. Whether this is a suitable approximation 
should be evaluated for each particular realization, but 
to our knowledge little work has been done on this sub- 
ject. The decay of the jth atom may then be described 
as 



decay 



duo 



6t (*)&>) 



+ p(u) ^2 (9m(w)\em)j(a>m\bj(w) + H.C.) 

m 

(100) 

where g m (uj) is the coupling constant of the mth excited 
state, and p(uS) is the density of states of the continuum. 
To arrive at the equations of motion for the atomic op- 
erators we first derive t he eq uations of motion for bj(u) 
using the Hamiltonian ( |1QQ[ ). We then formally solve 
this equation by integrating over time and substitute the 
result into the equation for |g m ) 7 (e m |. In the Markov 



approximation (Barnett and Radmore, 1997) Eq. (92) 



acquires the additional terms 



d 



lm/j \^m' 



(101) 

where we have ignored the Lamb shift. The decay and 
noise operators are given by 

7 m = 27r\g rn (uj rn )\ 2 p(uj m ) (102) 



and 



i\g m )j{a m >\ 



f j ( f \ _ «\ym/ 3 

am, em' V / 



/ dujg m >{uj)bj{uj,t = 0)e 



-i(w-w m /)t 



(103) 



with uo m being the transition frequency for the |e m ) to 
\a m ) transition. Note that if there is a difference between 
the population decay rate and twice the decay rate of the 
polarization, e.g., due to collisional broadening, the de- 
cay rate 7 m /2 appearing here should be the polarization 
decay rate. The correlation functions within the Markov 
approximation are 



(F j] F j> ) 

\ gm,em' gm" ,em"' I 



\ gm,em 



>F, 



tf't 



gm" ,em 



»,> 



o 

S(t - t')5 m ^ m >»5jj>{\g m )(g m »\). 

(104) 



Since most of the atoms are always in the state |0) the 
only important noise operators are F g o, e m- 

""^are 



The equations of motion that we derive in Sec. II. B 



due to the time evolution of az, and the time evolution 
of j+(r) used in the definition of a a in Eq. (21). Be- 
cause the decay simply adds a term in Eq. ( 101 ) which 
is similar to the detuning term in Eq. (92) we can ob- 



tain the equations of motion for these operators by using 
the replacement in Eq. (42). Note that for operators 



\gm){^m'\ the detuning and decay appears in the combi- 
nation A — ry/2 so that, e.g., for the ground state op- 
erators \gm)(gm'\ coupling to |# m )(e m "| one should use 
the minus sign in the substitution. On the other hand 
if the coupling is to operators \e m fr)(g m /\ one should use 
the plus sign. In order to conserve the commutation rela- 
tion for the operators a a and a^ one should also include 
the Langevin noise operators and the time derivative of 
(jx(r)) used in the definition of a^ (21). 

For the beam splitter interaction the time derivative of 
{jxif)) can be ignored because the strong classical field 
couples to an almost empty transition. The equation 
of motion can therefore be obtained simply by making 
the substitution (42). The equation of motion for the 



light (35) arises from the commutator of cll with the 



Hamiltonian ( |90| ). In the resulting equation of motion a^ 
couples to |0)7etand we should therefore use the minus 
sign in Eq. ( [42] ) . The atomic annihilation operator is 
proportional to|0)(l|, which couples to |0)(e| and |e)(l|. 
The contribution from |e)(l|, is however very small and 
has been neglected in Eq. ( 35 ) , and we should once again 
use the minus sign in the substitution (42). Ignoring the 
noise operators we arrive at the equations of motion (43). 



Including the noise operators gives the additional terms 



dz 
d_ 

dt 



a L (z,t)=... + -^^F(z,t) 



a A (z,t) 



vW 



(105) 



-F(z,t). 



2A-Z7 

The noise F(z,t) appearing here is defined by first defin- 
ing 



F(r,t) 



1 



y/n{?,t) 



Y, F ke{mr-r 3 ). 



(106) 



Integrating the equations of motion over the transverse 
coordinates we define F(z,t) in analogy with Eq. (29) 
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(ignoring the mode index m for a single mode). The 
resulting operator has the standard expectation value of 
vacuum noise (F^ '(z,t)F(z' \t')) = 0, (F(z, t)F\z' , t')) = 
S(z — z')5(t — t'). The addition of these noise operators 
ensures that the atom and light operators retains the 
correct commutation relations. 

For the parametric gain the situation is a little differ- 
ent. The light operator cll couples to the coherence |l)(e| 
and we should still use the substitution (42) with the mi- 
nus sign. The atomic coherence |0)(1| again couples to 
|0)(e| and |e)(l|, but now we can no longer neglect |e)(l|, 
which gives rise to the AC-Stark shift [the first term in 
the spin equation in Eq. (36)]. We should therefore use 
the plus sign in the substitution (42) for this term. 



To include the change of the mean spin (j x (r,t)} we 
include the decay of the densities no and n\ of atoms in 
states |0) and |1) given by 



d 
dt 
d 
dt 



n (r,t) 
ni(r,t)«0 



j\mt)i 

4A 2 + 7 2 



■n (r,t). 



(107) 



Since we also assume that the decay takes the atoms 
to some auxiliary state |a m ), an initially fully polar- 
ized state will remain fully polarized so that (j x (r,t)} « 
no(r,t)/2 as long as the interaction with the quantum 
field is weak (note, however, that we are describing a 
situation leading to superradiant scattering so that this 
approximation may break down very quickly). The con- 
tribution from the time derivative of the mean spin then 
exactly cancels the decay of a A (z, t) arising from the sub- 
stitution (42 ). Since there is no longer any decay, one also 



finds that the resulting equations of motion (53) do not 
contain any noise. 

Finally for the Faraday interaction the easiest way to 
proceed is again to combine the results for the beam split- 
ter interaction and the parametric gain. Similarly to the 
discussion in Appendix |IX.B the equations of motion for 
the Faraday interaction can therefore be obtained by sub- 
tracting the right hand side of Eq. (53) from the right 
hand side of Eq. (43) and dividing by \/2. Furthermore, 



we again ignore the change of the propagation caused by 
the index of refraction because this is accompanied by a 
similar change in the propagation of the classical light. 
The resulting equations in terms of the x and p operators 
are given in Eq. (54) without the noise operators. The 



noise operators give the additional terms 



d_ 
dz 

d_ 
dz 

d_ 

dt 

d_ 

dt 



x L (z,t) =... 
p L (z,t) =... 
x A (z,t) =... 
p A (z,t) =... 



gV^f 



v/4A 2 + 7 2 

gV^f 



V4A 2 + 7 2 
V4A 2 + 7 2 
V / 4A 2 + 7 2 



F x (z,t) 
F p (z,t) 
F x {z,t) 
F p (z,t). 



(108) 



Here we have for simplicity assumed Q and g to be real 
and have defined new noise operators 



F x (z,t)=^=(e i +F(z,t)+e- i +Ft(z,t)) 
F p (z, t) =-JL (e^F(z, t) - e-'W (z, t)) , 



(109) 



which have the standard commutation relation 
[F x (z, t), F p (z>, t')\ = i5{z - z')5(t - t'). 

E. Dimensionless equations of motion 

Casting the equations of motion for all three basic 
interactions in a dimensionless form allows to see that 



the constant n (56), naturally plays the role of the cou- 
pling constant for all protocols. We introduce the di- 
mensionless position and time coordinates s = z/L and 
v = /i(0, t)/h(0,T) running from to 1. The rescaled 
time variable simplifies the equations because it is pro- 
portional to the total integrated intensity of the field. 
In the weak saturation limit the dynamics is completely 
controlled by the incident number of photons in the clas- 
sical field. Changing the intensity will thus influence the 
temporal dynamics of the system, but the final state pri- 
marily depends on the total number of incident photons. 
When using such rescaled coordinates it is desirable to 
also change the field operators such that, e.g., an in- 
cident light field operator which is normalized in time 
[aL{t),ar L (t f )\ = S(t — t') is now normalized relative to 
the new time variable [cll(v), cl l (v')} = S(v — v'). This 
normalization is achieved with the rescaling 



cla(s = z/L) = VLcia(z) 



a L (v(t)) = V / ^ 0/ ,/ a L (t), 



(110) 



djft(t) 



where we have used the dimensionless coupling constant 
k = ^/ft(0, T), which appeared in the solution for the 
Faraday interaction, c.f., Eq. (57) in the far off-resonant 



limit A ^> 7. Below we omit the tilde on the operators 
on the left hand side. In these new rescaled variables the 



equations of motion for the beam splitter interaction (43) 
become 

— a L (s,v) = i 2 (2A-i ) aL ( s ' v " j ~ i ^^ a A{s,v) 



-a A (s, ,) = *.— 



2. e i0 

2-a A (s,v) -w—a L (s,v). 

(Ill) 



With the same rescaling for the parametric gain inter- 
action (53) we find 



d e i(f) + 

—a L (s,v) = -iK—a\(s,v) 

d VA . j* f 

— -a A (s, v) = -i—j-a A {s, v) - iK—a ] L (s, v), 



(112) 
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where we have again assumed that g is real. 

For the Faraday interaction the equations of motion 



(54) become 

d 7 2 d 

-Q-^L^V) = K PA {S,V) - 4A 2 +7 2 2 XL(g ^ ) 



d_ 
ds 
d_ 
dv 

dv 



Pl(s,v) 



7 



d 



4A 2 + 7 2 2 



Pl(s,v) 



(113) 



X A (s,v) = Kp L (s,v) - — X A (s,v) 



2d 



Pa(s,v) 



2d 



Pa(s,v), 



where we have ignored small corrections which vanish in 
the far detuned limit. 



F. Tensor Decomposition 

The coefficients determining the strength of the irre- 
ducible tensor components are 



a k (A) = (-) 1+F c k (2k + l) 



2F + 1 



Sr 



L F' 



5 F ,/A 



(2^+i){^;{H^T} 



where the expressions in curly brackets are 6j-symbols, 

A = Ap+i, $f' = A F+ i — Ap> and 



co = 1, 



c 2 



Cl 




y/WF(F + l)(2F - 1)(2F + 3) ' 



The last line is valid for F > 1/2, that is nuclear spin 
/ t^ 0, and has to be replaced by c 2 = for I = 0. 

In the asymptotic limit of large (blue) detuning, —A ^> 
5f', the sum in square brackets can be simplified by 
means of 

E (-) F '(2F' + l){^'{} 2 {^} = 



-F-l 



(_\{2J+2F+J'+I+k) ( J I F\ f J J k \ 
V ) X Fk J J X 11 J' ) 



to get 



a k = lim a fc (A) = (-) 2J+F+J ' +/+fc+1 c fc (2fc + l) 

A— )► — OO 



2F + 1 



( J I F\ ( J J k \ 
X Fk J ] X 1 1 J' S • 



From this expression it is evident that a 2 has to vanish 
because the triple {J, J, /c} = {1/2, 1/2, 2} does not sat- 
isfy the triangle inequality. For the particular case of the 
Cesium (I = 7/2) £> 2 -line at F = 4 -> F' = 3,4, 5 the 
asymptotic values of the non-vanishing coefficients are 
ao = 1/6 and a± = 1/24. 
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